Transient Simulations Unsaturated Zone
G. Bier
April, 2025
Introduction
Markdown to set up an unsaturated zone model and to analyze the impact of the van Genuchten parameters on Evapotranspiration, moisture distribution and groundwater recharge.
Cleaning the memory and loading required package(s)
rm(list = ls())
library(FVFE1D)
library(MIPLOT)Loading the soil data (Staring Series)
load(file='Staring.Rdata')
print(paste("Soil type: ",names(soil.set)," ",lapply(soil.set, function(x) x$name))) [1] "Soil type: B1 Non-loamy sand"
[2] "Soil type: B2 Loamy sand"
[3] "Soil type: B3 Very loamy sand"
[4] "Soil type: B4 Extremely loamy sand"
[5] "Soil type: B5 Coarse sand"
[6] "Soil type: B6 Boulder clay"
[7] "Soil type: B7 Sandy loam"
[8] "Soil type: B8 Silt loam"
[9] "Soil type: B9 Clayey loam"
[10] "Soil type: B10 Light clay"
[11] "Soil type: B11 Heavy clay"
[12] "Soil type: B12 Very heavy clay"
[13] "Soil type: B13 Loam"
[14] "Soil type: B14 Heavy loam"
[15] "Soil type: B15 Peaty sand"
[16] "Soil type: B16 Sandy peat"
[17] "Soil type: B17 Peaty clay"
[18] "Soil type: B18 Clayey peat"
[19] "Soil type: O1 Non-loamy sand"
[20] "Soil type: O2 Loamy sand"
[21] "Soil type: O3 Very loamy sand"
[22] "Soil type: O4 Extremely loamy sand"
[23] "Soil type: O5 Coarse sand"
[24] "Soil type: O6 Boulder clay"
[25] "Soil type: O7 River loam"
[26] "Soil type: O8 Sandy loam"
[27] "Soil type: O9 Silt loam"
[28] "Soil type: O10 Clayey loam"
[29] "Soil type: O11 Light clay"
[30] "Soil type: O12 Heavy clay"
[31] "Soil type: O13 Very heavy clay"
[32] "Soil type: O14 Loam"
[33] "Soil type: O15 Heavy loam"
[34] "Soil type: O16 Oligotrophic peat"
[35] "Soil type: O17 Eutrophic peat"
[36] "Soil type: O18 Peaty layer" print(paste("Ksat values: ",lapply(soil.set, function(y) y$name)," ",lapply(soil.set, function(x) x$ksat))) [1] "Ksat values: Non-loamy sand 23.41"
[2] "Ksat values: Loamy sand 12.52"
[3] "Ksat values: Very loamy sand 15.42"
[4] "Ksat values: Extremely loamy sand 29.22"
[5] "Ksat values: Coarse sand 52.91"
[6] "Ksat values: Boulder clay 100.69"
[7] "Ksat values: Sandy loam 14.07"
[8] "Ksat values: Silt loam 2.36"
[9] "Ksat values: Clayey loam 1.54"
[10] "Ksat values: Light clay 0.7"
[11] "Ksat values: Heavy clay 4.53"
[12] "Ksat values: Very heavy clay 5.37"
[13] "Ksat values: Loam 12.98"
[14] "Ksat values: Heavy loam 0.8"
[15] "Ksat values: Peaty sand 81.28"
[16] "Ksat values: Sandy peat 6.79"
[17] "Ksat values: Peaty clay 4.46"
[18] "Ksat values: Clayey peat 6.67"
[19] "Ksat values: Non-loamy sand 15.22"
[20] "Ksat values: Loamy sand 12.68"
[21] "Ksat values: Very loamy sand 10.87"
[22] "Ksat values: Extremely loamy sand 9.86"
[23] "Ksat values: Coarse sand 25"
[24] "Ksat values: Boulder clay 33.92"
[25] "Ksat values: River loam 39.1"
[26] "Ksat values: Sandy loam 9.08"
[27] "Ksat values: Silt loam 2.23"
[28] "Ksat values: Clayey loam 2.12"
[29] "Ksat values: Light clay 13.79"
[30] "Ksat values: Heavy clay 0.5"
[31] "Ksat values: Very heavy clay 4.37"
[32] "Ksat values: Loam 1.51"
[33] "Ksat values: Heavy loam 3.7"
[34] "Ksat values: Oligotrophic peat 1.07"
[35] "Ksat values: Eutrophic peat 2.93"
[36] "Ksat values: Peaty layer 43.45" Model setup
Internal flux in the unsaturated zone
The equation describing the internal flux in the unsaturated zone strongly resembles the basic Darcy equation to simulate flow in the saturated zone. Most important differences are the hydraulic conductivity and the state. Now the state is \(\psi\) instead of \(H\). Both are related: \(H = \psi + z\) where the \(z\) is the datum/reference (in the Netherlands: “NAP”).
\[ q=-k(\psi)\left [\frac{\partial \psi}{\partial z} + 1 \right] \]
Building the model for the Holterbeg
Be aware the length unit is in cm!
The model is based on:
- 1D vertical flow
- Richards/Darcy’s equation
- Top boundary flux based precipitation
- Bottom boundary flux based on Dirichlet or Cauchy type
- Spatial flux for the rooting function for transpiration
- Soil evaporation is neglected
Choose a soil type from the Staring series
From the Staring Series you can choose a typical soil type for ice pushed ridges. For example the soil type “Coarse sand” B5.
soiltype = soil.set$B1
print(unlist(soiltype))$soiltype
[1] "B1"
$theta.res
[1] 0.02
$theta.sat
[1] 0.43
$n
[1] 1.801
$alpha
[1] 0.0234
$ksat
[1] 23.41
$lambda
[1] 0
$name
[1] "Non-loamy sand"
$k.fun
function (v)
.approxfun(x, y, v, method, yleft, yright, f)
<environment: 0x00000200efd9b670>
$theta.fun
function (v)
.approxfun(x, y, v, method, yleft, yright, f)
<environment: 0x00000200efd9be88>
$c.fun
function (v)
.approxfun(x, y, v, method, yleft, yright, f)
<environment: 0x00000200efa32508>
$Richards.flux
function(x,state,gradstate)
{
return(-soil.set$B1$k.fun(state)*(gradstate+1))
}## domain of the model
top = 0.0 # at the surface
bottom = -200 # so 200 cm below surface!!
log_top = 0.1
log_bottom = log10(-bottom)
log_z = seq(from = log_top, to = log_bottom,by = 0.01)
nodes = -10^log_z
nodes = c(0.0,nodes)
domain=c(bottom,top)
###description of the internal flux##############
unsatdarcy.flux = function(z,psi,dpsidz)
{
return(-soiltype$k.fun(psi)*(dpsidz+1))
}
################### make the model #################
UnsatStat = newFLOW1D(domain,unsatdarcy.flux,"Unsaturated zone, 1D stationary")
####### discretisation #####################
# position the nodes 2 cm apart
#nodes = seq(from=domain[1],to=domain[2],by=1)
#set.discretisation(UnsatStat,nodes,"FV") #setting to FEquartic slows transient simulations way too much
set.discretisation(UnsatStat,nodes,"FEquadratic")
# An initialisation is certainly important (as the default initialisation by zeros does not make sense here).
# So we initialize with a no-flux profile:
Psi.equi = approxfun(c(domain[1],domain[2]),c(domain[2],domain[1]))
do.initialize(UnsatStat,Psi.equi)
# We want to avoid pressures to be too negative (<-1E-5)
Psi.in.range = function(z,Psi) # it has 2 arguments; the position and the state. Set the min. value to pF4.2 not pF5
{
# if(Psi<=-1*10^4.2) return(FALSE)
if(Psi <= -16000) return(FALSE) # also in the Feddes function -16000 which is 4.20412 and not 4.2
#if(Psi<=-1*10^5) return(FALSE)
else return(TRUE)
}
set.isacceptable(UnsatStat,Psi.in.range) #the set.isacceptable function limits the possible values that the state can become
############## Boundary conditions ##################
set.BC.fixedstate(UnsatStat,"left",0.0)
# At the bottom of the model, where the groundwater level is located, both pressure head and elevation head are zero.
ppt = 0.1
set.BC.fixedflux(UnsatStat,"right","ppt")
# 0.1cm/d=1mm/d. Although this flow is directed downwards, we define it positive because it goes into the domain.
######## checking and solving the model ####################
summary(UnsatStat)one dimensional flow model
name: Unsaturated zone, 1D stationary
model has no spatial fluxes
model has no point fluxes
left BC is of type state
right BC is of type given flux
numerical model of type FEquadratic1D
number of nodes is 222 solve.steps(UnsatStat,verboselevel=1)[1] "iteration 1 ; RMSM= 0.100091343973599 ; MAM= 1.48794139671791"
[1] "iteration 2 ; RMSM= 0.00323131833384584 ; MAM= 0.0303723562135154"
[1] "iteration 3 ; RMSM= 0.000778347115860355 ; MAM= 0.00629445455593794"
[1] "iteration 4 ; RMSM= 0.0000813461389181065 ; MAM= 0.000853309088570731"
[1] "stopped because of small RMSM= 0.00000129427348443462"
$RMSM
[1] 0.000001294273
$MAM
[1] 0.00001550802plot(UnsatStat,fluxplot = T,vertical = T)
Additional flow functions
In addition there will be functions for ponding, the lower boundary conditions (simulating the drainage effect) and transpiration of roots (Feddes function). The lower boundary condition is a Cauchy type BC: \[ Q_{drainge}= \frac{H_{drn.level} - H_{bottom}}{C} \]
drn.resistance = 50 #unit is in days
drn.level = -120 # so XX cm below surface level
drn.fun = function(psi)
{
# if((drn.level - (psi + bottom)) < 0) return((drn.level - (psi + bottom))/drn.resistance)
# else return(0)
return((drn.level - (psi + bottom))/drn.resistance)
}
Replacing the bottom fixed head boundary condition with the drainage type
set.BC.fluxstate(UnsatStat,"left", drn.fun)
summary(UnsatStat)one dimensional flow model
name: Unsaturated zone, 1D stationary
model has no spatial fluxes
model has no point fluxes
left BC is of type state-flux
right BC is of type given flux
numerical model of type FEquadratic1D
number of nodes is 222 drn.psi.equi = approxfun(c(domain[1],domain[2]),c((drn.level - bottom),drn.level))
do.initialize(UnsatStat,drn.psi.equi)
solve.steps(UnsatStat,verboselevel = 1)[1] "iteration 1 ; RMSM= 0.00454630834988694 ; MAM= 0.0673274601252841"
[1] "iteration 2 ; RMSM= 0.00369561174124103 ; MAM= 0.0498855800816715"
[1] "iteration 3 ; RMSM= 0.000419498101987421 ; MAM= 0.00398168606717239"
[1] "iteration 4 ; RMSM= 0.0000107526284007399 ; MAM= 0.000110670479414737"
[1] "stopped because of small RMSM= 0.000000000445958140107072"
$RMSM
[1] 0.0000000004459581
$MAM
[1] 0.000000004724472plot(UnsatStat,fluxplot = T,vertical = T)
##trying the new bottom boundary condition with a view different ppt's
ppt = 0.2
solve.steps(UnsatStat,verboselevel = 1)[1] "iteration 1 ; RMSM= 0.00502938503088121 ; MAM= 0.0606056170041377"
[1] "iteration 2 ; RMSM= 0.000568647268186233 ; MAM= 0.00581109257855628"
[1] "iteration 3 ; RMSM= 0.0000161578697832701 ; MAM= 0.000142062648655544"
[1] "stopped because of small RMSM= 0.000000000444692678300773"
$RMSM
[1] 0.0000000004446927
$MAM
[1] 0.000000004557229plot(UnsatStat,fluxplot = T,vertical = T)
#ppt = -0.2
# solve.steps(UnsatStat,verboselevel = 1)
# plot(UnsatStat,fluxplot = T,vertical = T)
# dataframe.balance(UnsatStat)
# dataframe.boundaries(UnsatStat)
###plotting functions showing the water table in the profile,
### the moisture distribution and transpiration through the roots
plot.wtable = function(model)
{
df = dataframe.states(model)
plot(df$state, df$x, type = "o",col = "blue",lwd=2.5,
main = "Psi(z), location water table", xlab = "psi in cm", ylab = "depth profile")
abline(h = drn.level, col = "brown", lwd= 3, lty = "dashed")
text(20,drn.level, "drainage level", col = "brown")
psi.fun = approxfun(df$state,df$x)
abline(h = psi.fun(0),col = "red", lwd = 2)
text(20,psi.fun(0),"water table" , col = "red")
grid()
}
plot.theta = function(model)
{
df = dataframe.states(model)
theta.psi = soiltype$theta.fun(df$state)
plot(theta.psi,df$x, type = "o", col = "green", lwd = 2.5,
main = "Theta(z), moisture content in profile", xlab = "theta", ylab = "depth profile")
abline(h = drn.level, col = "brown", lwd= 3, lty = "dashed")
text((mean(theta.psi)),drn.level, "drainage level", col = "brown")
grid()
}
plot.wtable(UnsatStat)
plot.theta(UnsatStat)
Adding the vegetation to the model
Evapotranspiration is an important environmental condition for processes near the soil surface, such as the development of vegetation. Evapotranspiration depends on:
- The potential evapotranspiration rate, which is determined by
weather conditions;
- The rooting distribution of the vegetation;
- The actual soil moisture contentin the rooting zone.
The total potential evapotranspiration rate of the vegetation is determined by weather conditions and will be given.
The rooting distribution may be described by a function
root.distr.fun(z) which gives the fraction of roots present
at each depth (giving 1 when integrated with depth). The rooting depth
for each vegetation type can be described in many ways. Here, at this
stage we will assume grass with 75% rooting in the first 20 cm and 25%
till 110 cm below surface. (Just my guess).
root.partition = c(0.0,-20.0,-20.01,-110.0,-110.01)
root.rates.fractions = c(0.75/20,0.75/20,0.25/90,0.25/90,0.0)
root.distr.fun = approxfun(root.partition,root.rates.fractions,rule = 2)
plot(root.partition,root.distr.fun(root.partition), type = "o",
col = "green", lwd = 2, main = "Root distribution function",
ylab = "Fraction of the root distribution (-)", xlab = "Rooting depth (cm) rel. to surface level")
grid()
###check for the integral
20*0.75/20 + 90*0.25/90[1] 1z.root = seq(-109.5,-0.5,by = 1)
sum(root.distr.fun(z.root))[1] 1The vegetation can only transpire optimally when there is sufficient moisture in the soil which is ‘easily’ available for the vegetation. The effect of soil moisture can be simulated with the so called Feddes (retired Prof. from WUR) function.
The code chunk below illustrates such a Feddes function. This function relates the root water uptake at a certain depth with the pressure (suction) at that depth. When the soil is too wet or too dry, plant transpiration and thus root uptake are reduced below their potential rates, and the Feddes reduction factor is smaller than 1.
options(scipen = 999) #avoid scientific notation (with E power) in the plot
#Feddes.psiabs = c(0.1,100,1000,16000) #typical psi values to define Feddes function. Note these are absolute psi's due to the plot. NOT in the function
Feddes.psiabs = c(0.1,100,1000,15000)
Feddes.red = c(0,1,1,0)
Feddes.fun= approxfun(-Feddes.psiabs,Feddes.red,rule=2)
plot(Feddes.psiabs,Feddes.fun(-Feddes.psiabs),type="l",lwd=3,xlim=rev(range(Feddes.psiabs)),log="x",ylab="Root uptake reduction",xlab="Suction (|psi|)")
abline(v=Feddes.psiabs,col="red")text(1,0.9,"Too wet", col="red")
text(300,0.7,"Optimal", col="red")text(5000,0.9,"Too dry",col="red")
grid()
Adding this spatial flux to the model
The effect of the vegetation, toots extracting water from the soil for transpiration depending on the root distribution, suction (\(\psi\)) and the potential evapotranpsiration.
Three aspects are considered:
- Rooting distribution
root.distr.fun
- Feddes function
Feddes.fun
- Potential evapotranspiration
Epot
The function to calculate this is :
Eact = root.distr.fun(z)*Feddes.fun(psi)*Epot
Eact = function(z,psi)
{
return(-root.distr.fun(z) * Feddes.fun(psi) * Epot)
}
add.spatialflux(UnsatStat, rate = Eact, name = "Transpiration")
current model
The current model consists of:
- precipitation applied to the upper boundary as a prescribed flux
ppt
- potential evapotranspiration
epotas a flux rate to determine Eact
- actual evapotranspiration
Eactas a spatial flux (along the roots)
- Cauchy type boundary condition at the lower boundary simulating the effect of a local drainage system
plot.roots = function(model)
{
df = dataframe.externalfluxes(model)
plot(df$`S{Transpiration}dx`,df$x, type = "o", col = "green", lwd = 2.5,
main = "Moisture extraction by roots in depth", xlab = "rate (cm/d)", ylab = "depth profile")
abline(h = drn.level, col = "brown", lwd = 3, lty = "dashed")
text(mean(df$`S{Transpiration}dx`),drn.level,"drainage level", col = "brown")
grid()
}
summary(UnsatStat)one dimensional flow model
name: Unsaturated zone, 1D stationary
model has 1 spatial flux, with name
Transpiration
model has no point fluxes
left BC is of type state-flux
right BC is of type given flux
numerical model of type FEquadratic1D
number of nodes is 222 ####testing the model with meteo data
ppt = 0.124
Epot = 0.05#0.2
do.initialize(UnsatStat, drn.psi.equi)
solve.steps(UnsatStat,verboselevel = 1)[1] "iteration 1 ; RMSM= 0.00562357446327675 ; MAM= 0.0834332892218055"
[1] "iteration 2 ; RMSM= 0.00412023136299819 ; MAM= 0.0591990193063127"
[1] "iteration 3 ; RMSM= 0.000360955447892504 ; MAM= 0.00471218175091925"
[1] "iteration 4 ; RMSM= 0.0000140871667402979 ; MAM= 0.000198318505135436"
[1] "stopped because of small RMSM= 0.000000000621385426734169"
$RMSM
[1] 0.0000000006213854
$MAM
[1] 0.000000008903264plot(UnsatStat,fluxplot = T)
plot.wtable(UnsatStat)
plot.theta(UnsatStat)
plot.roots(UnsatStat)
bal.asfixedflux = dataframe.balance(UnsatStat)
exter.asfixedflux = dataframe.externalfluxes(UnsatStat)In case of poding
When the amount of precipitation is larger than the soil can absorb, ponding (nl: plasvorming) will take place To simulate this, an additional Cauchy type boundary condition will be applied to the upper node.
\[ q_{ponding} = \frac{\psi_{surface}-0}{C_{top}} \] It will only be activated when the pressure at the surface (upper node) is larger than 0.0 cm
c.top = 50 #days
ponding = function(psi)
{
if (psi > 0)
{
return(-psi/c.top)
}else{
return(0)
}
}Adjusting the model
We need to replace the set.BC.fixedflux upper boundary
to an add.pointflux function and apply the ponding function
to the upper node with set.BC.fluxstate.
set.BC.fluxstate(model = UnsatStat, where = "right", func = ponding)
add.pointflux(model = UnsatStat, at = 0, value = "ppt", "precipitation")
summary(UnsatStat)one dimensional flow model
name: Unsaturated zone, 1D stationary
model has 1 spatial flux, with name
Transpiration
model has 1 point flux, with name
precipitation
left BC is of type state-flux
right BC is of type state-flux
numerical model of type FEquadratic1D
number of nodes is 222 solve.steps(UnsatStat, verboselevel = 1)[1] "stopped because of small RMSM= 0.0000000000000189448335251474"
$RMSM
[1] 0.00000000000001894483
$MAM
[1] 0.00000000000009915679plot.wtable(UnsatStat)
plot.theta(UnsatStat)
The transient model
When the dynamics are introduced into the model times series of different aspects could be analyzed. Introducing dynamics means actually that changes in time are considered. This impacts in the first place the amount of water stored in the model: \(\frac{\partial \theta}{\partial t}\).
The basic equation to be solved (Richard’s equation) is: \[ \frac{\partial \theta(\psi)}{\partial t} = -k(\psi)\left [\frac{\partial \psi}{\partial z} + 1 \right] \]
The moisture (storage) change over time (: \(\frac{\partial \theta(\psi)}{\partial t}\))
is implemented as a spatial flux of time and space Since we
are calculating \(\psi\) as the state
of the model, we are going to use the
soiltype$theta.fun(psi) function to determine the amount of
moisture content based on the current pressure head (\(\psi\))
This can be implemented numerically as:
\[
Q_{storage} = \frac{\theta \left (\psi^{t+\Delta t}\right ) - \theta
\left( \psi^t \right ) }{\Delta t}\Delta z
\] The new/current pressure head is the one which is evaluated at
the new time step: \(\psi^{t+\Delta
t}\) but since we also need to consider the previous pressure
head \(\psi^t\), we will save these
pressure just before we start with a new time step and call this
psi.old.
When more water is stored in the model, it means that less is available
for flow. This means that a minus sign should be introduced. This could
easily be accomplished by switching both pressure heads in the
aforementioned equation.
storage = function(z,state)
{
storage = (soiltype$theta.fun(psi.old(z)) - soiltype$theta.fun(state))/delta.t
}Be aware the storage function will be multiplied with
the nodal distance internally to come to the unit of L/T (cm/d).
Setting up the transient model
A copy of the stationary model is used to quickly set up the transient model
UnsatTrans = copy.model(UnsatStat)
add.spatialflux(UnsatTrans, rate = storage, name = "Storage")
do.initialize(UnsatTrans, init = drn.psi.equi)
psi.old = state.fun(UnsatTrans) #setting the strating state of the transient simulation to determine the first storage changeCreating a time loop for the calculations
A function to write the required out put in different data.frames for analysis and inspection later on.
save.data = function(model = modelname)
{
df.state = dataframe.states(model)
state.tmp <<- rbind(state.tmp,df.state$state) #adding the states (psi) to this data.frame
df.extern = dataframe.externalfluxes(model)
transp.tmp <<- rbind(transp.tmp,df.extern$`S{Transpiration}dx`) #for transpiration(roots) data
moist.tmp <<- rbind(moist.tmp,df.extern$`S{Storage}dx`)
#balance terms transp:[2,3],flow2sto:[3,2],sto2flow[3,3]
df.balance = dataframe.balance(model)
#the drainage/upward seepage at the bottom will be written in balance.dat
df.bound = dataframe.boundaries(model)
balance.tmp <<- rbind(balance.tmp,c(precip[current.time],evapot[current.time],
df.balance[2,3],df.balance[3,2],
df.balance[3,3],df.bound[1,4],df.bound[2,4]))
}Reading some meteo data from KNMI, precipitation and potential evapotranspiration In the file it reads: RH = Etmaalsom van de neerslag (in 0.1 mm) (-1 voor <0.05 mm) / Daily precipitation amount (in 0.1 mm) (-1 for <0.05 mm) EV24 = Referentiegewasverdamping (Makkink) (in 0.1 mm) / Potential evapotranspiration (Makkink) (in 0.1 mm)
meteo = read.table(file = "ppt_epot2020.txt", header = TRUE)
names(meteo) = c("date","ppt","epot")
#adjusting the unit to cm/day
precip = meteo$ppt/100
evapot = meteo$epot/100
#correcting for the -1 symbol and setting these values to 0
precip[which(precip < 0)] = 0
evapot[which(evapot < 0)] = 0Base run
begin.time = 1
current.time = begin.time
end.time = length(meteo$date)
delta.t = 1 #day, i.e. the time stepping
####containers to save intermediate results
state.tmp = c()
transp.tmp = c()
moist.tmp = c()
balance.tmp = c()
##set the initial state to equilibrium
do.initialize(UnsatTrans, init = drn.psi.equi)
start.run.time = Sys.time()
while(current.time <= end.time)
{
ppt = precip[as.integer(current.time)]
Epot = evapot[as.integer(current.time)]
psi.old = state.fun(UnsatTrans)
solve.steps(UnsatTrans,maxstatechangecrit = 1e-04,verboselevel = 1)
plot.wtable(UnsatTrans)
cat(paste("day nr.:",current.time, "\n"))
save.data(UnsatTrans)
current.time = current.time + delta.t
####plotting some stuff
# old.par = par(no.readonly=TRUE) #to save current plot settings and stuff
# layout = matrix(c(2,2),ncol = 2)#, byrow = TRUE)
# plot.roots(UnsatTrans)
# plot.theta(UnsatTrans)
# plot.wtable(UnsatTrans)
# par = old.par
#####saving the results
}[1] "stopped because of small RMSM= 0.00000322321410844525"
day nr.: 1
[1] "iteration 1 ; RMSM= 0.0182828895561381 ; MAM= 0.272386691455084"
[1] "iteration 2 ; RMSM= 0.013451485761806 ; MAM= 0.19911484630216"
[1] "iteration 3 ; RMSM= 0.00158347433531707 ; MAM= 0.0215915969539963"
[1] "iteration 4 ; RMSM= 0.0000199442197123952 ; MAM= 0.000167208936254437"
[1] "stopped because of small RMSM= 0.000000000411115759592144"
day nr.: 2
[1] "iteration 1 ; RMSM= 0.00091407662695114 ; MAM= 0.0105401176660775"
[1] "stopped because of small RMSM= 0.00000262659969892108"
day nr.: 3
[1] "iteration 1 ; RMSM= 0.00217355015216568 ; MAM= 0.0321369876938789"
[1] "iteration 2 ; RMSM= 0.0000462437775382547 ; MAM= 0.000667197626594354"
[1] "stopped because of small RMSM= 0.0000000109852168170647"
day nr.: 4
[1] "iteration 1 ; RMSM= 0.000074336801686789 ; MAM= 0.000715515884358295"
[1] "stopped because of small RMSM= 0.0000000391795042591672"
day nr.: 5
[1] "iteration 1 ; RMSM= 0.000521572505451668 ; MAM= 0.00753937130475631"
[1] "stopped because of small RMSM= 0.00000509753515243986"
day nr.: 6
[1] "iteration 1 ; RMSM= 0.0182675198053618 ; MAM= 0.272085352357907"
[1] "iteration 2 ; RMSM= 0.00943863502479566 ; MAM= 0.139567732872453"
[1] "iteration 3 ; RMSM= 0.000481180602958321 ; MAM= 0.00702202177191086"
[1] "stopped because of small RMSM= 0.000000982593167676514"
day nr.: 7
[1] "iteration 1 ; RMSM= 0.00168287839667422 ; MAM= 0.0215725177881581"
[1] "iteration 2 ; RMSM= 0.000166480153646534 ; MAM= 0.00236020628807465"
[1] "stopped because of small RMSM= 0.00000111915292686898"
day nr.: 8
[1] "iteration 1 ; RMSM= 0.00108870462609696 ; MAM= 0.0150131136787417"
[1] "stopped because of small RMSM= 0.00000460986461790549"
day nr.: 9
[1] "iteration 1 ; RMSM= 0.00449942216181567 ; MAM= 0.0666847571919904"
[1] "iteration 2 ; RMSM= 0.000221120238100151 ; MAM= 0.00319911940791098"
[1] "stopped because of small RMSM= 0.000000303411962925775"
day nr.: 10
[1] "iteration 1 ; RMSM= 0.0019307054794758 ; MAM= 0.0280741885903477"
[1] "iteration 2 ; RMSM= 0.00001903805809219 ; MAM= 0.000262070172301287"
[1] "stopped because of small RMSM= 0.00000000101207236306848"
day nr.: 11
[1] "iteration 1 ; RMSM= 0.0026292718284041 ; MAM= 0.0387582633731917"
[1] "iteration 2 ; RMSM= 0.000076186048044427 ; MAM= 0.00109094486394474"
[1] "stopped because of small RMSM= 0.0000000349454506747312"
day nr.: 12
[1] "iteration 1 ; RMSM= 0.0413298340234303 ; MAM= 0.615620101320226"
[1] "iteration 2 ; RMSM= 0.028976488110374 ; MAM= 0.431619902294365"
[1] "iteration 3 ; RMSM= 0.0150571350515181 ; MAM= 0.209283098174662"
[1] "iteration 4 ; RMSM= 0.00166741297120585 ; MAM= 0.0185360779771022"
[1] "iteration 5 ; RMSM= 0.0000529747269548393 ; MAM= 0.000683683932622969"
[1] "stopped because of small RMSM= 0.000000000451778817907105"
day nr.: 13
[1] "iteration 1 ; RMSM= 0.00868858865772752 ; MAM= 0.126533150099178"
[1] "iteration 2 ; RMSM= 0.00025291005473047 ; MAM= 0.0036214851469894"
[1] "stopped because of small RMSM= 0.000000135041670575866"
day nr.: 14
[1] "iteration 1 ; RMSM= 0.00613387747543616 ; MAM= 0.0909622207323768"
[1] "iteration 2 ; RMSM= 0.000478782594115609 ; MAM= 0.00703789815417716"
[1] "stopped because of small RMSM= 0.00000154698149809441"
day nr.: 15
[1] "iteration 1 ; RMSM= 0.000404144037735575 ; MAM= 0.00249193784200796"
[1] "stopped because of small RMSM= 0.00000642512358115432"
day nr.: 16
[1] "iteration 1 ; RMSM= 0.00313451525678383 ; MAM= 0.0460846405550748"
[1] "iteration 2 ; RMSM= 0.0000491709831992355 ; MAM= 0.000684025633769042"
[1] "stopped because of small RMSM= 0.00000000660183383423733"
day nr.: 17
[1] "iteration 1 ; RMSM= 0.00241011615234775 ; MAM= 0.0351519283550822"
[1] "iteration 2 ; RMSM= 0.0000588642801156408 ; MAM= 0.000842360066469334"
[1] "stopped because of small RMSM= 0.0000000179914082905065"
day nr.: 18
[1] "iteration 1 ; RMSM= 0.000281798916568899 ; MAM= 0.0032076643055098"
[1] "stopped because of small RMSM= 0.00000255283350143539"
day nr.: 19
[1] "iteration 1 ; RMSM= 0.0000968910434328737 ; MAM= 0.000467561305880791"
[1] "stopped because of small RMSM= 0.000000255674019476971"
day nr.: 20
[1] "iteration 1 ; RMSM= 0.00388263770749124 ; MAM= 0.0572980719534473"
[1] "iteration 2 ; RMSM= 0.000106296330472565 ; MAM= 0.00152780945867878"
[1] "stopped because of small RMSM= 0.0000000464277090251021"
day nr.: 21
[1] "iteration 1 ; RMSM= 0.00206334209197211 ; MAM= 0.0304448894127515"
[1] "iteration 2 ; RMSM= 0.0000489590982738735 ; MAM= 0.000700342305397905"
[1] "stopped because of small RMSM= 0.0000000129863108439463"
day nr.: 22
[1] "iteration 1 ; RMSM= 0.0000698731204624518 ; MAM= 0.000449616544120811"
[1] "stopped because of small RMSM= 0.000000365248076885782"
day nr.: 23
[1] "iteration 1 ; RMSM= 0.0000571370238377464 ; MAM= 0.000364389349088688"
[1] "stopped because of small RMSM= 0.0000000363299247346774"
day nr.: 24
[1] "iteration 1 ; RMSM= 0.0000516105963806834 ; MAM= 0.000465390650994726"
[1] "stopped because of small RMSM= 0.00000000977497127484989"
day nr.: 25
[1] "iteration 1 ; RMSM= 0.0405225412622045 ; MAM= 0.6032985012632"
[1] "iteration 2 ; RMSM= 0.0288964197246736 ; MAM= 0.429825069496442"
[1] "iteration 3 ; RMSM= 0.0198287708452361 ; MAM= 0.283146149041416"
[1] "iteration 4 ; RMSM= 0.002518306893147 ; MAM= 0.0274088942130284"
[1] "iteration 5 ; RMSM= 0.0000726249059974557 ; MAM= 0.000745579605956337"
[1] "stopped because of small RMSM= 0.00000000430149223930456"
day nr.: 26
[1] "iteration 1 ; RMSM= 0.00970031853707168 ; MAM= 0.139538338673209"
[1] "iteration 2 ; RMSM= 0.000339267788667207 ; MAM= 0.00472226505040557"
[1] "stopped because of small RMSM= 0.000000262405371654145"
day nr.: 27
[1] "iteration 1 ; RMSM= 0.000349045266863981 ; MAM= 0.00275470566651842"
[1] "stopped because of small RMSM= 0.00000284023937862365"
day nr.: 28
[1] "iteration 1 ; RMSM= 0.00111744520687165 ; MAM= 0.0151846753949046"
[1] "iteration 2 ; RMSM= 0.0000292087720668716 ; MAM= 0.000277365652223333"
[1] "stopped because of small RMSM= 0.00000000326263220101139"
day nr.: 29
[1] "iteration 1 ; RMSM= 0.00250085993601727 ; MAM= 0.0352812761662704"
[1] "iteration 2 ; RMSM= 0.000089142177300399 ; MAM= 0.00127743431794027"
[1] "stopped because of small RMSM= 0.0000000513336103391081"
day nr.: 30
[1] "iteration 1 ; RMSM= 0.00100052034886044 ; MAM= 0.0143739596392925"
[1] "stopped because of small RMSM= 0.00000357362234854109"
day nr.: 31
[1] "iteration 1 ; RMSM= 0.028457290761812 ; MAM= 0.421078693543642"
[1] "iteration 2 ; RMSM= 0.00380092920918391 ; MAM= 0.0534227001577237"
[1] "iteration 3 ; RMSM= 0.000122356322177937 ; MAM= 0.00143057910555378"
[1] "stopped because of small RMSM= 0.0000000156522014096964"
day nr.: 32
[1] "iteration 1 ; RMSM= 0.00409216452682703 ; MAM= 0.0604232833222159"
[1] "iteration 2 ; RMSM= 0.0000877181797132916 ; MAM= 0.00126539563394262"
[1] "stopped because of small RMSM= 0.000000022921149114888"
day nr.: 33
[1] "iteration 1 ; RMSM= 0.00763415874749484 ; MAM= 0.112233272534472"
[1] "iteration 2 ; RMSM= 0.00024440062189589 ; MAM= 0.003558961776293"
[1] "stopped because of small RMSM= 0.000000174041641628738"
day nr.: 34
[1] "iteration 1 ; RMSM= 0.0105718377795943 ; MAM= 0.156848364828656"
[1] "iteration 2 ; RMSM= 0.000958632797662269 ; MAM= 0.0139590145652104"
[1] "stopped because of small RMSM= 0.00000477293463834725"
day nr.: 35
[1] "iteration 1 ; RMSM= 0.000472152810730493 ; MAM= 0.00265143649898664"
[1] "iteration 2 ; RMSM= 0.0000114338315778698 ; MAM= 0.000135071581581658"
[1] "stopped because of small RMSM= 0.000000000605519154608615"
day nr.: 36
[1] "iteration 1 ; RMSM= 0.000320111364483511 ; MAM= 0.0022959798171323"
[1] "stopped because of small RMSM= 0.00000338471356862579"
day nr.: 37
[1] "iteration 1 ; RMSM= 0.000215276318191469 ; MAM= 0.00179794802863026"
[1] "stopped because of small RMSM= 0.000000236821639370442"
day nr.: 38
[1] "iteration 1 ; RMSM= 0.0299448569931019 ; MAM= 0.445989067636822"
[1] "iteration 2 ; RMSM= 0.0251808142862868 ; MAM= 0.37331601101589"
[1] "iteration 3 ; RMSM= 0.00317967261249034 ; MAM= 0.0434440198612641"
[1] "iteration 4 ; RMSM= 0.0000304966850174412 ; MAM= 0.000404508864041531"
[1] "stopped because of small RMSM= 0.00000000174496257250462"
day nr.: 39
[1] "iteration 1 ; RMSM= 0.00185806959161471 ; MAM= 0.0253167719417486"
[1] "stopped because of small RMSM= 0.00000646654941555628"
day nr.: 40
[1] "iteration 1 ; RMSM= 0.00491626026961275 ; MAM= 0.0728339587267981"
[1] "iteration 2 ; RMSM= 0.000220614225300766 ; MAM= 0.00313781735352968"
[1] "stopped because of small RMSM= 0.000000263619575947359"
day nr.: 41
[1] "iteration 1 ; RMSM= 0.00108524076981546 ; MAM= 0.0149572427854619"
[1] "iteration 2 ; RMSM= 0.0000251981679494888 ; MAM= 0.000305671567223194"
[1] "stopped because of small RMSM= 0.00000000532702022473557"
day nr.: 42
[1] "iteration 1 ; RMSM= 0.0581494345440423 ; MAM= 0.865294910117503"
[1] "iteration 2 ; RMSM= 0.0402038585211219 ; MAM= 0.597790777514407"
[1] "iteration 3 ; RMSM= 0.0204877132925039 ; MAM= 0.265398632463941"
[1] "iteration 4 ; RMSM= 0.00350202808909465 ; MAM= 0.0328715659713847"
[1] "iteration 5 ; RMSM= 0.000197030097241352 ; MAM= 0.00232386984304189"
[1] "stopped because of small RMSM= 0.0000000058930425212094"
day nr.: 43
[1] "iteration 1 ; RMSM= 0.0397640550493341 ; MAM= 0.578513332545982"
[1] "iteration 2 ; RMSM= 0.0076080779987671 ; MAM= 0.109541763057415"
[1] "iteration 3 ; RMSM= 0.000215927608465635 ; MAM= 0.00290491357962552"
[1] "stopped because of small RMSM= 0.0000000961111245932276"
day nr.: 44
[1] "iteration 1 ; RMSM= 0.00145558930658543 ; MAM= 0.0122389571107533"
[1] "iteration 2 ; RMSM= 0.0000586110762627069 ; MAM= 0.000373738010686136"
[1] "stopped because of small RMSM= 0.0000000239664512119916"
day nr.: 45
[1] "iteration 1 ; RMSM= 0.0421769573851926 ; MAM= 0.628004698066925"
[1] "iteration 2 ; RMSM= 0.0328733280193065 ; MAM= 0.486625944828488"
[1] "iteration 3 ; RMSM= 0.00359237928532264 ; MAM= 0.0512758377379411"
[1] "iteration 4 ; RMSM= 0.0000289453731715396 ; MAM= 0.000413025648488996"
[1] "stopped because of small RMSM= 0.00000000125233714313075"
day nr.: 46
[1] "iteration 1 ; RMSM= 0.010779776799406 ; MAM= 0.158643981677956"
[1] "iteration 2 ; RMSM= 0.000626869475820453 ; MAM= 0.00909605449757867"
[1] "stopped because of small RMSM= 0.0000023518897769476"
day nr.: 47
[1] "iteration 1 ; RMSM= 0.000455444625377607 ; MAM= 0.00431845039290085"
[1] "stopped because of small RMSM= 0.00000423402724071677"
day nr.: 48
[1] "iteration 1 ; RMSM= 0.000972694217735456 ; MAM= 0.0138417884579942"
[1] "stopped because of small RMSM= 0.00000450590603054167"
day nr.: 49
[1] "iteration 1 ; RMSM= 0.00112636984240674 ; MAM= 0.0117891667718262"
[1] "stopped because of small RMSM= 0.00000331743339837044"
day nr.: 50
[1] "iteration 1 ; RMSM= 0.00811022350528331 ; MAM= 0.119824186739318"
[1] "iteration 2 ; RMSM= 0.000623296861807146 ; MAM= 0.00908266782978785"
[1] "stopped because of small RMSM= 0.00000184463707520503"
day nr.: 51
[1] "iteration 1 ; RMSM= 0.0310783566504819 ; MAM= 0.462700892392869"
[1] "iteration 2 ; RMSM= 0.0170360227752911 ; MAM= 0.251445031043395"
[1] "iteration 3 ; RMSM= 0.000913256532237714 ; MAM= 0.0133265395129141"
[1] "stopped because of small RMSM= 0.00000207428575181735"
day nr.: 52
[1] "iteration 1 ; RMSM= 0.0960828476681933 ; MAM= 1.42928818279238"
[1] "iteration 2 ; RMSM= 0.0683823207875672 ; MAM= 1.01616336718645"
[1] "iteration 3 ; RMSM= 0.0460974445027747 ; MAM= 0.55904150486702"
[1] "iteration 4 ; RMSM= 0.00709214558572257 ; MAM= 0.0610982014524595"
[1] "iteration 5 ; RMSM= 0.000131701761902927 ; MAM= 0.000984091126056352"
[1] "stopped because of small RMSM= 0.00000000226057144289622"
day nr.: 53
[1] "iteration 1 ; RMSM= 0.0373268801597803 ; MAM= 0.536704799036815"
[1] "iteration 2 ; RMSM= 0.00380805535684331 ; MAM= 0.0505670261904031"
[1] "iteration 3 ; RMSM= 0.0000206554976824361 ; MAM= 0.000287991552654465"
[1] "stopped because of small RMSM= 0.000000000415312732901519"
day nr.: 54
[1] "iteration 1 ; RMSM= 0.0144582554582232 ; MAM= 0.205485453044034"
[1] "iteration 2 ; RMSM= 0.00164455316142716 ; MAM= 0.0233424967505742"
[1] "iteration 3 ; RMSM= 0.0000181622515416537 ; MAM= 0.000167891274194976"
[1] "stopped because of small RMSM= 0.000000000320519415305602"
day nr.: 55
[1] "iteration 1 ; RMSM= 0.0104817642402698 ; MAM= 0.148417898389586"
[1] "iteration 2 ; RMSM= 0.000357103689382116 ; MAM= 0.00514257662676731"
[1] "stopped because of small RMSM= 0.000000166954103997736"
day nr.: 56
[1] "iteration 1 ; RMSM= 0.00116090847413646 ; MAM= 0.0136780608851809"
[1] "stopped because of small RMSM= 0.00000505295802382615"
day nr.: 57
[1] "iteration 1 ; RMSM= 0.00466765845853577 ; MAM= 0.0664833426466207"
[1] "iteration 2 ; RMSM= 0.0000684722155795635 ; MAM= 0.000938334995107648"
[1] "stopped because of small RMSM= 0.0000000087098058889588"
day nr.: 58
[1] "iteration 1 ; RMSM= 0.0114606121358243 ; MAM= 0.168884044725215"
[1] "iteration 2 ; RMSM= 0.000801226526917987 ; MAM= 0.0117622827607976"
[1] "stopped because of small RMSM= 0.00000286297904306784"
day nr.: 59
[1] "iteration 1 ; RMSM= 0.00387523315031232 ; MAM= 0.0413424550313343"
[1] "iteration 2 ; RMSM= 0.000137967521663226 ; MAM= 0.00171892119995541"
[1] "stopped because of small RMSM= 0.000000115870269767348"
day nr.: 60
[1] "iteration 1 ; RMSM= 0.0359517204617781 ; MAM= 0.534722321434795"
[1] "iteration 2 ; RMSM= 0.0191109749195852 ; MAM= 0.282177711401045"
[1] "iteration 3 ; RMSM= 0.00118017797600924 ; MAM= 0.0173856975853356"
[1] "stopped because of small RMSM= 0.00000339476403421663"
day nr.: 61
[1] "iteration 1 ; RMSM= 0.0113022776939306 ; MAM= 0.16602424607089"
[1] "iteration 2 ; RMSM= 0.0008486381379158 ; MAM= 0.012282450628287"
[1] "stopped because of small RMSM= 0.00000297778302163283"
day nr.: 62
[1] "iteration 1 ; RMSM= 0.00246728810317353 ; MAM= 0.0245227883784592"
[1] "iteration 2 ; RMSM= 0.0000250669701099413 ; MAM= 0.000237874040802395"
[1] "stopped because of small RMSM= 0.00000000414330351698936"
day nr.: 63
[1] "iteration 1 ; RMSM= 0.0489287649466724 ; MAM= 0.728414202604777"
[1] "iteration 2 ; RMSM= 0.0438073092403503 ; MAM= 0.649272620889151"
[1] "iteration 3 ; RMSM= 0.00596421115915699 ; MAM= 0.0850097829215943"
[1] "iteration 4 ; RMSM= 0.000118320115061125 ; MAM= 0.00107866115501509"
[1] "stopped because of small RMSM= 0.0000000107846735435764"
day nr.: 64
[1] "iteration 1 ; RMSM= 0.0119416363373792 ; MAM= 0.174911409123773"
[1] "iteration 2 ; RMSM= 0.000599896140542837 ; MAM= 0.00861423512375781"
[1] "stopped because of small RMSM= 0.000000962534775303353"
day nr.: 65
[1] "iteration 1 ; RMSM= 0.00647084660882413 ; MAM= 0.0895068477713284"
[1] "iteration 2 ; RMSM= 0.000444242771547945 ; MAM= 0.0063250555455426"
[1] "stopped because of small RMSM= 0.00000132235662629663"
day nr.: 66
[1] "iteration 1 ; RMSM= 0.00195012173834752 ; MAM= 0.0222860388063776"
[1] "stopped because of small RMSM= 0.00000622449973671488"
day nr.: 67
[1] "iteration 1 ; RMSM= 0.000939548828272307 ; MAM= 0.0130416744967529"
[1] "stopped because of small RMSM= 0.00000430784140276577"
day nr.: 68
[1] "iteration 1 ; RMSM= 0.0613502233989791 ; MAM= 0.913713885126754"
[1] "iteration 2 ; RMSM= 0.0536709002509876 ; MAM= 0.795967283638241"
[1] "iteration 3 ; RMSM= 0.00714010448344881 ; MAM= 0.0972444989760173"
[1] "iteration 4 ; RMSM= 0.000329552911176694 ; MAM= 0.00396567561087457"
[1] "stopped because of small RMSM= 0.0000000693170722151981"
day nr.: 69
[1] "iteration 1 ; RMSM= 0.00391146029350082 ; MAM= 0.0319749945347325"
[1] "iteration 2 ; RMSM= 0.000123151000934344 ; MAM= 0.00137213486183096"
[1] "stopped because of small RMSM= 0.0000000080995436515593"
day nr.: 70
[1] "iteration 1 ; RMSM= 0.00547754974930795 ; MAM= 0.0512274112113084"
[1] "iteration 2 ; RMSM= 0.0000527944769257896 ; MAM= 0.000645756119928009"
[1] "stopped because of small RMSM= 0.00000000471391347539727"
day nr.: 71
[1] "iteration 1 ; RMSM= 0.00571382211943223 ; MAM= 0.058692299002062"
[1] "iteration 2 ; RMSM= 0.0000666165515023639 ; MAM= 0.000923830532231884"
[1] "stopped because of small RMSM= 0.000000165732956126755"
day nr.: 72
[1] "iteration 1 ; RMSM= 0.010097456152152 ; MAM= 0.143288407002839"
[1] "iteration 2 ; RMSM= 0.000892439011867346 ; MAM= 0.0127004039196407"
[1] "stopped because of small RMSM= 0.00000496223282126756"
day nr.: 73
[1] "iteration 1 ; RMSM= 0.00232282369252395 ; MAM= 0.021893960558014"
[1] "iteration 2 ; RMSM= 0.000243581028372565 ; MAM= 0.00262030967909299"
[1] "stopped because of small RMSM= 0.00000000295920644838636"
day nr.: 74
[1] "iteration 1 ; RMSM= 0.00289079475934319 ; MAM= 0.0122263764960786"
[1] "iteration 2 ; RMSM= 0.00125283187290468 ; MAM= 0.0124770205487209"
[1] "stopped because of small RMSM= 0.00000676253009548834"
day nr.: 75
[1] "iteration 1 ; RMSM= 0.00287294754221042 ; MAM= 0.00711157304518995"
[1] "iteration 2 ; RMSM= 0.000487264908535496 ; MAM= 0.00294317469394323"
[1] "stopped because of small RMSM= 0.00000250834925776388"
day nr.: 76
[1] "iteration 1 ; RMSM= 0.00349723678325041 ; MAM= 0.0413345209512455"
[1] "iteration 2 ; RMSM= 0.0000257300571093981 ; MAM= 0.000286678195626111"
[1] "stopped because of small RMSM= 0.000000000044375996959401"
day nr.: 77
[1] "iteration 1 ; RMSM= 0.00291162429834258 ; MAM= 0.0364839173713389"
[1] "stopped because of small RMSM= 0.00000768984248085162"
day nr.: 78
[1] "iteration 1 ; RMSM= 0.00175367716873295 ; MAM= 0.0221794751235877"
[1] "stopped because of small RMSM= 0.0000044473777770293"
day nr.: 79
[1] "iteration 1 ; RMSM= 0.000768943753952819 ; MAM= 0.00785395933315713"
[1] "stopped because of small RMSM= 0.00000420662368629997"
day nr.: 80
[1] "iteration 1 ; RMSM= 0.00169285980108814 ; MAM= 0.0199648403110811"
[1] "stopped because of small RMSM= 0.00000292098188261012"
day nr.: 81
[1] "iteration 1 ; RMSM= 0.00125114133751795 ; MAM= 0.0154772710466494"
[1] "stopped because of small RMSM= 0.00000334695108999061"
day nr.: 82
[1] "iteration 1 ; RMSM= 0.0004229702619332 ; MAM= 0.00537037555664788"
[1] "stopped because of small RMSM= 0.0000016236987280287"
day nr.: 83
[1] "iteration 1 ; RMSM= 0.00127711412439499 ; MAM= 0.0160141945053487"
[1] "stopped because of small RMSM= 0.00000278779654530596"
day nr.: 84
[1] "iteration 1 ; RMSM= 0.0002199948952121 ; MAM= 0.00127439512290908"
[1] "stopped because of small RMSM= 0.00000101200780431157"
day nr.: 85
[1] "iteration 1 ; RMSM= 0.000673902114443306 ; MAM= 0.00844290854630414"
[1] "stopped because of small RMSM= 0.00000226802332759756"
day nr.: 86
[1] "iteration 1 ; RMSM= 0.000273012969801386 ; MAM= 0.0025710892731908"
[1] "stopped because of small RMSM= 0.000000926468729860236"
day nr.: 87
[1] "iteration 1 ; RMSM= 0.000550509656150976 ; MAM= 0.00641993957274023"
[1] "stopped because of small RMSM= 0.000000555691349714486"
day nr.: 88
[1] "iteration 1 ; RMSM= 0.000455341285960915 ; MAM= 0.00539957717055736"
[1] "stopped because of small RMSM= 0.000000258234573774332"
day nr.: 89
[1] "iteration 1 ; RMSM= 0.000167560353317906 ; MAM= 0.000859601950669725"
[1] "stopped because of small RMSM= 0.00000160098950714768"
day nr.: 90
[1] "iteration 1 ; RMSM= 0.000340841553776057 ; MAM= 0.00383564756467017"
[1] "stopped because of small RMSM= 0.000000209790445046338"
day nr.: 91
[1] "iteration 1 ; RMSM= 0.000256539471109255 ; MAM= 0.00251866226026725"
[1] "stopped because of small RMSM= 0.000000408989345056322"
day nr.: 92
[1] "iteration 1 ; RMSM= 0.00211770405108111 ; MAM= 0.0313355989057188"
[1] "iteration 2 ; RMSM= 0.0000642049342248419 ; MAM= 0.000934197044198018"
[1] "stopped because of small RMSM= 0.0000000350775429381198"
day nr.: 93
[1] "iteration 1 ; RMSM= 0.00179524561689663 ; MAM= 0.0265685618360582"
[1] "iteration 2 ; RMSM= 0.0000654289187116251 ; MAM= 0.000943526073825353"
[1] "stopped because of small RMSM= 0.0000000451718548539042"
day nr.: 94
[1] "iteration 1 ; RMSM= 0.000261499401562116 ; MAM= 0.00164128069795623"
[1] "stopped because of small RMSM= 0.00000648205110840125"
day nr.: 95
[1] "iteration 1 ; RMSM= 0.000312788314932477 ; MAM= 0.00158935789903743"
[1] "stopped because of small RMSM= 0.00000622554434526112"
day nr.: 96
[1] "iteration 1 ; RMSM= 0.000317196380963883 ; MAM= 0.00126669153786411"
[1] "stopped because of small RMSM= 0.00000503774932267882"
day nr.: 97
[1] "iteration 1 ; RMSM= 0.000366948287953204 ; MAM= 0.00204262392465182"
[1] "stopped because of small RMSM= 0.00000521869747348959"
day nr.: 98
[1] "iteration 1 ; RMSM= 0.00039525840456695 ; MAM= 0.00217476389820934"
[1] "stopped because of small RMSM= 0.00000277853252331977"
day nr.: 99
[1] "iteration 1 ; RMSM= 0.000399724752722805 ; MAM= 0.00238584009861242"
[1] "stopped because of small RMSM= 0.00000542805019476261"
day nr.: 100
[1] "iteration 1 ; RMSM= 0.000449489036044472 ; MAM= 0.00240011900304821"
[1] "stopped because of small RMSM= 0.00000544287004293538"
day nr.: 101
[1] "iteration 1 ; RMSM= 0.000355854099523997 ; MAM= 0.00178005973693092"
[1] "stopped because of small RMSM= 0.00000324902814132643"
day nr.: 102
[1] "iteration 1 ; RMSM= 0.000285937307176073 ; MAM= 0.00125936878406011"
[1] "stopped because of small RMSM= 0.000000212225182948148"
day nr.: 103
[1] "iteration 1 ; RMSM= 0.000190034615438508 ; MAM= 0.00140109708202006"
[1] "stopped because of small RMSM= 0.000000873989551233943"
day nr.: 104
[1] "iteration 1 ; RMSM= 0.000378609098621325 ; MAM= 0.00210340801172448"
[1] "stopped because of small RMSM= 0.0000038760488256403"
day nr.: 105
[1] "iteration 1 ; RMSM= 0.000515882226061386 ; MAM= 0.00279174945297456"
[1] "stopped because of small RMSM= 0.00000978931389538086"
day nr.: 106
[1] "iteration 1 ; RMSM= 0.000472540918559116 ; MAM= 0.00252453406573001"
[1] "stopped because of small RMSM= 0.00000675277630034524"
day nr.: 107
[1] "iteration 1 ; RMSM= 0.0067850127145505 ; MAM= 0.100383392178941"
[1] "iteration 2 ; RMSM= 0.00613181478341326 ; MAM= 0.0900622212522981"
[1] "iteration 3 ; RMSM= 0.000926388571023517 ; MAM= 0.0123614617569509"
[1] "iteration 4 ; RMSM= 0.000013692312910836 ; MAM= 0.000135751115042715"
[1] "stopped because of small RMSM= 0.000000000978696500164978"
day nr.: 108
[1] "iteration 1 ; RMSM= 0.00294085659499999 ; MAM= 0.0424973386783033"
[1] "iteration 2 ; RMSM= 0.000321272621256435 ; MAM= 0.00444334966963866"
[1] "stopped because of small RMSM= 0.00000232647514941181"
day nr.: 109
[1] "iteration 1 ; RMSM= 0.000484367246883848 ; MAM= 0.00281952074278857"
[1] "iteration 2 ; RMSM= 0.0000113144630734527 ; MAM= 0.0000770535154266533"
[1] "stopped because of small RMSM= 0.00000000310561369921125"
day nr.: 110
[1] "iteration 1 ; RMSM= 0.000602200031195184 ; MAM= 0.00341755710191724"
[1] "iteration 2 ; RMSM= 0.0000193517622021593 ; MAM= 0.000132368287094183"
[1] "stopped because of small RMSM= 0.00000000360084549781604"
day nr.: 111
[1] "iteration 1 ; RMSM= 0.000701869979475278 ; MAM= 0.00356334986019732"
[1] "iteration 2 ; RMSM= 0.0000278308434695471 ; MAM= 0.000136519791970602"
[1] "stopped because of small RMSM= 0.000000208842107804125"
day nr.: 112
[1] "iteration 1 ; RMSM= 0.000698752615057275 ; MAM= 0.00305162929582906"
[1] "iteration 2 ; RMSM= 0.00003614092939775 ; MAM= 0.000261637957433829"
[1] "stopped because of small RMSM= 0.000000103803992454556"
day nr.: 113
[1] "iteration 1 ; RMSM= 0.00063486224875104 ; MAM= 0.00294812160245746"
[1] "iteration 2 ; RMSM= 0.0000282002989874091 ; MAM= 0.000257355307443884"
[1] "stopped because of small RMSM= 0.0000000116099926826205"
day nr.: 114
[1] "iteration 1 ; RMSM= 0.000307354850711933 ; MAM= 0.00159138547560161"
[1] "stopped because of small RMSM= 0.0000091927996906428"
day nr.: 115
[1] "iteration 1 ; RMSM= 0.000687264025846672 ; MAM= 0.00269951800382799"
[1] "iteration 2 ; RMSM= 0.0000367143303187267 ; MAM= 0.000245202595986592"
[1] "stopped because of small RMSM= 0.0000000323588742864625"
day nr.: 116
[1] "iteration 1 ; RMSM= 0.000719613657282257 ; MAM= 0.00270858955262827"
[1] "iteration 2 ; RMSM= 0.0000545240545441972 ; MAM= 0.000326109552291898"
[1] "stopped because of small RMSM= 0.000000520172134599785"
day nr.: 117
[1] "iteration 1 ; RMSM= 0.0126411412492505 ; MAM= 0.18783136066481"
[1] "iteration 2 ; RMSM= 0.0038740302968768 ; MAM= 0.0543622244863703"
[1] "iteration 3 ; RMSM= 0.00256297763189835 ; MAM= 0.0318181239413954"
[1] "iteration 4 ; RMSM= 0.00212625465441936 ; MAM= 0.0264749958447588"
[1] "iteration 5 ; RMSM= 0.00180748442771951 ; MAM= 0.0240122023982537"
[1] "iteration 6 ; RMSM= 0.00142758195302233 ; MAM= 0.0188724165441217"
[1] "iteration 7 ; RMSM= 0.0013009888067183 ; MAM= 0.0128793144157704"
[1] "iteration 8 ; RMSM= 0.000100790294752895 ; MAM= 0.000882697177382957"
[1] "stopped because of small RMSM= 0.0000000432028790430684"
day nr.: 118
[1] "iteration 1 ; RMSM= 0.0125415952059239 ; MAM= 0.170736502017439"
[1] "iteration 2 ; RMSM= 0.00917766929430939 ; MAM= 0.125101005689377"
[1] "iteration 3 ; RMSM= 0.00710370726862193 ; MAM= 0.0965221632447585"
[1] "iteration 4 ; RMSM= 0.00550996465034193 ; MAM= 0.0749279460874589"
[1] "iteration 5 ; RMSM= 0.00473330602778833 ; MAM= 0.0355864686936545"
[1] "iteration 6 ; RMSM= 0.000545879040998098 ; MAM= 0.00383508533997912"
[1] "stopped because of small RMSM= 0.000000626049665390807"
day nr.: 119
[1] "iteration 1 ; RMSM= 0.00552331210904174 ; MAM= 0.0783071083054791"
[1] "iteration 2 ; RMSM= 0.000663679064708799 ; MAM= 0.00880282789896623"
[1] "stopped because of small RMSM= 0.00000590042348698853"
day nr.: 120
[1] "iteration 1 ; RMSM= 0.00287314485988983 ; MAM= 0.0402421419475423"
[1] "iteration 2 ; RMSM= 0.000275394418688961 ; MAM= 0.00330966538408176"
[1] "stopped because of small RMSM= 0.00000270471602852146"
day nr.: 121
[1] "iteration 1 ; RMSM= 0.0238828692666991 ; MAM= 0.351833068838063"
[1] "iteration 2 ; RMSM= 0.0174883903256304 ; MAM= 0.256203303146643"
[1] "iteration 3 ; RMSM= 0.0128767054885526 ; MAM= 0.186496880543439"
[1] "iteration 4 ; RMSM= 0.0095867860489772 ; MAM= 0.0944037837499594"
[1] "iteration 5 ; RMSM= 0.000989643392455504 ; MAM= 0.00642214433002131"
[1] "iteration 6 ; RMSM= 0.0000155866349766706 ; MAM= 0.000134443684729812"
[1] "stopped because of small RMSM= 0.0000000000502389411000808"
day nr.: 122
[1] "iteration 1 ; RMSM= 0.0169345697636649 ; MAM= 0.245009774365003"
[1] "iteration 2 ; RMSM= 0.00317745966047317 ; MAM= 0.0456114222095279"
[1] "iteration 3 ; RMSM= 0.0000872038828392025 ; MAM= 0.00115476386132891"
[1] "stopped because of small RMSM= 0.0000000376047007785093"
day nr.: 123
[1] "iteration 1 ; RMSM= 0.000542716893101139 ; MAM= 0.00474541013244283"
[1] "iteration 2 ; RMSM= 0.0000191180520154983 ; MAM= 0.000246665480416464"
[1] "stopped because of small RMSM= 0.00000000437935036196015"
day nr.: 124
[1] "iteration 1 ; RMSM= 0.00115172815513572 ; MAM= 0.017067271725746"
[1] "iteration 2 ; RMSM= 0.0000763673189720043 ; MAM= 0.00111980356479899"
[1] "stopped because of small RMSM= 0.000000147336609456161"
day nr.: 125
[1] "iteration 1 ; RMSM= 0.000329130460337613 ; MAM= 0.00256233223849119"
[1] "iteration 2 ; RMSM= 0.0000142162983242168 ; MAM= 0.0000991681252059604"
[1] "stopped because of small RMSM= 0.0000000972192521001707"
day nr.: 126
[1] "iteration 1 ; RMSM= 0.00040970188364354 ; MAM= 0.00316929646315478"
[1] "iteration 2 ; RMSM= 0.000022986557890359 ; MAM= 0.000163487035961828"
[1] "stopped because of small RMSM= 0.000000123346520924638"
day nr.: 127
[1] "iteration 1 ; RMSM= 0.00049252340549499 ; MAM= 0.00304624835500844"
[1] "iteration 2 ; RMSM= 0.0000260337374531884 ; MAM= 0.000204837634138803"
[1] "stopped because of small RMSM= 0.00000000465341238744351"
day nr.: 128
[1] "iteration 1 ; RMSM= 0.000623718581397345 ; MAM= 0.00345526266604615"
[1] "iteration 2 ; RMSM= 0.0000567453188004339 ; MAM= 0.000372771726925659"
[1] "stopped because of small RMSM= 0.0000001661498114598"
day nr.: 129
[1] "iteration 1 ; RMSM= 0.000160648421744323 ; MAM= 0.00191202987017883"
[1] "stopped because of small RMSM= 0.000000629373905535258"
day nr.: 130
[1] "iteration 1 ; RMSM= 0.000383385402172673 ; MAM= 0.00186010054038181"
[1] "stopped because of small RMSM= 0.0000024847919644016"
day nr.: 131
[1] "iteration 1 ; RMSM= 0.000202100986600875 ; MAM= 0.00290440505645366"
[1] "stopped because of small RMSM= 0.00000316477210841294"
day nr.: 132
[1] "iteration 1 ; RMSM= 0.0000977142507105659 ; MAM= 0.000584923747534341"
[1] "stopped because of small RMSM= 0.000000929781025902471"
day nr.: 133
[1] "iteration 1 ; RMSM= 0.000353220983319026 ; MAM= 0.00168112538536302"
[1] "stopped because of small RMSM= 0.00000549180898558123"
day nr.: 134
[1] "iteration 1 ; RMSM= 0.000243768453252029 ; MAM= 0.0013084171100272"
[1] "stopped because of small RMSM= 0.00000242905195201315"
day nr.: 135
[1] "iteration 1 ; RMSM= 0.000229484635680608 ; MAM= 0.00122496156418322"
[1] "stopped because of small RMSM= 0.00000568653841812584"
day nr.: 136
[1] "iteration 1 ; RMSM= 0.000381286254504528 ; MAM= 0.00178266927231892"
[1] "iteration 2 ; RMSM= 0.0000126908413796292 ; MAM= 0.0000854544347360393"
[1] "stopped because of small RMSM= 0.00000000335341634260187"
day nr.: 137
[1] "iteration 1 ; RMSM= 0.000860207135613677 ; MAM= 0.00344826707372725"
[1] "iteration 2 ; RMSM= 0.000121415315171535 ; MAM= 0.000694591776421449"
[1] "stopped because of small RMSM= 0.00000169691557419127"
day nr.: 138
[1] "iteration 1 ; RMSM= 0.000627949941961352 ; MAM= 0.00242942717928358"
[1] "iteration 2 ; RMSM= 0.0000916283484939895 ; MAM= 0.000422317770813828"
[1] "stopped because of small RMSM= 0.00000221529152087858"
day nr.: 139
[1] "iteration 1 ; RMSM= 0.000431068304982015 ; MAM= 0.00169919672415716"
[1] "iteration 2 ; RMSM= 0.0000563830372178635 ; MAM= 0.000283970798215662"
[1] "stopped because of small RMSM= 0.00000175247588954713"
day nr.: 140
[1] "iteration 1 ; RMSM= 0.000948817892857944 ; MAM= 0.00329079280920944"
[1] "iteration 2 ; RMSM= 0.000284108431312747 ; MAM= 0.00190501364256815"
[1] "iteration 3 ; RMSM= 0.0000675870139042552 ; MAM= 0.000670171187132337"
[1] "stopped because of small RMSM= 0.00000535893018087409"
day nr.: 141
[1] "iteration 1 ; RMSM= 0.000419914559277641 ; MAM= 0.00272388379534009"
[1] "iteration 2 ; RMSM= 0.0000584292799984979 ; MAM= 0.000555827988383644"
[1] "stopped because of small RMSM= 0.000000577230698492013"
day nr.: 142
[1] "iteration 1 ; RMSM= 0.00116364826424905 ; MAM= 0.00986304457743537"
[1] "iteration 2 ; RMSM= 0.000464113442950074 ; MAM= 0.00587340384580623"
[1] "iteration 3 ; RMSM= 0.000141924694713485 ; MAM= 0.00165963654502599"
[1] "stopped because of small RMSM= 0.000000474813116952884"
day nr.: 143
[1] "iteration 1 ; RMSM= 0.0334679125618319 ; MAM= 0.498507418792786"
[1] "iteration 2 ; RMSM= 0.0328811841740668 ; MAM= 0.48976626978423"
[1] "iteration 3 ; RMSM= 0.0319264182959581 ; MAM= 0.475537228171687"
[1] "iteration 4 ; RMSM= 0.0221114819620225 ; MAM= 0.212787743069558"
[1] "iteration 5 ; RMSM= 0.0220589707077938 ; MAM= 0.210833147073307"
[1] "iteration 6 ; RMSM= 0.0220333752919799 ; MAM= 0.209907295028066"
[1] "iteration 7 ; RMSM= 0.0216332457953376 ; MAM= 0.206767185820693"
[1] "iteration 8 ; RMSM= 0.0213847935159715 ; MAM= 0.202143138853316"
[1] "iteration 9 ; RMSM= 0.0209573813348291 ; MAM= 0.199078898287108"
[1] "iteration 10 ; RMSM= 0.0205096190727622 ; MAM= 0.194768566421505"
[1] "iteration 11 ; RMSM= 0.0202331763018993 ; MAM= 0.19172719464607"
[1] "iteration 12 ; RMSM= 0.0198674752653809 ; MAM= 0.187623752743794"
[1] "iteration 13 ; RMSM= 0.0194507722414294 ; MAM= 0.183721302646999"
[1] "iteration 14 ; RMSM= 0.0191790833548109 ; MAM= 0.18144417585227"
[1] "iteration 15 ; RMSM= 0.018809107300258 ; MAM= 0.177787703228428"
[1] "iteration 16 ; RMSM= 0.0181988293858517 ; MAM= 0.16202981750349"
[1] "iteration 17 ; RMSM= 0.0167322796828425 ; MAM= 0.142483166875718"
[1] "iteration 18 ; RMSM= 0.0151847684571389 ; MAM= 0.126731207201459"
[1] "iteration 19 ; RMSM= 0.0143285372012006 ; MAM= 0.11920519031782"
[1] "iteration 20 ; RMSM= 0.0135382925249031 ; MAM= 0.120393018380255"
[1] "iteration 21 ; RMSM= 0.0127889072923085 ; MAM= 0.120971467803161"
[1] "iteration 22 ; RMSM= 0.0107743228890709 ; MAM= 0.100390792286874"
[1] "iteration 23 ; RMSM= 0.00745590735327784 ; MAM= 0.0714343057516303"
[1] "iteration 24 ; RMSM= 0.00574950423829575 ; MAM= 0.0622131719247819"
[1] "iteration 25 ; RMSM= 0.00444497110351685 ; MAM= 0.0455792029465228"
[1] "iteration 26 ; RMSM= 0.00330256150421884 ; MAM= 0.033322742004"
[1] "iteration 27 ; RMSM= 0.00253955925041868 ; MAM= 0.026089596470043"
[1] "iteration 28 ; RMSM= 0.00236766344610807 ; MAM= 0.0158656238839016"
[1] "iteration 29 ; RMSM= 0.0000760927984409066 ; MAM= 0.000465483098671837"
[1] "stopped because of small RMSM= 0.00000000460175128297654"
day nr.: 144
[1] "iteration 1 ; RMSM= 0.0118783114957754 ; MAM= 0.167902223570255"
[1] "iteration 2 ; RMSM= 0.00425445387818911 ; MAM= 0.0594745904428678"
[1] "iteration 3 ; RMSM= 0.000885330460907957 ; MAM= 0.0110899517931335"
[1] "iteration 4 ; RMSM= 0.0000264844480167505 ; MAM= 0.000284216329993675"
[1] "stopped because of small RMSM= 0.0000000141268132791374"
day nr.: 145
[1] "iteration 1 ; RMSM= 0.001207715207448 ; MAM= 0.00863730607139746"
[1] "iteration 2 ; RMSM= 0.000247460692848719 ; MAM= 0.00169775218199105"
[1] "iteration 3 ; RMSM= 0.0000161282426029198 ; MAM= 0.000144650994306982"
[1] "stopped because of small RMSM= 0.0000000208388823456745"
day nr.: 146
[1] "iteration 1 ; RMSM= 0.000586534128951818 ; MAM= 0.00342718188866835"
[1] "iteration 2 ; RMSM= 0.0000737239274611012 ; MAM= 0.000569242739215998"
[1] "stopped because of small RMSM= 0.00000360275444584343"
day nr.: 147
[1] "iteration 1 ; RMSM= 0.000643765677997476 ; MAM= 0.00408678754473094"
[1] "iteration 2 ; RMSM= 0.000109422925450594 ; MAM= 0.000845323330929599"
[1] "iteration 3 ; RMSM= 0.0000156458314461991 ; MAM= 0.000141616676870395"
[1] "stopped because of small RMSM= 0.000000202412579556588"
day nr.: 148
[1] "iteration 1 ; RMSM= 0.000902116577262944 ; MAM= 0.00457880092688251"
[1] "iteration 2 ; RMSM= 0.000588847052022243 ; MAM= 0.00626044461364771"
[1] "iteration 3 ; RMSM= 0.000110550459688012 ; MAM= 0.00142312595599035"
[1] "iteration 4 ; RMSM= 0.0000107365762350864 ; MAM= 0.000126700116846866"
[1] "stopped because of small RMSM= 0.0000000102138159680339"
day nr.: 149
[1] "iteration 1 ; RMSM= 0.00119227588669708 ; MAM= 0.0094197936651684"
[1] "iteration 2 ; RMSM= 0.000903821567683713 ; MAM= 0.00694407772769815"
[1] "iteration 3 ; RMSM= 0.000700790981872334 ; MAM= 0.00555111232318028"
[1] "iteration 4 ; RMSM= 0.000516150822315424 ; MAM= 0.00358075916750366"
[1] "iteration 5 ; RMSM= 0.000397770684927298 ; MAM= 0.00298332572921865"
[1] "iteration 6 ; RMSM= 0.000292896630737783 ; MAM= 0.00150769206954317"
[1] "iteration 7 ; RMSM= 0.000228023187886962 ; MAM= 0.00153773791157504"
[1] "iteration 8 ; RMSM= 0.000190130249218982 ; MAM= 0.00179771013176298"
[1] "iteration 9 ; RMSM= 0.000124441906897455 ; MAM= 0.00138568362263802"
[1] "stopped because of small RMSM= 0.00000481644345151869"
day nr.: 150
[1] "iteration 1 ; RMSM= 0.000913521838880025 ; MAM= 0.00349595315284634"
[1] "iteration 2 ; RMSM= 0.000875233577494666 ; MAM= 0.00479653471247059"
[1] "iteration 3 ; RMSM= 0.000725669929990676 ; MAM= 0.00584094254598604"
[1] "iteration 4 ; RMSM= 0.000695035688529998 ; MAM= 0.00593655307575686"
[1] "iteration 5 ; RMSM= 0.00051540699983613 ; MAM= 0.00359673467996992"
[1] "iteration 6 ; RMSM= 0.000406714348112204 ; MAM= 0.00334847707225272"
[1] "iteration 7 ; RMSM= 0.000404953203848065 ; MAM= 0.00400814239720093"
[1] "iteration 8 ; RMSM= 0.000321257560239975 ; MAM= 0.00336021307318326"
[1] "iteration 9 ; RMSM= 0.000229282780779957 ; MAM= 0.0022124287681638"
[1] "iteration 10 ; RMSM= 0.000193292303001049 ; MAM= 0.00213118492740271"
[1] "iteration 11 ; RMSM= 0.000161528941608945 ; MAM= 0.00191502209769756"
[1] "iteration 12 ; RMSM= 0.000123570249109736 ; MAM= 0.00147208585934414"
[1] "iteration 13 ; RMSM= 0.000116167327729306 ; MAM= 0.00138499849807028"
[1] "iteration 14 ; RMSM= 0.000109383687823233 ; MAM= 0.00130688946298427"
[1] "iteration 15 ; RMSM= 0.0000806137646020617 ; MAM= 0.000872024640445349"
[1] "stopped because of small RMSM= 0.00000140510452339548"
day nr.: 151
[1] "iteration 1 ; RMSM= 0.000785522817798539 ; MAM= 0.00715569986498202"
[1] "iteration 2 ; RMSM= 0.000783560850717068 ; MAM= 0.00726709318595769"
[1] "iteration 3 ; RMSM= 0.000739306242407582 ; MAM= 0.00699519985062943"
[1] "iteration 4 ; RMSM= 0.000735704803853296 ; MAM= 0.00694525762635925"
[1] "iteration 5 ; RMSM= 0.000734987402584034 ; MAM= 0.00693851873014939"
[1] "iteration 6 ; RMSM= 0.000734911452594397 ; MAM= 0.00693896485330658"
[1] "iteration 7 ; RMSM= 0.000666453811166355 ; MAM= 0.00769015985745903"
[1] "iteration 8 ; RMSM= 0.000520204314589678 ; MAM= 0.00608370451911575"
[1] "iteration 9 ; RMSM= 0.000380116141648747 ; MAM= 0.00358339664889524"
[1] "iteration 10 ; RMSM= 0.000259039887334991 ; MAM= 0.00214612993885119"
[1] "iteration 11 ; RMSM= 0.0000393851049596368 ; MAM= 0.000485868972663506"
[1] "stopped because of small RMSM= 0.000000236257967255932"
day nr.: 152
[1] "iteration 1 ; RMSM= 0.000750617231611513 ; MAM= 0.00372554921776709"
[1] "iteration 2 ; RMSM= 0.000719096211823562 ; MAM= 0.00486668780943331"
[1] "iteration 3 ; RMSM= 0.000575621846949395 ; MAM= 0.0046932477272141"
[1] "iteration 4 ; RMSM= 0.000463552367438635 ; MAM= 0.00434680442789999"
[1] "iteration 5 ; RMSM= 0.000361048523946496 ; MAM= 0.00263815299065483"
[1] "iteration 6 ; RMSM= 0.000255895923132272 ; MAM= 0.00183810452291742"
[1] "iteration 7 ; RMSM= 0.0000304930150306496 ; MAM= 0.000350168343682546"
[1] "stopped because of small RMSM= 0.000000337748726174802"
day nr.: 153
[1] "iteration 1 ; RMSM= 0.000350395058660713 ; MAM= 0.00228932913602282"
[1] "iteration 2 ; RMSM= 0.000263447722934419 ; MAM= 0.00146986997164998"
[1] "iteration 3 ; RMSM= 0.000154440452000909 ; MAM= 0.0016678283681146"
[1] "stopped because of small RMSM= 0.0000093894719689717"
day nr.: 154
[1] "iteration 1 ; RMSM= 0.0508818939619092 ; MAM= 0.757864155006502"
[1] "iteration 2 ; RMSM= 0.0500223340018617 ; MAM= 0.745057746390497"
[1] "iteration 3 ; RMSM= 0.0489875342767247 ; MAM= 0.729633193122589"
[1] "iteration 4 ; RMSM= 0.0482298726546274 ; MAM= 0.718320643407103"
[1] "iteration 5 ; RMSM= 0.0469622270672982 ; MAM= 0.699361128724963"
[1] "iteration 6 ; RMSM= 0.0465487540337261 ; MAM= 0.693287061255787"
[1] "iteration 7 ; RMSM= 0.0464385255220307 ; MAM= 0.691655213893405"
[1] "iteration 8 ; RMSM= 0.0459933648212895 ; MAM= 0.681237643608991"
[1] "iteration 9 ; RMSM= 0.0458210011173342 ; MAM= 0.678843863747698"
[1] "iteration 10 ; RMSM= 0.0456297148317436 ; MAM= 0.676002879254304"
[1] "iteration 11 ; RMSM= 0.0454647285030637 ; MAM= 0.673737754720944"
[1] "iteration 12 ; RMSM= 0.0453021137257887 ; MAM= 0.671225396149532"
[1] "iteration 13 ; RMSM= 0.0451412489084841 ; MAM= 0.668920885959854"
[1] "iteration 14 ; RMSM= 0.0449725164364424 ; MAM= 0.666392135920859"
[1] "iteration 15 ; RMSM= 0.0448063904648188 ; MAM= 0.663839757381457"
[1] "iteration 16 ; RMSM= 0.0446454141161294 ; MAM= 0.661380201032446"
[1] "iteration 17 ; RMSM= 0.0444926188481911 ; MAM= 0.659024913412378"
[1] "iteration 18 ; RMSM= 0.0441231722374785 ; MAM= 0.647649276717309"
[1] "iteration 19 ; RMSM= 0.0438767073868042 ; MAM= 0.644464633038193"
[1] "iteration 20 ; RMSM= 0.0437055983657003 ; MAM= 0.642107912949262"
[1] "iteration 21 ; RMSM= 0.0435211354936256 ; MAM= 0.639447600045331"
[1] "iteration 22 ; RMSM= 0.0433628916265997 ; MAM= 0.637073518590765"
[1] "iteration 23 ; RMSM= 0.0428239048948062 ; MAM= 0.624852397702469"
[1] "iteration 24 ; RMSM= 0.0422990835559792 ; MAM= 0.615935923140838"
[1] "iteration 25 ; RMSM= 0.0418225735828516 ; MAM= 0.604058754013822"
[1] "iteration 26 ; RMSM= 0.0415641540213937 ; MAM= 0.593402833370711"
[1] "iteration 27 ; RMSM= 0.0409420671730674 ; MAM= 0.584670889816248"
[1] "iteration 28 ; RMSM= 0.0400286611677649 ; MAM= 0.573943037665685"
[1] "iteration 29 ; RMSM= 0.0393931489107301 ; MAM= 0.565600649623245"
[1] "iteration 30 ; RMSM= 0.0388119260150261 ; MAM= 0.557039809069216"
[1] "iteration 31 ; RMSM= 0.0381312003364749 ; MAM= 0.548834162485671"
[1] "iteration 32 ; RMSM= 0.0375826202076221 ; MAM= 0.539972878919999"
[1] "iteration 33 ; RMSM= 0.0369440473553603 ; MAM= 0.532258381054769"
[1] "iteration 34 ; RMSM= 0.0361891801761996 ; MAM= 0.521171059799005"
[1] "iteration 35 ; RMSM= 0.0355885388372281 ; MAM= 0.513711741664268"
[1] "iteration 36 ; RMSM= 0.0350441743542886 ; MAM= 0.504703521401237"
[1] "iteration 37 ; RMSM= 0.034463626630765 ; MAM= 0.497438527229129"
[1] "iteration 38 ; RMSM= 0.0339602048868598 ; MAM= 0.490172225240796"
[1] "iteration 39 ; RMSM= 0.033348006567438 ; MAM= 0.481064562092785"
[1] "iteration 40 ; RMSM= 0.03282598666948 ; MAM= 0.47405412831079"
[1] "iteration 41 ; RMSM= 0.0323979343205911 ; MAM= 0.46707373686146"
[1] "iteration 42 ; RMSM= 0.031802204406727 ; MAM= 0.458874139315281"
[1] "iteration 43 ; RMSM= 0.0312228786275634 ; MAM= 0.45052346829971"
[1] "iteration 44 ; RMSM= 0.0307795070635421 ; MAM= 0.443863531653487"
[1] "iteration 45 ; RMSM= 0.0302587156266469 ; MAM= 0.436569404039519"
[1] "iteration 46 ; RMSM= 0.0298228674266782 ; MAM= 0.430069403772383"
[1] "iteration 47 ; RMSM= 0.0293022755171071 ; MAM= 0.422374399844016"
[1] "iteration 48 ; RMSM= 0.0288420185342847 ; MAM= 0.416074666746489"
[1] "iteration 49 ; RMSM= 0.0283670306431731 ; MAM= 0.409175850341893"
[1] "iteration 50 ; RMSM= 0.0278917396507715 ; MAM= 0.384486745884718"
[1] "iteration 51 ; RMSM= 0.0277679412138551 ; MAM= 0.363663391123137"
[1] "iteration 52 ; RMSM= 0.0248658620289867 ; MAM= 0.34066749755981"
[1] "iteration 53 ; RMSM= 0.0228023343035071 ; MAM= 0.318388063155704"
[1] "iteration 54 ; RMSM= 0.0224142400981912 ; MAM= 0.313543466664952"
[1] "iteration 55 ; RMSM= 0.0212091060644436 ; MAM= 0.295825743064297"
[1] "iteration 56 ; RMSM= 0.0197601503185173 ; MAM= 0.274212851146886"
[1] "iteration 57 ; RMSM= 0.0186169142267194 ; MAM= 0.258450152255065"
[1] "iteration 58 ; RMSM= 0.0175765824689275 ; MAM= 0.243239809385887"
[1] "iteration 59 ; RMSM= 0.0162093013087114 ; MAM= 0.225608730449395"
[1] "iteration 60 ; RMSM= 0.0152324722236851 ; MAM= 0.212324729577355"
[1] "iteration 61 ; RMSM= 0.0142910213932841 ; MAM= 0.199643898675854"
[1] "iteration 62 ; RMSM= 0.0132515008910795 ; MAM= 0.185139321946221"
[1] "iteration 63 ; RMSM= 0.0125238944559173 ; MAM= 0.174122794513403"
[1] "iteration 64 ; RMSM= 0.0117462015799149 ; MAM= 0.163663275543655"
[1] "iteration 65 ; RMSM= 0.0110442013220659 ; MAM= 0.153765279711354"
[1] "iteration 66 ; RMSM= 0.0105921219917639 ; MAM= 0.108792953063659"
[1] "iteration 67 ; RMSM= 0.00727709851345172 ; MAM= 0.0832340076645283"
[1] "iteration 68 ; RMSM= 0.00540133072978354 ; MAM= 0.0615862302853956"
[1] "iteration 69 ; RMSM= 0.0040892459661292 ; MAM= 0.0469463358317082"
[1] "iteration 70 ; RMSM= 0.00304705539814658 ; MAM= 0.0355607891994554"
[1] "iteration 71 ; RMSM= 0.00226424216809602 ; MAM= 0.0266871184591692"
[1] "iteration 72 ; RMSM= 0.0019075129510388 ; MAM= 0.00884275553524477"
[1] "iteration 73 ; RMSM= 0.0000964703140284629 ; MAM= 0.00119628008356867"
[1] "stopped because of small RMSM= 0.00000102029878030696"
day nr.: 155
[1] "iteration 1 ; RMSM= 0.00693482089402016 ; MAM= 0.079621258728569"
[1] "iteration 2 ; RMSM= 0.00678140577175682 ; MAM= 0.0768828359469653"
[1] "iteration 3 ; RMSM= 0.00670041931430327 ; MAM= 0.0757475401789987"
[1] "iteration 4 ; RMSM= 0.0066000206178492 ; MAM= 0.074619308855213"
[1] "iteration 5 ; RMSM= 0.00650240067784715 ; MAM= 0.0735022095774905"
[1] "iteration 6 ; RMSM= 0.00641025642781682 ; MAM= 0.0723954791289071"
[1] "iteration 7 ; RMSM= 0.00631675077977382 ; MAM= 0.0715709605583665"
[1] "iteration 8 ; RMSM= 0.00622414547317994 ; MAM= 0.0704608265556003"
[1] "iteration 9 ; RMSM= 0.00613240446755696 ; MAM= 0.0693649674637301"
[1] "iteration 10 ; RMSM= 0.00609978148756092 ; MAM= 0.0644406594110323"
[1] "iteration 11 ; RMSM= 0.00602146454840197 ; MAM= 0.0640368146906235"
[1] "iteration 12 ; RMSM= 0.0058782564028861 ; MAM= 0.0612858422117268"
[1] "iteration 13 ; RMSM= 0.00556905709441557 ; MAM= 0.0573759328000667"
[1] "iteration 14 ; RMSM= 0.00512121511409856 ; MAM= 0.0496833188750045"
[1] "iteration 15 ; RMSM= 0.00496110818268771 ; MAM= 0.046965143663562"
[1] "iteration 16 ; RMSM= 0.0046977931537167 ; MAM= 0.0470170249416291"
[1] "iteration 17 ; RMSM= 0.00458764477868918 ; MAM= 0.0439820079625662"
[1] "iteration 18 ; RMSM= 0.00430168278978204 ; MAM= 0.0410301747696757"
[1] "iteration 19 ; RMSM= 0.00397086104228797 ; MAM= 0.03643963508076"
[1] "iteration 20 ; RMSM= 0.00379713015592072 ; MAM= 0.0343147729981137"
[1] "iteration 21 ; RMSM= 0.00357676663209446 ; MAM= 0.0322785569465966"
[1] "iteration 22 ; RMSM= 0.00331511573528085 ; MAM= 0.0304070276049019"
[1] "iteration 23 ; RMSM= 0.00322851224634323 ; MAM= 0.0297908316139678"
[1] "iteration 24 ; RMSM= 0.00302419962515867 ; MAM= 0.0278705712964045"
[1] "iteration 25 ; RMSM= 0.00284331587915636 ; MAM= 0.0260498355864135"
[1] "iteration 26 ; RMSM= 0.00262169094634501 ; MAM= 0.023284979332033"
[1] "iteration 27 ; RMSM= 0.00251599967601724 ; MAM= 0.0218800725223746"
[1] "iteration 28 ; RMSM= 0.00236129002141721 ; MAM= 0.0205520292872062"
[1] "iteration 29 ; RMSM= 0.0022138213643353 ; MAM= 0.0192973759069478"
[1] "iteration 30 ; RMSM= 0.00208597276548296 ; MAM= 0.0181145665201681"
[1] "iteration 31 ; RMSM= 0.00197627289790278 ; MAM= 0.0170007970012075"
[1] "iteration 32 ; RMSM= 0.00181742212697773 ; MAM= 0.015952825036365"
[1] "iteration 33 ; RMSM= 0.00173462838967268 ; MAM= 0.014967164758779"
[1] "iteration 34 ; RMSM= 0.00164002571840865 ; MAM= 0.0144777567254143"
[1] "iteration 35 ; RMSM= 0.00155095247426355 ; MAM= 0.0135688237796094"
[1] "iteration 36 ; RMSM= 0.00146305582492168 ; MAM= 0.0130385201988291"
[1] "iteration 37 ; RMSM= 0.00142740636077939 ; MAM= 0.00853144617600476"
[1] "iteration 38 ; RMSM= 0.00113778650491048 ; MAM= 0.00650054156241786"
[1] "iteration 39 ; RMSM= 0.000863964807057533 ; MAM= 0.00491971706133232"
[1] "iteration 40 ; RMSM= 0.000695804219490821 ; MAM= 0.00371017975796917"
[1] "iteration 41 ; RMSM= 0.000469332924333298 ; MAM= 0.0027925086657572"
[1] "iteration 42 ; RMSM= 0.000358797986908962 ; MAM= 0.0020992932242585"
[1] "iteration 43 ; RMSM= 0.00027745162745262 ; MAM= 0.00157704859678764"
[1] "iteration 44 ; RMSM= 0.000215539047029125 ; MAM= 0.00118412245104738"
[1] "iteration 45 ; RMSM= 0.000166283136268972 ; MAM= 0.000888823047115572"
[1] "iteration 46 ; RMSM= 0.000125244589835149 ; MAM= 0.000667003230819474"
[1] "iteration 47 ; RMSM= 0.000120202559247503 ; MAM= 0.000604091737971768"
[1] "stopped because of small RMSM= 0.0000000173862477824128"
day nr.: 156
[1] "iteration 1 ; RMSM= 0.0310097414413448 ; MAM= 0.2914384187024"
[1] "iteration 2 ; RMSM= 0.0104804733391486 ; MAM= 0.101124036921275"
[1] "iteration 3 ; RMSM= 0.00217300473681587 ; MAM= 0.0245595180451223"
[1] "iteration 4 ; RMSM= 0.000087935558626338 ; MAM= 0.000969436833495582"
[1] "stopped because of small RMSM= 0.0000000737082078866378"
day nr.: 157
[1] "iteration 1 ; RMSM= 0.00568240281302171 ; MAM= 0.0781759081806796"
[1] "iteration 2 ; RMSM= 0.00371483797680724 ; MAM= 0.0188954152561691"
[1] "iteration 3 ; RMSM= 0.000303979084634823 ; MAM= 0.00220318577136283"
[1] "stopped because of small RMSM= 0.00000454786540622365"
day nr.: 158
[1] "iteration 1 ; RMSM= 0.00718315171710381 ; MAM= 0.102169334304605"
[1] "iteration 2 ; RMSM= 0.00111300595309451 ; MAM= 0.0158728429088527"
[1] "iteration 3 ; RMSM= 0.0000190224613822397 ; MAM= 0.000256833884942756"
[1] "stopped because of small RMSM= 0.00000000332508741735287"
day nr.: 159
[1] "iteration 1 ; RMSM= 0.000788796351584323 ; MAM= 0.00947963833885546"
[1] "iteration 2 ; RMSM= 0.0000616762342596091 ; MAM= 0.000816639150665025"
[1] "stopped because of small RMSM= 0.00000024140044647929"
day nr.: 160
[1] "iteration 1 ; RMSM= 0.000195240218226917 ; MAM= 0.00103734922625346"
[1] "stopped because of small RMSM= 0.00000544654512300499"
day nr.: 161
[1] "iteration 1 ; RMSM= 0.00884010999708547 ; MAM= 0.131481786509428"
[1] "iteration 2 ; RMSM= 0.00702926645407606 ; MAM= 0.102047115310213"
[1] "iteration 3 ; RMSM= 0.000932851574922212 ; MAM= 0.0126480378343778"
[1] "stopped because of small RMSM= 0.00000910544213890519"
day nr.: 162
[1] "iteration 1 ; RMSM= 0.00947175872061551 ; MAM= 0.140521134602891"
[1] "iteration 2 ; RMSM= 0.00585094008132826 ; MAM= 0.0770963253043292"
[1] "iteration 3 ; RMSM= 0.000714931015780659 ; MAM= 0.00657864436096989"
[1] "iteration 4 ; RMSM= 0.0000128649948239943 ; MAM= 0.00014669249696625"
[1] "stopped because of small RMSM= 0.000000000091952301126998"
day nr.: 163
[1] "iteration 1 ; RMSM= 0.00073248635606419 ; MAM= 0.00551207346914689"
[1] "iteration 2 ; RMSM= 0.0000133283309459636 ; MAM= 0.000128924671855657"
[1] "stopped because of small RMSM= 0.000000000431667842163378"
day nr.: 164
[1] "iteration 1 ; RMSM= 0.0366808671195872 ; MAM= 0.478024100950425"
[1] "iteration 2 ; RMSM= 0.0286284858902545 ; MAM= 0.349609533400483"
[1] "iteration 3 ; RMSM= 0.0227644613697726 ; MAM= 0.264834292200052"
[1] "iteration 4 ; RMSM= 0.0171693780733093 ; MAM= 0.192195514395809"
[1] "iteration 5 ; RMSM= 0.0135061699641831 ; MAM= 0.146168522274295"
[1] "iteration 6 ; RMSM= 0.0100712863801874 ; MAM= 0.110366851061903"
[1] "iteration 7 ; RMSM= 0.00852847629313906 ; MAM= 0.0599671851240904"
[1] "iteration 8 ; RMSM= 0.000443447064425746 ; MAM= 0.00276163180388456"
[1] "stopped because of small RMSM= 0.0000000335598005154624"
day nr.: 165
[1] "iteration 1 ; RMSM= 0.0459432447584807 ; MAM= 0.663489574037176"
[1] "iteration 2 ; RMSM= 0.0134013007855495 ; MAM= 0.191773007041488"
[1] "iteration 3 ; RMSM= 0.00121510281462626 ; MAM= 0.0168458629819681"
[1] "stopped because of small RMSM= 0.00000629027140506458"
day nr.: 166
[1] "iteration 1 ; RMSM= 0.00128129972364793 ; MAM= 0.0114243587274515"
[1] "iteration 2 ; RMSM= 0.0000820056706891092 ; MAM= 0.000589173703850809"
[1] "stopped because of small RMSM= 0.000000155164800670354"
day nr.: 167
[1] "iteration 1 ; RMSM= 0.0298780717797025 ; MAM= 0.444433628430877"
[1] "iteration 2 ; RMSM= 0.0297753992284275 ; MAM= 0.435046549712696"
[1] "iteration 3 ; RMSM= 0.00435989614981675 ; MAM= 0.0592156406423971"
[1] "iteration 4 ; RMSM= 0.000111658981006877 ; MAM= 0.00134543323125658"
[1] "stopped because of small RMSM= 0.000000022707396096194"
day nr.: 168
[1] "iteration 1 ; RMSM= 0.00241113031329165 ; MAM= 0.0332443441771415"
[1] "iteration 2 ; RMSM= 0.000023410508825448 ; MAM= 0.000321642517373377"
[1] "stopped because of small RMSM= 0.00000000143665401939815"
day nr.: 169
[1] "iteration 1 ; RMSM= 0.00805410165654433 ; MAM= 0.116736034187105"
[1] "iteration 2 ; RMSM= 0.00115428133804264 ; MAM= 0.0169521868742846"
[1] "iteration 3 ; RMSM= 0.0000181743801835646 ; MAM= 0.000256815536784329"
[1] "stopped because of small RMSM= 0.00000000223290424445602"
day nr.: 170
[1] "iteration 1 ; RMSM= 0.00035558991652263 ; MAM= 0.00298636756757292"
[1] "iteration 2 ; RMSM= 0.000010607228543303 ; MAM= 0.0000867104096548502"
[1] "stopped because of small RMSM= 0.00000000513828512936226"
day nr.: 171
[1] "iteration 1 ; RMSM= 0.000366414847455451 ; MAM= 0.00478504351316556"
[1] "stopped because of small RMSM= 0.00000594307648943598"
day nr.: 172
[1] "iteration 1 ; RMSM= 0.000599075597723549 ; MAM= 0.00821651715150426"
[1] "iteration 2 ; RMSM= 0.0000370208636619409 ; MAM= 0.000488179570235442"
[1] "stopped because of small RMSM= 0.000000103225479547273"
day nr.: 173
[1] "iteration 1 ; RMSM= 0.000555687774294877 ; MAM= 0.0038853071727265"
[1] "iteration 2 ; RMSM= 0.0000620998861000901 ; MAM= 0.000440600828443223"
[1] "stopped because of small RMSM= 0.000000126306703844227"
day nr.: 174
[1] "iteration 1 ; RMSM= 0.000902872604761643 ; MAM= 0.005234990710444"
[1] "iteration 2 ; RMSM= 0.000101009534319855 ; MAM= 0.000671102989982804"
[1] "stopped because of small RMSM= 0.000000771945468787192"
day nr.: 175
[1] "iteration 1 ; RMSM= 0.000866539770011078 ; MAM= 0.00493841638792718"
[1] "iteration 2 ; RMSM= 0.000117973618646697 ; MAM= 0.000687733515997296"
[1] "stopped because of small RMSM= 0.000000543571216102459"
day nr.: 176
[1] "iteration 1 ; RMSM= 0.0102429455335309 ; MAM= 0.145654490141439"
[1] "iteration 2 ; RMSM= 0.00700964921521779 ; MAM= 0.0969903781649964"
[1] "iteration 3 ; RMSM= 0.00661571432766635 ; MAM= 0.0762263985139045"
[1] "iteration 4 ; RMSM= 0.00168308506664497 ; MAM= 0.0145282417920589"
[1] "iteration 5 ; RMSM= 0.0000341163312301496 ; MAM= 0.00028029193111162"
[1] "stopped because of small RMSM= 0.00000000260127276387325"
day nr.: 177
[1] "iteration 1 ; RMSM= 0.00294553373487913 ; MAM= 0.0408515033016974"
[1] "iteration 2 ; RMSM= 0.000106490009594287 ; MAM= 0.00137276062477926"
[1] "stopped because of small RMSM= 0.0000000991654752872895"
day nr.: 178
[1] "iteration 1 ; RMSM= 0.00140227128579759 ; MAM= 0.0189198206530007"
[1] "iteration 2 ; RMSM= 0.0000607059520012787 ; MAM= 0.000803009269198164"
[1] "stopped because of small RMSM= 0.0000000756581890133288"
day nr.: 179
[1] "iteration 1 ; RMSM= 0.00106865822794875 ; MAM= 0.0156963096167691"
[1] "iteration 2 ; RMSM= 0.0000914749237890393 ; MAM= 0.00134575275868196"
[1] "stopped because of small RMSM= 0.000000496090802299802"
day nr.: 180
[1] "iteration 1 ; RMSM= 0.0559473084253771 ; MAM= 0.833559774102235"
[1] "iteration 2 ; RMSM= 0.00991992685429065 ; MAM= 0.138883242860988"
[1] "iteration 3 ; RMSM= 0.00933038754743452 ; MAM= 0.129103365197225"
[1] "iteration 4 ; RMSM= 0.00890577444839769 ; MAM= 0.123825082635496"
[1] "iteration 5 ; RMSM= 0.00856876138013873 ; MAM= 0.0955236272719401"
[1] "iteration 6 ; RMSM= 0.00675349209286398 ; MAM= 0.0744425593652136"
[1] "iteration 7 ; RMSM= 0.00564095458041491 ; MAM= 0.0571825048736012"
[1] "iteration 8 ; RMSM= 0.0043608439074886 ; MAM= 0.0454595733261645"
[1] "iteration 9 ; RMSM= 0.00333806833397136 ; MAM= 0.0348471487723653"
[1] "iteration 10 ; RMSM= 0.00246610907448327 ; MAM= 0.0272167790658655"
[1] "iteration 11 ; RMSM= 0.00195516741165994 ; MAM= 0.0209316050234067"
[1] "iteration 12 ; RMSM= 0.00146673136973702 ; MAM= 0.0159589061290306"
[1] "iteration 13 ; RMSM= 0.00141811142583506 ; MAM= 0.0107357942838169"
[1] "stopped because of small RMSM= 0.00000997188151644126"
day nr.: 181
[1] "iteration 1 ; RMSM= 0.0213371662491624 ; MAM= 0.255089884141682"
[1] "iteration 2 ; RMSM= 0.0161793432594993 ; MAM= 0.198880385325078"
[1] "iteration 3 ; RMSM= 0.0125361871267678 ; MAM= 0.149616524630777"
[1] "iteration 4 ; RMSM= 0.0100271639899027 ; MAM= 0.117296131174782"
[1] "iteration 5 ; RMSM= 0.0075756237428942 ; MAM= 0.0905070258443368"
[1] "iteration 6 ; RMSM= 0.00745115678780508 ; MAM= 0.0503090209145973"
[1] "iteration 7 ; RMSM= 0.000391298633277508 ; MAM= 0.00221812806844013"
[1] "stopped because of small RMSM= 0.0000000413904675963074"
day nr.: 182
[1] "iteration 1 ; RMSM= 0.0150206025431347 ; MAM= 0.204484306789956"
[1] "iteration 2 ; RMSM= 0.00126338255764543 ; MAM= 0.0161576516240825"
[1] "stopped because of small RMSM= 0.00000402265773403829"
day nr.: 183
[1] "iteration 1 ; RMSM= 0.00928318164262648 ; MAM= 0.133110553400733"
[1] "iteration 2 ; RMSM= 0.00136333323385816 ; MAM= 0.0199767699559452"
[1] "iteration 3 ; RMSM= 0.0000227670625508004 ; MAM= 0.000326180006360987"
[1] "stopped because of small RMSM= 0.00000000323004269338993"
day nr.: 184
[1] "iteration 1 ; RMSM= 0.0371621430674296 ; MAM= 0.553253104480902"
[1] "iteration 2 ; RMSM= 0.0257920721534757 ; MAM= 0.38377048916026"
[1] "iteration 3 ; RMSM= 0.0131768203816223 ; MAM= 0.159549210677712"
[1] "iteration 4 ; RMSM= 0.00133131557743669 ; MAM= 0.0118816136528611"
[1] "iteration 5 ; RMSM= 0.0000137110461284002 ; MAM= 0.00013195352008144"
[1] "stopped because of small RMSM= 0.0000000000798158984019651"
day nr.: 185
[1] "iteration 1 ; RMSM= 0.00621673935076787 ; MAM= 0.0882172952880034"
[1] "iteration 2 ; RMSM= 0.000145297239702201 ; MAM= 0.00198716610658922"
[1] "stopped because of small RMSM= 0.0000000517780627959071"
day nr.: 186
[1] "iteration 1 ; RMSM= 0.00121486681829451 ; MAM= 0.0156606923166808"
[1] "iteration 2 ; RMSM= 0.0000358802760471577 ; MAM= 0.00048061269949784"
[1] "stopped because of small RMSM= 0.00000000207528274637223"
day nr.: 187
[1] "iteration 1 ; RMSM= 0.000426103252604514 ; MAM= 0.00459643121046546"
[1] "stopped because of small RMSM= 0.00000138117336467326"
day nr.: 188
[1] "iteration 1 ; RMSM= 0.0305406081479264 ; MAM= 0.451812549288658"
[1] "iteration 2 ; RMSM= 0.00472963109920119 ; MAM= 0.0635013561781227"
[1] "iteration 3 ; RMSM= 0.0000894541301037808 ; MAM= 0.000866229618018965"
[1] "stopped because of small RMSM= 0.000000164414071833768"
day nr.: 189
[1] "iteration 1 ; RMSM= 0.0050318309671213 ; MAM= 0.0561709230339511"
[1] "iteration 2 ; RMSM= 0.000448272175374427 ; MAM= 0.00307378875154563"
[1] "stopped because of small RMSM= 0.00000863518536782725"
day nr.: 190
[1] "iteration 1 ; RMSM= 0.0156969148676112 ; MAM= 0.226970498713178"
[1] "iteration 2 ; RMSM= 0.00189885955627315 ; MAM= 0.0276570196756623"
[1] "iteration 3 ; RMSM= 0.0000183566465657598 ; MAM= 0.000258992140184724"
[1] "stopped because of small RMSM= 0.000000000911345152266654"
day nr.: 191
[1] "iteration 1 ; RMSM= 0.00207894377340432 ; MAM= 0.0285284221744607"
[1] "iteration 2 ; RMSM= 0.00017159854160938 ; MAM= 0.00221119000852998"
[1] "stopped because of small RMSM= 0.000000868871203109686"
day nr.: 192
[1] "iteration 1 ; RMSM= 0.000517895801479739 ; MAM= 0.00378207030541084"
[1] "iteration 2 ; RMSM= 0.0000112435279639017 ; MAM= 0.0000945089560742893"
[1] "stopped because of small RMSM= 0.0000000079697923583525"
day nr.: 193
[1] "iteration 1 ; RMSM= 0.000382052312252907 ; MAM= 0.0015719951217022"
[1] "iteration 2 ; RMSM= 0.0000170539911518467 ; MAM= 0.000123081287347969"
[1] "stopped because of small RMSM= 0.00000000365254794742245"
day nr.: 194
[1] "iteration 1 ; RMSM= 0.0137364213202791 ; MAM= 0.203509385947778"
[1] "iteration 2 ; RMSM= 0.00758597412569679 ; MAM= 0.112077008794303"
[1] "iteration 3 ; RMSM= 0.000667189294247227 ; MAM= 0.00935852186904218"
[1] "stopped because of small RMSM= 0.00000185133647668546"
day nr.: 195
[1] "iteration 1 ; RMSM= 0.00276816873162496 ; MAM= 0.0410188783278895"
[1] "iteration 2 ; RMSM= 0.000105466837356353 ; MAM= 0.00149922685629442"
[1] "stopped because of small RMSM= 0.0000000880318006935239"
day nr.: 196
[1] "iteration 1 ; RMSM= 0.0188980144885109 ; MAM= 0.281397597508061"
[1] "iteration 2 ; RMSM= 0.0144244210182992 ; MAM= 0.213244449117063"
[1] "iteration 3 ; RMSM= 0.00135575998721758 ; MAM= 0.0186474470383044"
[1] "stopped because of small RMSM= 0.00000753860235997357"
day nr.: 197
[1] "iteration 1 ; RMSM= 0.00893038213051676 ; MAM= 0.131135795877075"
[1] "iteration 2 ; RMSM= 0.000906355913419014 ; MAM= 0.0131884256937781"
[1] "stopped because of small RMSM= 0.00000579808624594024"
day nr.: 198
[1] "iteration 1 ; RMSM= 0.000491003533927336 ; MAM= 0.00345088913725398"
[1] "iteration 2 ; RMSM= 0.0000166959697860213 ; MAM= 0.000150232543305119"
[1] "stopped because of small RMSM= 0.0000000120903459233088"
day nr.: 199
[1] "iteration 1 ; RMSM= 0.000181288742025261 ; MAM= 0.00101721432016946"
[1] "stopped because of small RMSM= 0.0000076354988981625"
day nr.: 200
[1] "iteration 1 ; RMSM= 0.000339980566948811 ; MAM= 0.00228034672588089"
[1] "iteration 2 ; RMSM= 0.0000152672249688587 ; MAM= 0.0000890035536413591"
[1] "stopped because of small RMSM= 0.000000002954438735894"
day nr.: 201
[1] "iteration 1 ; RMSM= 0.000427711749809744 ; MAM= 0.00269559445374135"
[1] "iteration 2 ; RMSM= 0.0000147452224431319 ; MAM= 0.000104245352758166"
[1] "stopped because of small RMSM= 0.00000000100320669903969"
day nr.: 202
[1] "iteration 1 ; RMSM= 0.000366174632329878 ; MAM= 0.00201076218226767"
[1] "stopped because of small RMSM= 0.0000040833792862168"
day nr.: 203
[1] "iteration 1 ; RMSM= 0.000521847538137023 ; MAM= 0.00338223147141006"
[1] "iteration 2 ; RMSM= 0.0000167775131852889 ; MAM= 0.000115389491981525"
[1] "stopped because of small RMSM= 0.000000000986453708797568"
day nr.: 204
[1] "iteration 1 ; RMSM= 0.0173159649117574 ; MAM= 0.256839901413743"
[1] "iteration 2 ; RMSM= 0.0120053414240673 ; MAM= 0.177412579076464"
[1] "iteration 3 ; RMSM= 0.00692184937925453 ; MAM= 0.0720684407685036"
[1] "iteration 4 ; RMSM= 0.000909219335232044 ; MAM= 0.00893956933084067"
[1] "iteration 5 ; RMSM= 0.0000601546343598549 ; MAM= 0.000750403474501005"
[1] "stopped because of small RMSM= 0.00000000222080315516475"
day nr.: 205
[1] "iteration 1 ; RMSM= 0.00297328784362155 ; MAM= 0.0265136844243267"
[1] "iteration 2 ; RMSM= 0.000142569627338043 ; MAM= 0.001621883713258"
[1] "stopped because of small RMSM= 0.0000012495865230732"
day nr.: 206
[1] "iteration 1 ; RMSM= 0.000420282110986784 ; MAM= 0.00330081053422263"
[1] "iteration 2 ; RMSM= 0.0000142383775052174 ; MAM= 0.000156435696784168"
[1] "stopped because of small RMSM= 0.0000000000548385086140268"
day nr.: 207
[1] "iteration 1 ; RMSM= 0.0031904621974196 ; MAM= 0.0468008779878265"
[1] "iteration 2 ; RMSM= 0.000199097444243708 ; MAM= 0.00276898787451429"
[1] "stopped because of small RMSM= 0.000000248614985698209"
day nr.: 208
[1] "iteration 1 ; RMSM= 0.00117284280306958 ; MAM= 0.017250941537094"
[1] "iteration 2 ; RMSM= 0.0000431796886126057 ; MAM= 0.000611729180542619"
[1] "stopped because of small RMSM= 0.0000000395100586647929"
day nr.: 209
[1] "iteration 1 ; RMSM= 0.00014438852263441 ; MAM= 0.000937011914452548"
[1] "stopped because of small RMSM= 0.00000364869401287001"
day nr.: 210
[1] "iteration 1 ; RMSM= 0.000364409585032337 ; MAM= 0.00232776699836332"
[1] "iteration 2 ; RMSM= 0.0000265055518588941 ; MAM= 0.000169392283594041"
[1] "stopped because of small RMSM= 0.00000000579935147743993"
day nr.: 211
[1] "iteration 1 ; RMSM= 0.000626018827899476 ; MAM= 0.00369557435131121"
[1] "iteration 2 ; RMSM= 0.0000499747174049401 ; MAM= 0.000289464384002763"
[1] "stopped because of small RMSM= 0.00000012914772615661"
day nr.: 212
[1] "iteration 1 ; RMSM= 0.000466636200599093 ; MAM= 0.00500610944358043"
[1] "stopped because of small RMSM= 0.0000024825210749069"
day nr.: 213
[1] "iteration 1 ; RMSM= 0.000887521069542615 ; MAM= 0.0114607188429047"
[1] "iteration 2 ; RMSM= 0.0000315278542034348 ; MAM= 0.000436103221569838"
[1] "stopped because of small RMSM= 0.0000000234194397241733"
day nr.: 214
[1] "iteration 1 ; RMSM= 0.000217807996074756 ; MAM= 0.00109325758543977"
[1] "stopped because of small RMSM= 0.000000457454727275687"
day nr.: 215
[1] "iteration 1 ; RMSM= 0.000640324289776116 ; MAM= 0.00347573587319469"
[1] "iteration 2 ; RMSM= 0.0000317634644799087 ; MAM= 0.000224589092616363"
[1] "stopped because of small RMSM= 0.0000000397481736241514"
day nr.: 216
[1] "iteration 1 ; RMSM= 0.000904162020740761 ; MAM= 0.00468308317852059"
[1] "iteration 2 ; RMSM= 0.0000697257378500072 ; MAM= 0.000484454975381861"
[1] "stopped because of small RMSM= 0.0000000249862893280471"
day nr.: 217
[1] "iteration 1 ; RMSM= 0.000847529632789622 ; MAM= 0.0042468704207131"
[1] "iteration 2 ; RMSM= 0.0000746986277260752 ; MAM= 0.000487015531256442"
[1] "stopped because of small RMSM= 0.0000000509372632353836"
day nr.: 218
[1] "iteration 1 ; RMSM= 0.000902140987081853 ; MAM= 0.00424925392652928"
[1] "iteration 2 ; RMSM= 0.0000920989286124954 ; MAM= 0.000652966801107119"
[1] "stopped because of small RMSM= 0.00000027706941134574"
day nr.: 219
[1] "iteration 1 ; RMSM= 0.000976037509571947 ; MAM= 0.00436913607634107"
[1] "iteration 2 ; RMSM= 0.000124583800800423 ; MAM= 0.000928312711697958"
[1] "stopped because of small RMSM= 0.00000126946214372669"
day nr.: 220
[1] "iteration 1 ; RMSM= 0.000950305853238975 ; MAM= 0.00354092311732418"
[1] "iteration 2 ; RMSM= 0.000104924296762528 ; MAM= 0.000523072497117754"
[1] "stopped because of small RMSM= 0.00000244235115312125"
day nr.: 221
[1] "iteration 1 ; RMSM= 0.000818925388502404 ; MAM= 0.00276683489704746"
[1] "iteration 2 ; RMSM= 0.000112985270254613 ; MAM= 0.000575725987185744"
[1] "stopped because of small RMSM= 0.00000495020933864549"
day nr.: 222
[1] "iteration 1 ; RMSM= 0.000745890920097489 ; MAM= 0.00270868710665567"
[1] "iteration 2 ; RMSM= 0.000170443626854365 ; MAM= 0.00107513354194388"
[1] "iteration 3 ; RMSM= 0.0000570936787209307 ; MAM= 0.000612171333340912"
[1] "stopped because of small RMSM= 0.000000345222848343858"
day nr.: 223
[1] "iteration 1 ; RMSM= 0.000883217338249451 ; MAM= 0.00659532301459982"
[1] "iteration 2 ; RMSM= 0.000438456764729958 ; MAM= 0.00528793223868559"
[1] "iteration 3 ; RMSM= 0.000109031374457191 ; MAM= 0.0011506291895391"
[1] "stopped because of small RMSM= 0.00000162797714619533"
day nr.: 224
[1] "iteration 1 ; RMSM= 0.0211853472183861 ; MAM= 0.315413977768936"
[1] "iteration 2 ; RMSM= 0.0163101049420066 ; MAM= 0.237773019568993"
[1] "iteration 3 ; RMSM= 0.0160265473375665 ; MAM= 0.235013972774068"
[1] "iteration 4 ; RMSM= 0.0158324435029292 ; MAM= 0.232242466591011"
[1] "iteration 5 ; RMSM= 0.0156617535935019 ; MAM= 0.230045174414424"
[1] "iteration 6 ; RMSM= 0.0153398852294146 ; MAM= 0.224964227286387"
[1] "iteration 7 ; RMSM= 0.0151196338346262 ; MAM= 0.221860469224141"
[1] "iteration 8 ; RMSM= 0.0149414747363414 ; MAM= 0.219446924240317"
[1] "iteration 9 ; RMSM= 0.0146702379790886 ; MAM= 0.215409778391453"
[1] "iteration 10 ; RMSM= 0.0143887757116614 ; MAM= 0.205750784557513"
[1] "iteration 11 ; RMSM= 0.0140666964475478 ; MAM= 0.197408385750051"
[1] "iteration 12 ; RMSM= 0.0131092368503333 ; MAM= 0.186795654038416"
[1] "iteration 13 ; RMSM= 0.012202028673309 ; MAM= 0.174343964434186"
[1] "iteration 14 ; RMSM= 0.0114302777659192 ; MAM= 0.164655033117302"
[1] "iteration 15 ; RMSM= 0.0107991283520728 ; MAM= 0.156288931753166"
[1] "iteration 16 ; RMSM= 0.0100539979503235 ; MAM= 0.145358250438999"
[1] "iteration 17 ; RMSM= 0.00805279785577359 ; MAM= 0.102655226513231"
[1] "iteration 18 ; RMSM= 0.00530528862435067 ; MAM= 0.069288354670133"
[1] "iteration 19 ; RMSM= 0.00420829852780637 ; MAM= 0.0547896574941817"
[1] "iteration 20 ; RMSM= 0.00335171530227514 ; MAM= 0.0429516729831673"
[1] "iteration 21 ; RMSM= 0.00276348700584979 ; MAM= 0.0162910707267805"
[1] "iteration 22 ; RMSM= 0.000359368036416108 ; MAM= 0.00204822888850592"
[1] "stopped because of small RMSM= 0.000000726309605023152"
day nr.: 225
[1] "iteration 1 ; RMSM= 0.00824183329897833 ; MAM= 0.114317775659092"
[1] "iteration 2 ; RMSM= 0.00303719629134137 ; MAM= 0.0402169016605431"
[1] "iteration 3 ; RMSM= 0.000610759341771511 ; MAM= 0.00735184290965763"
[1] "iteration 4 ; RMSM= 0.0000287931164732147 ; MAM= 0.000309377718523133"
[1] "stopped because of small RMSM= 0.0000000466876462668653"
day nr.: 226
[1] "iteration 1 ; RMSM= 0.000722734747336444 ; MAM= 0.0063187220852701"
[1] "iteration 2 ; RMSM= 0.000161219794568183 ; MAM= 0.00122997405216459"
[1] "stopped because of small RMSM= 0.00000950924468664727"
day nr.: 227
[1] "iteration 1 ; RMSM= 0.0248933458601506 ; MAM= 0.369534096564938"
[1] "iteration 2 ; RMSM= 0.0230618857821424 ; MAM= 0.342410768386324"
[1] "iteration 3 ; RMSM= 0.0217324728362409 ; MAM= 0.322726034073249"
[1] "iteration 4 ; RMSM= 0.0168179371954715 ; MAM= 0.201648876539782"
[1] "iteration 5 ; RMSM= 0.0116235454170397 ; MAM= 0.150042979883718"
[1] "iteration 6 ; RMSM= 0.00858465963618943 ; MAM= 0.106296828494875"
[1] "iteration 7 ; RMSM= 0.00653204668138017 ; MAM= 0.0821407812764807"
[1] "iteration 8 ; RMSM= 0.00644945208655771 ; MAM= 0.0320989333080772"
[1] "iteration 9 ; RMSM= 0.000844140360957931 ; MAM= 0.00394428642872027"
[1] "stopped because of small RMSM= 0.00000468459202791357"
day nr.: 228
[1] "iteration 1 ; RMSM= 0.0108573629485653 ; MAM= 0.151284914186911"
[1] "iteration 2 ; RMSM= 0.00286350102573675 ; MAM= 0.0389450459753754"
[1] "iteration 3 ; RMSM= 0.000286780930407055 ; MAM= 0.003618880421036"
[1] "stopped because of small RMSM= 0.00000198983508985184"
day nr.: 229
[1] "iteration 1 ; RMSM= 0.0010427675873351 ; MAM= 0.0105734428652557"
[1] "iteration 2 ; RMSM= 0.000168727096866819 ; MAM= 0.00164604229713414"
[1] "stopped because of small RMSM= 0.00000517688000879512"
day nr.: 230
[1] "iteration 1 ; RMSM= 0.000503071853212576 ; MAM= 0.00351930818418639"
[1] "iteration 2 ; RMSM= 0.0000361727374666508 ; MAM= 0.000324946317272258"
[1] "stopped because of small RMSM= 0.000000356864130555552"
day nr.: 231
[1] "iteration 1 ; RMSM= 0.00350289596050633 ; MAM= 0.0516940842777332"
[1] "iteration 2 ; RMSM= 0.00248935402835091 ; MAM= 0.0367138866894062"
[1] "iteration 3 ; RMSM= 0.00161505933844785 ; MAM= 0.0229349246511521"
[1] "iteration 4 ; RMSM= 0.000174940335656763 ; MAM= 0.00210015654485102"
[1] "stopped because of small RMSM= 0.000000934049883245952"
day nr.: 232
[1] "iteration 1 ; RMSM= 0.00163734982941902 ; MAM= 0.0161478064278506"
[1] "iteration 2 ; RMSM= 0.000319857237419163 ; MAM= 0.00218419619648935"
[1] "stopped because of small RMSM= 0.0000052292510840779"
day nr.: 233
[1] "iteration 1 ; RMSM= 0.0141365051917636 ; MAM= 0.183161639234772"
[1] "iteration 2 ; RMSM= 0.0103986206738952 ; MAM= 0.137929276297581"
[1] "iteration 3 ; RMSM= 0.00747251343536088 ; MAM= 0.0976583801413627"
[1] "iteration 4 ; RMSM= 0.00542522293957908 ; MAM= 0.0668211200121588"
[1] "iteration 5 ; RMSM= 0.00435801785344098 ; MAM= 0.0495804588960573"
[1] "iteration 6 ; RMSM= 0.00393613937223137 ; MAM= 0.0229388293172826"
[1] "iteration 7 ; RMSM= 0.00028342699412069 ; MAM= 0.00162076432961633"
[1] "stopped because of small RMSM= 0.000000139480478793741"
day nr.: 234
[1] "iteration 1 ; RMSM= 0.0113481646854021 ; MAM= 0.159638217537048"
[1] "iteration 2 ; RMSM= 0.00329846714925989 ; MAM= 0.0464251761694857"
[1] "iteration 3 ; RMSM= 0.000275071278857954 ; MAM= 0.00352166784932413"
[1] "stopped because of small RMSM= 0.00000142904234058318"
day nr.: 235
[1] "iteration 1 ; RMSM= 0.0100140552226113 ; MAM= 0.148716979378743"
[1] "iteration 2 ; RMSM= 0.00716909215407166 ; MAM= 0.106274354732171"
[1] "iteration 3 ; RMSM= 0.00389663166238944 ; MAM= 0.0536285929754063"
[1] "iteration 4 ; RMSM= 0.000443099130713972 ; MAM= 0.00460035162279823"
[1] "stopped because of small RMSM= 0.00000554016598584693"
day nr.: 236
[1] "iteration 1 ; RMSM= 0.00256469082630785 ; MAM= 0.029018419829399"
[1] "iteration 2 ; RMSM= 0.00238845938427627 ; MAM= 0.0164017447827254"
[1] "iteration 3 ; RMSM= 0.000164249336038126 ; MAM= 0.00113905167830268"
[1] "stopped because of small RMSM= 0.000000223843755717051"
day nr.: 237
[1] "iteration 1 ; RMSM= 0.0269795916227063 ; MAM= 0.395107415781203"
[1] "iteration 2 ; RMSM= 0.0202862843523862 ; MAM= 0.290475556405883"
[1] "iteration 3 ; RMSM= 0.0147373025855199 ; MAM= 0.207212405037495"
[1] "iteration 4 ; RMSM= 0.0143910222010679 ; MAM= 0.0969207872897035"
[1] "iteration 5 ; RMSM= 0.0023186573521802 ; MAM= 0.0146774411929326"
[1] "iteration 6 ; RMSM= 0.0000309285018178087 ; MAM= 0.000203273436248708"
[1] "stopped because of small RMSM= 0.00000000034430226208913"
day nr.: 238
[1] "iteration 1 ; RMSM= 0.0163509016930457 ; MAM= 0.231818065701678"
[1] "iteration 2 ; RMSM= 0.00229150057836991 ; MAM= 0.0322903151744721"
[1] "iteration 3 ; RMSM= 0.0000317018634662506 ; MAM= 0.000406694239251049"
[1] "stopped because of small RMSM= 0.00000000348384290789629"
day nr.: 239
[1] "iteration 1 ; RMSM= 0.0202907313951917 ; MAM= 0.29198774975955"
[1] "iteration 2 ; RMSM= 0.00250726011953735 ; MAM= 0.0347196635633625"
[1] "iteration 3 ; RMSM= 0.0000319153254979797 ; MAM= 0.000300955892716448"
[1] "stopped because of small RMSM= 0.00000000125022935838808"
day nr.: 240
[1] "iteration 1 ; RMSM= 0.00892186903131347 ; MAM= 0.129950257427915"
[1] "iteration 2 ; RMSM= 0.000962918137373543 ; MAM= 0.0139877996877759"
[1] "stopped because of small RMSM= 0.00000719262025804628"
day nr.: 241
[1] "iteration 1 ; RMSM= 0.0179369177511387 ; MAM= 0.266569892512735"
[1] "iteration 2 ; RMSM= 0.0109675117919296 ; MAM= 0.158519492665943"
[1] "iteration 3 ; RMSM= 0.0010638243935379 ; MAM= 0.0141707635629602"
[1] "stopped because of small RMSM= 0.00000462059765395344"
day nr.: 242
[1] "iteration 1 ; RMSM= 0.00574189039406472 ; MAM= 0.0832207094277401"
[1] "iteration 2 ; RMSM= 0.000363562531819227 ; MAM= 0.00525274904519206"
[1] "stopped because of small RMSM= 0.00000100492717236918"
day nr.: 243
[1] "iteration 1 ; RMSM= 0.00136458882369678 ; MAM= 0.020185748964194"
[1] "iteration 2 ; RMSM= 0.0000705300524724156 ; MAM= 0.0010390143452317"
[1] "stopped because of small RMSM= 0.000000091400288269981"
day nr.: 244
[1] "iteration 1 ; RMSM= 0.000190599807119507 ; MAM= 0.00155837668857547"
[1] "stopped because of small RMSM= 0.00000349122182247476"
day nr.: 245
[1] "iteration 1 ; RMSM= 0.0157761328773462 ; MAM= 0.214383069513689"
[1] "iteration 2 ; RMSM= 0.0140073612231978 ; MAM= 0.161352091108717"
[1] "iteration 3 ; RMSM= 0.0112638772067801 ; MAM= 0.113146174264533"
[1] "iteration 4 ; RMSM= 0.0087985496699275 ; MAM= 0.0799652678508269"
[1] "iteration 5 ; RMSM= 0.00700702491692191 ; MAM= 0.0668684224420286"
[1] "iteration 6 ; RMSM= 0.00520443656875937 ; MAM= 0.0490489446693991"
[1] "iteration 7 ; RMSM= 0.00517030319477796 ; MAM= 0.0437029246289227"
[1] "iteration 8 ; RMSM= 0.00032682838868538 ; MAM= 0.00269846148640918"
[1] "stopped because of small RMSM= 0.0000000672205202943773"
day nr.: 246
[1] "iteration 1 ; RMSM= 0.0122139996692372 ; MAM= 0.17317283270865"
[1] "iteration 2 ; RMSM= 0.000422931794685529 ; MAM= 0.00589014332367144"
[1] "stopped because of small RMSM= 0.000000522691842249871"
day nr.: 247
[1] "iteration 1 ; RMSM= 0.00197669959561837 ; MAM= 0.0247173984785953"
[1] "iteration 2 ; RMSM= 0.0000460978633068474 ; MAM= 0.000569483391327785"
[1] "stopped because of small RMSM= 0.00000000364509581523404"
day nr.: 248
[1] "iteration 1 ; RMSM= 0.00645388626101892 ; MAM= 0.0949809965994474"
[1] "iteration 2 ; RMSM= 0.00051958297779486 ; MAM= 0.00758041522044281"
[1] "stopped because of small RMSM= 0.00000207866005535543"
day nr.: 249
[1] "iteration 1 ; RMSM= 0.000730118959496214 ; MAM= 0.00815045702491883"
[1] "iteration 2 ; RMSM= 0.0000238074138152149 ; MAM= 0.000296284708257523"
[1] "stopped because of small RMSM= 0.0000000105027878081638"
day nr.: 250
[1] "iteration 1 ; RMSM= 0.000088573910714474 ; MAM= 0.000599339300987364"
[1] "stopped because of small RMSM= 0.000000459877136195682"
day nr.: 251
[1] "iteration 1 ; RMSM= 0.00059239576926552 ; MAM= 0.00862080642585969"
[1] "stopped because of small RMSM= 0.00000292936600341997"
day nr.: 252
[1] "iteration 1 ; RMSM= 0.00100122746419689 ; MAM= 0.0148210247909653"
[1] "iteration 2 ; RMSM= 0.0000244218469850215 ; MAM= 0.000355347792999028"
[1] "stopped because of small RMSM= 0.00000000735505160314502"
day nr.: 253
[1] "iteration 1 ; RMSM= 0.0000903870758973563 ; MAM= 0.000593907612755755"
[1] "stopped because of small RMSM= 0.00000202087623118348"
day nr.: 254
[1] "iteration 1 ; RMSM= 0.0000698997583693528 ; MAM= 0.000512097938198133"
[1] "stopped because of small RMSM= 0.000000207838623396488"
day nr.: 255
[1] "iteration 1 ; RMSM= 0.000128726405821452 ; MAM= 0.00063253774173995"
[1] "stopped because of small RMSM= 0.00000374934180692087"
day nr.: 256
[1] "iteration 1 ; RMSM= 0.000233511679267041 ; MAM= 0.0017858077360542"
[1] "stopped because of small RMSM= 0.00000474220158662343"
day nr.: 257
[1] "iteration 1 ; RMSM= 0.000273619021595894 ; MAM= 0.00184757885110604"
[1] "stopped because of small RMSM= 0.00000496522028507265"
day nr.: 258
[1] "iteration 1 ; RMSM= 0.000240753036782958 ; MAM= 0.00133084376777357"
[1] "stopped because of small RMSM= 0.00000119611169987048"
day nr.: 259
[1] "iteration 1 ; RMSM= 0.000244239148198212 ; MAM= 0.0016298328445182"
[1] "stopped because of small RMSM= 0.00000277674618875539"
day nr.: 260
[1] "iteration 1 ; RMSM= 0.000312400337503448 ; MAM= 0.00192069593508454"
[1] "stopped because of small RMSM= 0.00000383183948586525"
day nr.: 261
[1] "iteration 1 ; RMSM= 0.000291550966372066 ; MAM= 0.00160106325204308"
[1] "stopped because of small RMSM= 0.00000291521579738126"
day nr.: 262
[1] "iteration 1 ; RMSM= 0.000289837915316682 ; MAM= 0.00164273644060296"
[1] "stopped because of small RMSM= 0.00000214679203766018"
day nr.: 263
[1] "iteration 1 ; RMSM= 0.000245081899444774 ; MAM= 0.0011774644835427"
[1] "stopped because of small RMSM= 0.000000715384725397229"
day nr.: 264
[1] "iteration 1 ; RMSM= 0.000227000571194654 ; MAM= 0.00117335754691029"
[1] "stopped because of small RMSM= 0.000000377728715291482"
day nr.: 265
[1] "iteration 1 ; RMSM= 0.0127365963019392 ; MAM= 0.189215488315796"
[1] "iteration 2 ; RMSM= 0.00870311420601204 ; MAM= 0.128987905143461"
[1] "iteration 3 ; RMSM= 0.00472500241870482 ; MAM= 0.0574490523634174"
[1] "iteration 4 ; RMSM= 0.000529196205515012 ; MAM= 0.00534256517516535"
[1] "iteration 5 ; RMSM= 0.0000128316073759612 ; MAM= 0.000142973229352694"
[1] "stopped because of small RMSM= 0.000000000751123845842333"
day nr.: 266
[1] "iteration 1 ; RMSM= 0.00236860151077932 ; MAM= 0.0216816338269844"
[1] "iteration 2 ; RMSM= 0.000203531645024639 ; MAM= 0.00131527414765192"
[1] "stopped because of small RMSM= 0.00000267964983587408"
day nr.: 267
[1] "iteration 1 ; RMSM= 0.00183305555024428 ; MAM= 0.0257089793055889"
[1] "iteration 2 ; RMSM= 0.000049483250538438 ; MAM= 0.00060919705724459"
[1] "stopped because of small RMSM= 0.0000000114938249554383"
day nr.: 268
[1] "iteration 1 ; RMSM= 0.0539606108305613 ; MAM= 0.655168112838614"
[1] "iteration 2 ; RMSM= 0.0430809657632906 ; MAM= 0.448682246023699"
[1] "iteration 3 ; RMSM= 0.0306596416620167 ; MAM= 0.286406636535079"
[1] "iteration 4 ; RMSM= 0.0217049011504147 ; MAM= 0.212443208867506"
[1] "iteration 5 ; RMSM= 0.0152629154427392 ; MAM= 0.136031882304847"
[1] "iteration 6 ; RMSM= 0.0110313833283424 ; MAM= 0.10081896070057"
[1] "iteration 7 ; RMSM= 0.00855791444512155 ; MAM= 0.0655266340160399"
[1] "iteration 8 ; RMSM= 0.000509832674110284 ; MAM= 0.00326999318807086"
[1] "stopped because of small RMSM= 0.0000000470535749456772"
day nr.: 269
[1] "iteration 1 ; RMSM= 0.0372682419235083 ; MAM= 0.535352273086731"
[1] "iteration 2 ; RMSM= 0.0051300352158345 ; MAM= 0.0734215173814008"
[1] "iteration 3 ; RMSM= 0.0000782808789128077 ; MAM= 0.00108463465311742"
[1] "stopped because of small RMSM= 0.0000000110644999187748"
day nr.: 270
[1] "iteration 1 ; RMSM= 0.0020931627458541 ; MAM= 0.0262835317565817"
[1] "iteration 2 ; RMSM= 0.0000683797082466586 ; MAM= 0.000581357396985116"
[1] "stopped because of small RMSM= 0.0000000264766900324813"
day nr.: 271
[1] "iteration 1 ; RMSM= 0.00802490263039734 ; MAM= 0.117904082329876"
[1] "iteration 2 ; RMSM= 0.0002551662482225 ; MAM= 0.00366151203677045"
[1] "stopped because of small RMSM= 0.000000175833294601088"
day nr.: 272
[1] "iteration 1 ; RMSM= 0.0110617179640411 ; MAM= 0.16328966545787"
[1] "iteration 2 ; RMSM= 0.00117718580091374 ; MAM= 0.0172085766547218"
[1] "stopped because of small RMSM= 0.00000783995381384104"
day nr.: 273
[1] "iteration 1 ; RMSM= 0.0132332501708194 ; MAM= 0.194037455179519"
[1] "iteration 2 ; RMSM= 0.000776284214831653 ; MAM= 0.0110949990624938"
[1] "stopped because of small RMSM= 0.00000162397817363302"
day nr.: 274
[1] "iteration 1 ; RMSM= 0.00667205320654051 ; MAM= 0.0987893434621181"
[1] "iteration 2 ; RMSM= 0.000459327770092094 ; MAM= 0.00665003271312081"
[1] "stopped because of small RMSM= 0.00000129087965019961"
day nr.: 275
[1] "iteration 1 ; RMSM= 0.0291011937033244 ; MAM= 0.432953076184261"
[1] "iteration 2 ; RMSM= 0.0201693975182233 ; MAM= 0.298201596363094"
[1] "iteration 3 ; RMSM= 0.00209250106261372 ; MAM= 0.0288050476663315"
[1] "iteration 4 ; RMSM= 0.000015899906247264 ; MAM= 0.00015655254155289"
[1] "stopped because of small RMSM= 0.000000000245505803384065"
day nr.: 276
[1] "iteration 1 ; RMSM= 0.00748975484887064 ; MAM= 0.10988113020257"
[1] "iteration 2 ; RMSM= 0.000339869984506973 ; MAM= 0.00492314406646771"
[1] "stopped because of small RMSM= 0.000000437945066718044"
day nr.: 277
[1] "iteration 1 ; RMSM= 0.000139620942960992 ; MAM= 0.00110185823462292"
[1] "stopped because of small RMSM= 0.000000807810628125344"
day nr.: 278
[1] "iteration 1 ; RMSM= 0.034544910465039 ; MAM= 0.510152309332995"
[1] "iteration 2 ; RMSM= 0.00484115326538444 ; MAM= 0.0674031661785223"
[1] "iteration 3 ; RMSM= 0.0000638353518896932 ; MAM= 0.000805615551972005"
[1] "stopped because of small RMSM= 0.00000000576933289171683"
day nr.: 279
[1] "iteration 1 ; RMSM= 0.00364070437688068 ; MAM= 0.0527369294379093"
[1] "iteration 2 ; RMSM= 0.0000533328286379203 ; MAM= 0.000766068927705088"
[1] "stopped because of small RMSM= 0.00000000703356556375847"
day nr.: 280
[1] "iteration 1 ; RMSM= 0.0635112650551568 ; MAM= 0.945879133471334"
[1] "iteration 2 ; RMSM= 0.0448351542600806 ; MAM= 0.667326451063605"
[1] "iteration 3 ; RMSM= 0.0198849817685081 ; MAM= 0.221524582729351"
[1] "iteration 4 ; RMSM= 0.00269800400699088 ; MAM= 0.027538907672195"
[1] "iteration 5 ; RMSM= 0.0000619606198869165 ; MAM= 0.00060500740297087"
[1] "stopped because of small RMSM= 0.000000000577535456828387"
day nr.: 281
[1] "iteration 1 ; RMSM= 0.0380731204790036 ; MAM= 0.552490599447382"
[1] "iteration 2 ; RMSM= 0.0065925320721936 ; MAM= 0.0963706722240774"
[1] "iteration 3 ; RMSM= 0.000147403019100168 ; MAM= 0.00211197116701903"
[1] "stopped because of small RMSM= 0.000000212841355881047"
day nr.: 282
[1] "iteration 1 ; RMSM= 0.0731932648681806 ; MAM= 1.08915395596631"
[1] "iteration 2 ; RMSM= 0.0606348897294198 ; MAM= 0.898202353918626"
[1] "iteration 3 ; RMSM= 0.00697954698610895 ; MAM= 0.0940043311376391"
[1] "iteration 4 ; RMSM= 0.000107218894488674 ; MAM= 0.000978556247325235"
[1] "stopped because of small RMSM= 0.00000000609740174852008"
day nr.: 283
[1] "iteration 1 ; RMSM= 0.0403146316354989 ; MAM= 0.590537164777351"
[1] "iteration 2 ; RMSM= 0.00795728123350355 ; MAM= 0.117235475698651"
[1] "iteration 3 ; RMSM= 0.000241560353430916 ; MAM= 0.00344411345466486"
[1] "stopped because of small RMSM= 0.000000126234322361534"
day nr.: 284
[1] "iteration 1 ; RMSM= 0.00235376820778041 ; MAM= 0.018940307670191"
[1] "iteration 2 ; RMSM= 0.000127527672820356 ; MAM= 0.00109273273302565"
[1] "stopped because of small RMSM= 0.000000432600763953225"
day nr.: 285
[1] "iteration 1 ; RMSM= 0.00106086783397575 ; MAM= 0.00715621660356948"
[1] "iteration 2 ; RMSM= 0.000030706680027812 ; MAM= 0.000183299924256641"
[1] "stopped because of small RMSM= 0.0000000626345928563934"
day nr.: 286
[1] "iteration 1 ; RMSM= 0.000652557241025871 ; MAM= 0.00404296133524057"
[1] "stopped because of small RMSM= 0.00000898228902914463"
day nr.: 287
[1] "iteration 1 ; RMSM= 0.00074257309427954 ; MAM= 0.00667056266796795"
[1] "stopped because of small RMSM= 0.00000158840530792291"
day nr.: 288
[1] "iteration 1 ; RMSM= 0.000467498544002579 ; MAM= 0.00463094475178139"
[1] "stopped because of small RMSM= 0.000000860033310985211"
day nr.: 289
[1] "iteration 1 ; RMSM= 0.000302564541810282 ; MAM= 0.00317903826598284"
[1] "stopped because of small RMSM= 0.000000432075750153085"
day nr.: 290
[1] "iteration 1 ; RMSM= 0.00051105280753552 ; MAM= 0.0056851988192835"
[1] "stopped because of small RMSM= 0.000000551244121379233"
day nr.: 291
[1] "iteration 1 ; RMSM= 0.0004427982455154 ; MAM= 0.00533418891071635"
[1] "stopped because of small RMSM= 0.0000003052293764436"
day nr.: 292
[1] "iteration 1 ; RMSM= 0.000180048995283497 ; MAM= 0.00186092687367354"
[1] "stopped because of small RMSM= 0.0000000776558260399807"
day nr.: 293
[1] "iteration 1 ; RMSM= 0.0327881365665333 ; MAM= 0.488120991877256"
[1] "iteration 2 ; RMSM= 0.0304426245451847 ; MAM= 0.445555419221569"
[1] "iteration 3 ; RMSM= 0.00517228463062785 ; MAM= 0.0734169894360507"
[1] "iteration 4 ; RMSM= 0.000217603032751908 ; MAM= 0.00261033081689288"
[1] "stopped because of small RMSM= 0.0000000558251955561502"
day nr.: 294
[1] "iteration 1 ; RMSM= 0.0142459764753175 ; MAM= 0.210411369236774"
[1] "iteration 2 ; RMSM= 0.00139440757977897 ; MAM= 0.0201023541586796"
[1] "iteration 3 ; RMSM= 0.0000118112939037994 ; MAM= 0.00016215249485437"
[1] "stopped because of small RMSM= 0.000000000421042127521272"
day nr.: 295
[1] "iteration 1 ; RMSM= 0.000406744558914703 ; MAM= 0.00321624463684846"
[1] "stopped because of small RMSM= 0.00000942230714910851"
day nr.: 296
[1] "iteration 1 ; RMSM= 0.000190345171988043 ; MAM= 0.00121689814799466"
[1] "stopped because of small RMSM= 0.000000621446229237433"
day nr.: 297
[1] "iteration 1 ; RMSM= 0.00717307564189011 ; MAM= 0.105864785886696"
[1] "iteration 2 ; RMSM= 0.000302705914554912 ; MAM= 0.00437824608052245"
[1] "stopped because of small RMSM= 0.000000348744066267921"
day nr.: 298
[1] "iteration 1 ; RMSM= 0.00349306972101482 ; MAM= 0.0515878039133579"
[1] "iteration 2 ; RMSM= 0.000108241537168127 ; MAM= 0.00153665305815337"
[1] "stopped because of small RMSM= 0.0000000647459895619324"
day nr.: 299
[1] "iteration 1 ; RMSM= 0.0151341098234442 ; MAM= 0.222712899747388"
[1] "iteration 2 ; RMSM= 0.00102523970280595 ; MAM= 0.0147744908524989"
[1] "stopped because of small RMSM= 0.00000359822881564933"
day nr.: 300
[1] "iteration 1 ; RMSM= 0.0533248618698116 ; MAM= 0.794137320682459"
[1] "iteration 2 ; RMSM= 0.036592524315627 ; MAM= 0.544157329155353"
[1] "iteration 3 ; RMSM= 0.0174747930926758 ; MAM= 0.218924048784815"
[1] "iteration 4 ; RMSM= 0.0027944924141034 ; MAM= 0.0291866363798503"
[1] "iteration 5 ; RMSM= 0.0000783676218548998 ; MAM= 0.000915073177506054"
[1] "stopped because of small RMSM= 0.000000000763745074253186"
day nr.: 301
[1] "iteration 1 ; RMSM= 0.00433898740851549 ; MAM= 0.0497503879184348"
[1] "iteration 2 ; RMSM= 0.000148416512096851 ; MAM= 0.001364087178235"
[1] "stopped because of small RMSM= 0.0000000180551335025896"
day nr.: 302
[1] "iteration 1 ; RMSM= 0.0146002833947132 ; MAM= 0.212644958070913"
[1] "iteration 2 ; RMSM= 0.00123185025110261 ; MAM= 0.0179904024597792"
[1] "stopped because of small RMSM= 0.00000959348178619891"
day nr.: 303
[1] "iteration 1 ; RMSM= 0.00188169187522953 ; MAM= 0.0192917907581845"
[1] "stopped because of small RMSM= 0.0000033966760467917"
day nr.: 304
[1] "iteration 1 ; RMSM= 0.000441923404328802 ; MAM= 0.0029498306019981"
[1] "stopped because of small RMSM= 0.00000217031637437115"
day nr.: 305
[1] "iteration 1 ; RMSM= 0.00163501289781572 ; MAM= 0.0176588333978526"
[1] "stopped because of small RMSM= 0.000000471575127032077"
day nr.: 306
[1] "iteration 1 ; RMSM= 0.0093291171653216 ; MAM= 0.137514262580624"
[1] "iteration 2 ; RMSM= 0.000715473511964713 ; MAM= 0.0103815348043897"
[1] "stopped because of small RMSM= 0.00000283140596092175"
day nr.: 307
[1] "iteration 1 ; RMSM= 0.000712360104226268 ; MAM= 0.00515871414208396"
[1] "iteration 2 ; RMSM= 0.0000184444997394901 ; MAM= 0.000167917049521821"
[1] "stopped because of small RMSM= 0.00000000126262507058922"
day nr.: 308
[1] "iteration 1 ; RMSM= 0.000500185602665332 ; MAM= 0.00300939222134848"
[1] "stopped because of small RMSM= 0.00000408912169714781"
day nr.: 309
[1] "iteration 1 ; RMSM= 0.000636335045718744 ; MAM= 0.00623264929225076"
[1] "stopped because of small RMSM= 0.00000208740639348807"
day nr.: 310
[1] "iteration 1 ; RMSM= 0.000895745853517001 ; MAM= 0.00985909541554342"
[1] "stopped because of small RMSM= 0.000000672647335214803"
day nr.: 311
[1] "iteration 1 ; RMSM= 0.000174629660872578 ; MAM= 0.000796815245510646"
[1] "stopped because of small RMSM= 0.000000929879007168971"
day nr.: 312
[1] "iteration 1 ; RMSM= 0.000721949487345237 ; MAM= 0.00800000983577107"
[1] "stopped because of small RMSM= 0.000000745198590211259"
day nr.: 313
[1] "iteration 1 ; RMSM= 0.000845700062593151 ; MAM= 0.0095153522615993"
[1] "stopped because of small RMSM= 0.000000402658017489086"
day nr.: 314
[1] "iteration 1 ; RMSM= 0.000131860353053701 ; MAM= 0.00151910756096441"
[1] "stopped because of small RMSM= 0.000000184647169222877"
day nr.: 315
[1] "iteration 1 ; RMSM= 0.0213772568822742 ; MAM= 0.314641466012132"
[1] "iteration 2 ; RMSM= 0.00278305842877407 ; MAM= 0.039388421361336"
[1] "iteration 3 ; RMSM= 0.0000259962076913618 ; MAM= 0.000353255728591184"
[1] "stopped because of small RMSM= 0.00000000139726394983913"
day nr.: 316
[1] "iteration 1 ; RMSM= 0.000740537618357882 ; MAM= 0.009408722447768"
[1] "stopped because of small RMSM= 0.00000314251997833316"
day nr.: 317
[1] "iteration 1 ; RMSM= 0.00134547384555461 ; MAM= 0.0198399377816803"
[1] "iteration 2 ; RMSM= 0.0000234510809654065 ; MAM= 0.000306551257427519"
[1] "stopped because of small RMSM= 0.00000000319224802276502"
day nr.: 318
[1] "iteration 1 ; RMSM= 0.039530010324458 ; MAM= 0.588840083373537"
[1] "iteration 2 ; RMSM= 0.0393842704060014 ; MAM= 0.574246502763867"
[1] "iteration 3 ; RMSM= 0.00541766495533203 ; MAM= 0.0734526498075319"
[1] "iteration 4 ; RMSM= 0.000205647295542174 ; MAM= 0.00271352766688659"
[1] "stopped because of small RMSM= 0.0000000257500350672023"
day nr.: 319
[1] "iteration 1 ; RMSM= 0.0166877186475735 ; MAM= 0.24505049759061"
[1] "iteration 2 ; RMSM= 0.00141658593354786 ; MAM= 0.020549007201964"
[1] "stopped because of small RMSM= 0.00000634448615972944"
day nr.: 320
[1] "iteration 1 ; RMSM= 0.000801822518888078 ; MAM= 0.0101263111240657"
[1] "iteration 2 ; RMSM= 0.0000262645858504528 ; MAM= 0.000300075536428454"
[1] "stopped because of small RMSM= 0.00000000315295591545336"
day nr.: 321
[1] "iteration 1 ; RMSM= 0.000673824698352558 ; MAM= 0.00872313851540028"
[1] "iteration 2 ; RMSM= 0.0000136993268124817 ; MAM= 0.00014262900118858"
[1] "stopped because of small RMSM= 0.00000000132867725865332"
day nr.: 322
[1] "iteration 1 ; RMSM= 0.00102606272011859 ; MAM= 0.0148650077128816"
[1] "stopped because of small RMSM= 0.00000657133552235984"
day nr.: 323
[1] "iteration 1 ; RMSM= 0.00132602372909507 ; MAM= 0.0191953132283572"
[1] "iteration 2 ; RMSM= 0.000020848980658582 ; MAM= 0.00028493262685968"
[1] "stopped because of small RMSM= 0.0000000026539080074105"
day nr.: 324
[1] "iteration 1 ; RMSM= 0.000111174521275274 ; MAM= 0.000757698019714814"
[1] "stopped because of small RMSM= 0.000000244538924384389"
day nr.: 325
[1] "iteration 1 ; RMSM= 0.00633259155746277 ; MAM= 0.0933144183884949"
[1] "iteration 2 ; RMSM= 0.000252065559974362 ; MAM= 0.00362653643675148"
[1] "stopped because of small RMSM= 0.000000254953549772169"
day nr.: 326
[1] "iteration 1 ; RMSM= 0.0033479028553855 ; MAM= 0.0494965090097556"
[1] "iteration 2 ; RMSM= 0.000110898543179292 ; MAM= 0.00158139134638647"
[1] "stopped because of small RMSM= 0.0000000655850463274961"
day nr.: 327
[1] "iteration 1 ; RMSM= 0.00021231872429172 ; MAM= 0.00204131069301089"
[1] "stopped because of small RMSM= 0.000000286523188098023"
day nr.: 328
[1] "iteration 1 ; RMSM= 0.0000745774245642777 ; MAM= 0.000453668025632376"
[1] "stopped because of small RMSM= 0.0000000872273879670082"
day nr.: 329
[1] "iteration 1 ; RMSM= 0.0129140438696921 ; MAM= 0.190282311245288"
[1] "iteration 2 ; RMSM= 0.000976555247269296 ; MAM= 0.0142110872527554"
[1] "stopped because of small RMSM= 0.00000419278265082677"
day nr.: 330
[1] "iteration 1 ; RMSM= 0.000152796631971142 ; MAM= 0.00100239336459561"
[1] "stopped because of small RMSM= 0.00000113583611432917"
day nr.: 331
[1] "iteration 1 ; RMSM= 0.00193569196360127 ; MAM= 0.0281376541838539"
[1] "iteration 2 ; RMSM= 0.0000453061412179905 ; MAM= 0.000637403376438362"
[1] "stopped because of small RMSM= 0.0000000122606814708547"
day nr.: 332
[1] "iteration 1 ; RMSM= 0.000113364178193839 ; MAM= 0.000616669710410307"
[1] "stopped because of small RMSM= 0.000000930137948727568"
day nr.: 333
[1] "iteration 1 ; RMSM= 0.0339803179647569 ; MAM= 0.506190627396055"
[1] "iteration 2 ; RMSM= 0.0235986711026466 ; MAM= 0.351441785459134"
[1] "iteration 3 ; RMSM= 0.0105895044166201 ; MAM= 0.148761633647943"
[1] "iteration 4 ; RMSM= 0.0011988454096914 ; MAM= 0.0132668177594496"
[1] "stopped because of small RMSM= 0.00000670546341843385"
day nr.: 334
[1] "iteration 1 ; RMSM= 0.00153736187697858 ; MAM= 0.0200664864817137"
[1] "iteration 2 ; RMSM= 0.0000105203712348608 ; MAM= 0.000107391582343352"
[1] "stopped because of small RMSM= 0.0000000000745337069531795"
day nr.: 335
[1] "iteration 1 ; RMSM= 0.00745312505850025 ; MAM= 0.110582283499832"
[1] "iteration 2 ; RMSM= 0.000613069939751973 ; MAM= 0.00899537797423415"
[1] "stopped because of small RMSM= 0.00000217992597268832"
day nr.: 336
[1] "iteration 1 ; RMSM= 0.0202210326171381 ; MAM= 0.297086908472702"
[1] "iteration 2 ; RMSM= 0.00138240865016867 ; MAM= 0.0198684266592369"
[1] "stopped because of small RMSM= 0.00000506557315409523"
day nr.: 337
[1] "iteration 1 ; RMSM= 0.0061975199350258 ; MAM= 0.091826593313461"
[1] "iteration 2 ; RMSM= 0.000315038832645917 ; MAM= 0.00456081140042863"
[1] "stopped because of small RMSM= 0.000000403537775757022"
day nr.: 338
[1] "iteration 1 ; RMSM= 0.000995566407938156 ; MAM= 0.0137114778472644"
[1] "iteration 2 ; RMSM= 0.0000216793550776352 ; MAM= 0.000273223847153395"
[1] "stopped because of small RMSM= 0.00000000415780010710169"
day nr.: 339
[1] "iteration 1 ; RMSM= 0.0011339255098161 ; MAM= 0.0165194218436889"
[1] "stopped because of small RMSM= 0.00000663223155512427"
day nr.: 340
[1] "iteration 1 ; RMSM= 0.0365055470143573 ; MAM= 0.54380486494333"
[1] "iteration 2 ; RMSM= 0.0280023427616389 ; MAM= 0.408034302621838"
[1] "iteration 3 ; RMSM= 0.00387438748405728 ; MAM= 0.0532019180316183"
[1] "iteration 4 ; RMSM= 0.0000545658463403655 ; MAM= 0.000582643370924507"
[1] "stopped because of small RMSM= 0.00000000379522378253428"
day nr.: 341
[1] "iteration 1 ; RMSM= 0.0189368976958816 ; MAM= 0.279599813741316"
[1] "iteration 2 ; RMSM= 0.00236120347404085 ; MAM= 0.034572675427374"
[1] "iteration 3 ; RMSM= 0.0000237288200319234 ; MAM= 0.000324877977258931"
[1] "stopped because of small RMSM= 0.00000000125042658590755"
day nr.: 342
[1] "iteration 1 ; RMSM= 0.000621287191762222 ; MAM= 0.00598153348323857"
[1] "iteration 2 ; RMSM= 0.0000128257400079345 ; MAM= 0.0000889881471175871"
[1] "stopped because of small RMSM= 0.0000000183192197169606"
day nr.: 343
[1] "iteration 1 ; RMSM= 0.000355411574547495 ; MAM= 0.00251461397154737"
[1] "stopped because of small RMSM= 0.00000466631726340654"
day nr.: 344
[1] "iteration 1 ; RMSM= 0.00542471381131799 ; MAM= 0.0798264036176448"
[1] "iteration 2 ; RMSM= 0.000242849026903823 ; MAM= 0.00351217892154351"
[1] "stopped because of small RMSM= 0.000000223760454146897"
day nr.: 345
[1] "iteration 1 ; RMSM= 0.000294125459500818 ; MAM= 0.00411670706841608"
[1] "stopped because of small RMSM= 0.000000397135501756948"
day nr.: 346
[1] "iteration 1 ; RMSM= 0.000121948559872251 ; MAM= 0.00131113493252996"
[1] "stopped because of small RMSM= 0.00000012516705106435"
day nr.: 347
[1] "iteration 1 ; RMSM= 0.0330408295959283 ; MAM= 0.49226148323243"
[1] "iteration 2 ; RMSM= 0.0270782170473625 ; MAM= 0.401034112736857"
[1] "iteration 3 ; RMSM= 0.00318464699358725 ; MAM= 0.0438324434347224"
[1] "iteration 4 ; RMSM= 0.0000805738321122913 ; MAM= 0.000711781893255825"
[1] "stopped because of small RMSM= 0.0000000026760707120122"
day nr.: 348
[1] "iteration 1 ; RMSM= 0.00739696282895862 ; MAM= 0.104404325430115"
[1] "iteration 2 ; RMSM= 0.000184223585311893 ; MAM= 0.00246873632767564"
[1] "stopped because of small RMSM= 0.0000000644114325305668"
day nr.: 349
[1] "iteration 1 ; RMSM= 0.00508806837689697 ; MAM= 0.0751233583591981"
[1] "iteration 2 ; RMSM= 0.000303886533443326 ; MAM= 0.00439580103672773"
[1] "stopped because of small RMSM= 0.000000467226114805962"
day nr.: 350
[1] "iteration 1 ; RMSM= 0.0116611655226077 ; MAM= 0.171134397648452"
[1] "iteration 2 ; RMSM= 0.000624237651318234 ; MAM= 0.00902367009091132"
[1] "stopped because of small RMSM= 0.00000114163550832599"
day nr.: 351
[1] "iteration 1 ; RMSM= 0.00588345431501028 ; MAM= 0.0867130334217397"
[1] "iteration 2 ; RMSM= 0.000291873540934448 ; MAM= 0.00419219464617334"
[1] "stopped because of small RMSM= 0.000000448342977367924"
day nr.: 352
[1] "iteration 1 ; RMSM= 0.000338947837698477 ; MAM= 0.00270238123782091"
[1] "stopped because of small RMSM= 0.00000438518411739834"
day nr.: 353
[1] "iteration 1 ; RMSM= 0.000288291709513186 ; MAM= 0.0022535060597359"
[1] "stopped because of small RMSM= 0.00000110330284570586"
day nr.: 354
[1] "iteration 1 ; RMSM= 0.0534144113467585 ; MAM= 0.789453017200251"
[1] "iteration 2 ; RMSM= 0.0497173729851875 ; MAM= 0.449807549811983"
[1] "iteration 3 ; RMSM= 0.0187021538238595 ; MAM= 0.121468792526579"
[1] "iteration 4 ; RMSM= 0.00223498197847115 ; MAM= 0.0161938375779642"
[1] "iteration 5 ; RMSM= 0.0000151460983867945 ; MAM= 0.000144904714197924"
[1] "stopped because of small RMSM= 0.0000000000375544372621136"
day nr.: 355
[1] "iteration 1 ; RMSM= 0.021098647578484 ; MAM= 0.298839668505958"
[1] "iteration 2 ; RMSM= 0.000936817626109154 ; MAM= 0.012878628071518"
[1] "stopped because of small RMSM= 0.00000123994718404481"
day nr.: 356
[1] "iteration 1 ; RMSM= 0.0652297687025515 ; MAM= 0.970587490506639"
[1] "iteration 2 ; RMSM= 0.036880179119359 ; MAM= 0.517862640367068"
[1] "iteration 3 ; RMSM= 0.00372023769677041 ; MAM= 0.0444812022620105"
[1] "iteration 4 ; RMSM= 0.0000193363296672207 ; MAM= 0.000210221127263654"
[1] "stopped because of small RMSM= 0.000000000282599573899535"
day nr.: 357
[1] "iteration 1 ; RMSM= 0.0265179937569624 ; MAM= 0.382652368751865"
[1] "iteration 2 ; RMSM= 0.00200182696327529 ; MAM= 0.0287925757909128"
[1] "iteration 3 ; RMSM= 0.0000214456794651731 ; MAM= 0.000262001702440839"
[1] "stopped because of small RMSM= 0.00000000033691922324411"
day nr.: 358
[1] "iteration 1 ; RMSM= 0.0128431104461694 ; MAM= 0.17850634143087"
[1] "iteration 2 ; RMSM= 0.00141002903471654 ; MAM= 0.0200113897253808"
[1] "iteration 3 ; RMSM= 0.0000126034933064742 ; MAM= 0.000159442222242828"
[1] "stopped because of small RMSM= 0.000000000412845141650055"
day nr.: 359
[1] "iteration 1 ; RMSM= 0.00161281069939139 ; MAM= 0.0144163791601684"
[1] "iteration 2 ; RMSM= 0.0000455534552416907 ; MAM= 0.000361663344430219"
[1] "stopped because of small RMSM= 0.0000000169753706186998"
day nr.: 360
[1] "iteration 1 ; RMSM= 0.0141588288675232 ; MAM= 0.209260122170753"
[1] "iteration 2 ; RMSM= 0.000627221582853378 ; MAM= 0.00916870417825677"
[1] "stopped because of small RMSM= 0.000000886871133317532"
day nr.: 361
[1] "iteration 1 ; RMSM= 0.0065349650409892 ; MAM= 0.0963955615883269"
[1] "iteration 2 ; RMSM= 0.000237255264387203 ; MAM= 0.00338300352118304"
[1] "stopped because of small RMSM= 0.000000179899658716836"
day nr.: 362
[1] "iteration 1 ; RMSM= 0.00128582283314374 ; MAM= 0.0153139076608193"
[1] "iteration 2 ; RMSM= 0.0000347973544165324 ; MAM= 0.000407920404563375"
[1] "stopped because of small RMSM= 0.00000000472975498952675"
day nr.: 363
[1] "iteration 1 ; RMSM= 0.00181098033304034 ; MAM= 0.0155599388730577"
[1] "iteration 2 ; RMSM= 0.0000177648834557851 ; MAM= 0.000201003873013784"
[1] "stopped because of small RMSM= 0.00000000186691984508656"
day nr.: 364
[1] "iteration 1 ; RMSM= 0.00171621135240449 ; MAM= 0.0185511288144851"
[1] "stopped because of small RMSM= 0.00000264073281719587"
day nr.: 365 
end.run.time = Sys.time()
cat(paste("simulation time is :",end.run.time - start.run.time))simulation time is : 7.53177636464437## to make dataframes and csv's from the saved data
balance.res = data.frame(balance.tmp)
time.series = seq(begin.time,end.time,by = delta.t) #add this to the balanc.res data.frame for use with IJPlot
balance.res = cbind(time.series,balance.res)
names(balance.res) = c("Daynr","Precip","Epot","Eact","Sto2flow","Flow2sto","Recharge","Ponding")
write.csv(x = balance.res, file = "Balance.csv")
state.res = data.frame(state.tmp)
names(state.res) = as.character(nodes)
write.csv(x = state.res, file = "States.csv")
transp.res = data.frame(transp.tmp)
names(transp.res) = as.character(nodes)
write.csv(x = transp.res, file = "Transp.csv")
moist.res = data.frame(moist.tmp)
names(moist.res) = as.character(nodes)
write.csv(x = moist.res, file = "Moisture.csv")
## creating the base output data.frames
balance.base = balance.res
state.base = state.res
transp.base = transp.res
moist.base = moist.res
Viewing the results
First a base plot function to have a look at the water balance terms in time Next we use MIPLOT package from Paul Torfs (former colleague) to look at the moisture, pressure and root uptake distribution in the model in time. These are thus time series
plot.trans.balance = function() # this is a very basic simple and generic R plot
{
time.series = c(1:length(balance.res[,1]))
flux.range = range(balance.res[,2:8])
plot(x = time.series , y = balance.res$Eact, ylim = flux.range, xlab = "time (d)",
ylab = "flux rate cm/d", main = "Water balance terms", type = "l",
lwd = 2, col = "red")
lines(x = time.series, y = balance.res$Sto2flow, lwd = 2,col = "lightblue")
lines(x = time.series, y = balance.res$Flow2sto, lwd = 2, col = "lightgreen")
lines(x = time.series, y = balance.res$Recharge, lwd = 2, col = "magenta")
lines(x = time.series, y = balance.res$Precip, lwd = 2, col = 'blue')
lines(x = time.series, y = balance.res$Ponding, lwd = 2, col = 'orange')
lines(x = time.series, y = balance.res$Epot, lwd = 2, col = 'tomato3')
grid()
}
plot.trans.balance() # this a very basic simple and generic R plot
#below the IJPLOT package is used to plot different type of data
balance.plot = new_miplot(name = "balance",fxright = 0.90, rightmargin = 0.0)
color.set = rainbow(n = length(balance.res[1,]))
add_line(balance.plot, data = balance.res,dataxname = names(balance.res[1]), datayname = names(balance.res[2]), color = color.set[2],width = 2)
add_line(balance.plot, data = balance.res,dataxname = names(balance.res[1]), datayname = names(balance.res[3]), color = color.set[3],width = 2)
add_line(balance.plot, data = balance.res,dataxname = names(balance.res[1]),datayname = names(balance.res[4]), color = color.set[4],width = 2)
add_line(balance.plot, data = balance.res,dataxname = names(balance.res[1]),datayname = names(balance.res[5]), color = color.set[5],width = 2)
add_line(balance.plot, data = balance.res,dataxname = names(balance.res[1]),datayname = names(balance.res[6]), color = color.set[6],width = 2)
add_line(balance.plot, data = balance.res,dataxname = names(balance.res[1]),datayname = names(balance.res[7]), color = color.set[7],width = 2)
add_line(balance.plot, data = balance.res,dataxname = names(balance.res[1]),datayname = names(balance.res[8]), color = color.set[8],width = 2)
add_xlabel(balance.plot, text = "time steps (d)", fontsize = 12)
add_ylabel(balance.plot, text = "flux rate (cm/d)", fontsize = 12, offset = 25)
add_title(balance.plot, text = "Water balance terms")
add_vline(balance.plot, x = seq(min(balance.res[,1]),max(balance.res[,1]),length.out = 5), dash=c(2,2), color = "grey")
add_hline(balance.plot, y = seq(min(balance.res[,2:8]),max(balance.res[,2:8]),length.out = 5),dash=c(2,2), color = "grey")
legend.plot = new_miplot(name = "legend",yfixed = TRUE, xfixed = TRUE,fxleft = 0.90,fybottom = 0.7,
bottommargin = 0, leftmargin = 0, horaxisnumdiv = 0, veraxisnumdiv = 0,
imposedmaxx= max(nchar(names(balance.res))))
for (i in 2:length(balance.res[1,]))
{
add_points(mip = legend.plot,x = 0, y = length(balance.res[1,])-i, marker = "square",size = 8, fillcolor = color.set[i])
add_text(mip = legend.plot, atx= 0, aty = length(balance.res[1,])-i,text = paste(names(balance.res[i])),fontsize = 12,
pos = "E",offset = 5)
}
###for testing
write.table(moist.res,file = "moisture.dat")
write.table(transp.res, file = "transp.dat")
###
draw(balance.plot,legend.plot)Registered S3 method overwritten by 'htmlwidgets':
method from
print.htmlwidget tools:rstudioflux.plot = new_ijplot(name = "flux", fxright = 0.90, rightmargin = 0.0)Error in new_ijplot(name = "flux", fxright = 0.9, rightmargin = 0) :
could not find function "new_ijplot"LS0tDQp0aXRsZTogIlRyYW5zaWVudCBTaW11bGF0aW9ucyBVbnNhdHVyYXRlZCBab25lIg0KYXV0aG9yOiAiRy4gQmllciINCmRhdGU6ICJBcHJpbCwgMjAyNSINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNCmBgYHtyIHNldHVwLCBpbmNsdWRlPUZBTFNFfQ0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KGVjaG8gPSBUUlVFKQ0KYGBgDQoNCg0KIyBJbnRyb2R1Y3Rpb24NCg0KTWFya2Rvd24gdG8gc2V0IHVwIGFuIHVuc2F0dXJhdGVkIHpvbmUgbW9kZWwgYW5kIHRvIGFuYWx5emUgdGhlIGltcGFjdCBvZiB0aGUgdmFuIEdlbnVjaHRlbiBwYXJhbWV0ZXJzIG9uIEV2YXBvdHJhbnNwaXJhdGlvbiwgIG1vaXN0dXJlIGRpc3RyaWJ1dGlvbiBhbmQgZ3JvdW5kd2F0ZXIgcmVjaGFyZ2UuDQoNCg0KIyMjIENsZWFuaW5nIHRoZSBtZW1vcnkgYW5kIGxvYWRpbmcgcmVxdWlyZWQgcGFja2FnZShzKSANCmBgYHtyfQ0Kcm0obGlzdCA9IGxzKCkpDQpsaWJyYXJ5KEZWRkUxRCkNCmxpYnJhcnkoTUlQTE9UKQ0KYGBgDQoNCg0KIyMjIExvYWRpbmcgdGhlIHNvaWwgZGF0YSAoU3RhcmluZyBTZXJpZXMpDQoNCmBgYHtyfQ0KbG9hZChmaWxlPSdTdGFyaW5nLlJkYXRhJykNCnByaW50KHBhc3RlKCJTb2lsIHR5cGU6ICIsbmFtZXMoc29pbC5zZXQpLCIgIixsYXBwbHkoc29pbC5zZXQsIGZ1bmN0aW9uKHgpIHgkbmFtZSkpKQ0KDQpwcmludChwYXN0ZSgiS3NhdCB2YWx1ZXM6ICIsbGFwcGx5KHNvaWwuc2V0LCBmdW5jdGlvbih5KSB5JG5hbWUpLCIgICIsbGFwcGx5KHNvaWwuc2V0LCBmdW5jdGlvbih4KSB4JGtzYXQpKSkNCmBgYA0KDQoNCg0KIyBNb2RlbCBzZXR1cA0KDQojIyBJbnRlcm5hbCBmbHV4IGluIHRoZSB1bnNhdHVyYXRlZCB6b25lICANCg0KVGhlIGVxdWF0aW9uIGRlc2NyaWJpbmcgdGhlIGludGVybmFsIGZsdXggaW4gdGhlIHVuc2F0dXJhdGVkIHpvbmUgc3Ryb25nbHkgcmVzZW1ibGVzIHRoZSBiYXNpYyBEYXJjeSBlcXVhdGlvbiB0byBzaW11bGF0ZSBmbG93IGluIHRoZSBzYXR1cmF0ZWQgem9uZS4gTW9zdCBpbXBvcnRhbnQgZGlmZmVyZW5jZXMgYXJlIHRoZSBoeWRyYXVsaWMgY29uZHVjdGl2aXR5IGFuZCB0aGUgc3RhdGUuDQpOb3cgdGhlIHN0YXRlIGlzICRccHNpJCBpbnN0ZWFkIG9mICRIJC4gQm90aCBhcmUgcmVsYXRlZDogJEggPSBccHNpICsgeiQgd2hlcmUgdGhlICR6JCBpcyB0aGUgZGF0dW0vcmVmZXJlbmNlIChpbiB0aGUgTmV0aGVybGFuZHM6ICJOQVAiKS4NCg0KJCQNCnE9LWsoXHBzaSlcbGVmdCBbXGZyYWN7XHBhcnRpYWwgXHBzaX17XHBhcnRpYWwgen0gKyAxIFxyaWdodF0NCiQkDQoNCiMjIEJ1aWxkaW5nIHRoZSBtb2RlbCBmb3IgdGhlIEhvbHRlcmJlZw0KDQoqKkJlIGF3YXJlKiogdGhlIGxlbmd0aCB1bml0IGlzIGluIGNtIQ0KDQpUaGUgbW9kZWwgaXMgYmFzZWQgb246ICANCg0KKiAgMUQgdmVydGljYWwgZmxvdyANCiogIFJpY2hhcmRzL0RhcmN5J3MgZXF1YXRpb24gIA0KKiAgVG9wIGJvdW5kYXJ5IGZsdXggYmFzZWQgcHJlY2lwaXRhdGlvbg0KKiAgQm90dG9tIGJvdW5kYXJ5IGZsdXggYmFzZWQgb24gRGlyaWNobGV0IG9yIENhdWNoeSB0eXBlDQoqICBTcGF0aWFsIGZsdXggZm9yIHRoZSByb290aW5nIGZ1bmN0aW9uIGZvciB0cmFuc3BpcmF0aW9uDQoqICBTb2lsIGV2YXBvcmF0aW9uIGlzIG5lZ2xlY3RlZA0KDQojIyBDaG9vc2UgYSBzb2lsIHR5cGUgZnJvbSB0aGUgU3RhcmluZyBzZXJpZXMNCg0KRnJvbSB0aGUgU3RhcmluZyBTZXJpZXMgeW91IGNhbiBjaG9vc2UgYSB0eXBpY2FsIHNvaWwgdHlwZSBmb3IgaWNlIHB1c2hlZCByaWRnZXMuDQpGb3IgZXhhbXBsZSB0aGUgc29pbCB0eXBlICJDb2Fyc2Ugc2FuZCIgQjUuIA0KDQoNCmBgYHtyfQ0KDQpzb2lsdHlwZSA9IHNvaWwuc2V0JEIxDQoNCnByaW50KHVubGlzdChzb2lsdHlwZSkpDQoNCiMjIGRvbWFpbiBvZiB0aGUgbW9kZWwNCnRvcCA9IDAuMCAgIyBhdCB0aGUgc3VyZmFjZQ0KYm90dG9tID0gLTIwMCAjIHNvIDIwMCBjbSBiZWxvdyBzdXJmYWNlISENCg0KbG9nX3RvcCA9IDAuMQ0KbG9nX2JvdHRvbSA9IGxvZzEwKC1ib3R0b20pDQpsb2dfeiA9IHNlcShmcm9tID0gbG9nX3RvcCwgdG8gPSBsb2dfYm90dG9tLGJ5ID0gMC4wMSkNCm5vZGVzID0gLTEwXmxvZ196DQpub2RlcyA9IGMoMC4wLG5vZGVzKQ0KDQoNCmRvbWFpbj1jKGJvdHRvbSx0b3ApDQoNCg0KIyMjZGVzY3JpcHRpb24gb2YgdGhlIGludGVybmFsIGZsdXgjIyMjIyMjIyMjIyMjIw0KdW5zYXRkYXJjeS5mbHV4ID0gZnVuY3Rpb24oeixwc2ksZHBzaWR6KQ0Kew0KICByZXR1cm4oLXNvaWx0eXBlJGsuZnVuKHBzaSkqKGRwc2lkeisxKSkNCn0NCg0KIyMjIyMjIyMjIyMjIyMjIyMjIyBtYWtlIHRoZSBtb2RlbCAjIyMjIyMjIyMjIyMjIyMjIw0KVW5zYXRTdGF0ID0gbmV3RkxPVzFEKGRvbWFpbix1bnNhdGRhcmN5LmZsdXgsIlVuc2F0dXJhdGVkIHpvbmUsIDFEIHN0YXRpb25hcnkiKQ0KDQojIyMjIyMjIGRpc2NyZXRpc2F0aW9uICMjIyMjIyMjIyMjIyMjIyMjIyMjIw0KIyBwb3NpdGlvbiB0aGUgbm9kZXMgMiBjbSBhcGFydA0KI25vZGVzID0gc2VxKGZyb209ZG9tYWluWzFdLHRvPWRvbWFpblsyXSxieT0xKQ0KI3NldC5kaXNjcmV0aXNhdGlvbihVbnNhdFN0YXQsbm9kZXMsIkZWIikgI3NldHRpbmcgdG8gRkVxdWFydGljIHNsb3dzIHRyYW5zaWVudCBzaW11bGF0aW9ucyB3YXkgdG9vIG11Y2gNCnNldC5kaXNjcmV0aXNhdGlvbihVbnNhdFN0YXQsbm9kZXMsIkZFcXVhZHJhdGljIikgDQojIEFuIGluaXRpYWxpc2F0aW9uIGlzIGNlcnRhaW5seSBpbXBvcnRhbnQgKGFzIHRoZSBkZWZhdWx0IGluaXRpYWxpc2F0aW9uIGJ5IHplcm9zIGRvZXMgbm90IG1ha2Ugc2Vuc2UgaGVyZSkuIA0KIyBTbyB3ZSBpbml0aWFsaXplIHdpdGggYSBuby1mbHV4IHByb2ZpbGU6DQpQc2kuZXF1aSA9IGFwcHJveGZ1bihjKGRvbWFpblsxXSxkb21haW5bMl0pLGMoZG9tYWluWzJdLGRvbWFpblsxXSkpDQpkby5pbml0aWFsaXplKFVuc2F0U3RhdCxQc2kuZXF1aSkNCg0KIyBXZSB3YW50IHRvIGF2b2lkIHByZXNzdXJlcyB0byBiZSB0b28gbmVnYXRpdmUgKDwtMUUtNSkNClBzaS5pbi5yYW5nZSA9IGZ1bmN0aW9uKHosUHNpKSAjIGl0IGhhcyAyIGFyZ3VtZW50czsgdGhlIHBvc2l0aW9uIGFuZCB0aGUgc3RhdGUuIFNldCB0aGUgbWluLiB2YWx1ZSB0byBwRjQuMiBub3QgcEY1DQp7DQogIyBpZihQc2k8PS0xKjEwXjQuMikgcmV0dXJuKEZBTFNFKQ0KICBpZihQc2kgPD0gLTE2MDAwKSByZXR1cm4oRkFMU0UpICMgYWxzbyBpbiB0aGUgRmVkZGVzIGZ1bmN0aW9uIC0xNjAwMCB3aGljaCBpcyA0LjIwNDEyIGFuZCBub3QgNC4yDQogICAjaWYoUHNpPD0tMSoxMF41KSByZXR1cm4oRkFMU0UpDQogIGVsc2UgcmV0dXJuKFRSVUUpDQp9DQpzZXQuaXNhY2NlcHRhYmxlKFVuc2F0U3RhdCxQc2kuaW4ucmFuZ2UpICN0aGUgc2V0LmlzYWNjZXB0YWJsZSBmdW5jdGlvbiBsaW1pdHMgdGhlIHBvc3NpYmxlIHZhbHVlcyB0aGF0IHRoZSBzdGF0ZSBjYW4gYmVjb21lDQoNCiMjIyMjIyMjIyMjIyMjIEJvdW5kYXJ5IGNvbmRpdGlvbnMgIyMjIyMjIyMjIyMjIyMjIyMjDQpzZXQuQkMuZml4ZWRzdGF0ZShVbnNhdFN0YXQsImxlZnQiLDAuMCkNCiMgQXQgdGhlIGJvdHRvbSBvZiB0aGUgbW9kZWwsIHdoZXJlIHRoZSBncm91bmR3YXRlciBsZXZlbCBpcyBsb2NhdGVkLCBib3RoIHByZXNzdXJlIGhlYWQgYW5kIGVsZXZhdGlvbiBoZWFkIGFyZSB6ZXJvLiANCg0KcHB0ID0gMC4xDQpzZXQuQkMuZml4ZWRmbHV4KFVuc2F0U3RhdCwicmlnaHQiLCJwcHQiKSAgDQojIDAuMWNtL2Q9MW1tL2QuIEFsdGhvdWdoIHRoaXMgZmxvdyBpcyBkaXJlY3RlZCBkb3dud2FyZHMsIHdlIGRlZmluZSBpdCBwb3NpdGl2ZSBiZWNhdXNlIGl0IGdvZXMgaW50byB0aGUgZG9tYWluLg0KDQoNCiMjIyMjIyMjIGNoZWNraW5nIGFuZCBzb2x2aW5nIHRoZSBtb2RlbCAjIyMjIyMjIyMjIyMjIyMjIyMjIw0Kc3VtbWFyeShVbnNhdFN0YXQpDQpzb2x2ZS5zdGVwcyhVbnNhdFN0YXQsdmVyYm9zZWxldmVsPTEpDQpwbG90KFVuc2F0U3RhdCxmbHV4cGxvdCA9IFQsdmVydGljYWwgPSBUKQ0KDQpgYGANCg0KIyMgQWRkaXRpb25hbCBmbG93IGZ1bmN0aW9ucw0KDQpJbiBhZGRpdGlvbiB0aGVyZSB3aWxsIGJlIGZ1bmN0aW9ucyBmb3IgcG9uZGluZywgdGhlIGxvd2VyIGJvdW5kYXJ5IGNvbmRpdGlvbnMgKHNpbXVsYXRpbmcgdGhlIGRyYWluYWdlIGVmZmVjdCkgYW5kIHRyYW5zcGlyYXRpb24gb2Ygcm9vdHMgKEZlZGRlcyBmdW5jdGlvbikuDQpUaGUgbG93ZXIgYm91bmRhcnkgY29uZGl0aW9uIGlzIGEgQ2F1Y2h5IHR5cGUgQkM6DQokJA0KUV97ZHJhaW5nZX09IFxmcmFje0hfe2Rybi5sZXZlbH0gLSBIX3tib3R0b219fXtDfQ0KJCQNCg0KYGBge3J9DQpkcm4ucmVzaXN0YW5jZSA9IDUwICN1bml0IGlzIGluIGRheXMNCmRybi5sZXZlbCA9IC0xMjAgIyBzbyBYWCBjbSBiZWxvdyBzdXJmYWNlIGxldmVsDQpkcm4uZnVuID0gZnVuY3Rpb24ocHNpKQ0Kew0KIyAgaWYoKGRybi5sZXZlbCAtIChwc2kgKyBib3R0b20pKSA8IDApIHJldHVybigoZHJuLmxldmVsIC0gKHBzaSArIGJvdHRvbSkpL2Rybi5yZXNpc3RhbmNlKQ0KIyAgZWxzZSByZXR1cm4oMCkNCiAgcmV0dXJuKChkcm4ubGV2ZWwgLSAocHNpICsgYm90dG9tKSkvZHJuLnJlc2lzdGFuY2UpDQp9DQoNCg0KYGBgDQoNCiMjIyBSZXBsYWNpbmcgdGhlIGJvdHRvbSBmaXhlZCBoZWFkIGJvdW5kYXJ5IGNvbmRpdGlvbiB3aXRoIHRoZSBkcmFpbmFnZSB0eXBlDQoNCmBgYHtyfQ0Kc2V0LkJDLmZsdXhzdGF0ZShVbnNhdFN0YXQsImxlZnQiLCBkcm4uZnVuKQ0Kc3VtbWFyeShVbnNhdFN0YXQpDQoNCmRybi5wc2kuZXF1aSA9IGFwcHJveGZ1bihjKGRvbWFpblsxXSxkb21haW5bMl0pLGMoKGRybi5sZXZlbCAtIGJvdHRvbSksZHJuLmxldmVsKSkNCiAgICAgICAgICAgICAgICAgICAgICAgICANCmRvLmluaXRpYWxpemUoVW5zYXRTdGF0LGRybi5wc2kuZXF1aSkNCnNvbHZlLnN0ZXBzKFVuc2F0U3RhdCx2ZXJib3NlbGV2ZWwgPSAxKQ0KcGxvdChVbnNhdFN0YXQsZmx1eHBsb3QgPSBULHZlcnRpY2FsID0gVCkNCg0KIyN0cnlpbmcgdGhlIG5ldyBib3R0b20gYm91bmRhcnkgY29uZGl0aW9uIHdpdGggYSB2aWV3IGRpZmZlcmVudCBwcHQncw0KcHB0ID0gMC4yDQpzb2x2ZS5zdGVwcyhVbnNhdFN0YXQsdmVyYm9zZWxldmVsID0gMSkNCnBsb3QoVW5zYXRTdGF0LGZsdXhwbG90ID0gVCx2ZXJ0aWNhbCA9IFQpDQoNCiNwcHQgPSAtMC4yDQojIHNvbHZlLnN0ZXBzKFVuc2F0U3RhdCx2ZXJib3NlbGV2ZWwgPSAxKQ0KIyBwbG90KFVuc2F0U3RhdCxmbHV4cGxvdCA9IFQsdmVydGljYWwgPSBUKQ0KIyBkYXRhZnJhbWUuYmFsYW5jZShVbnNhdFN0YXQpDQojIGRhdGFmcmFtZS5ib3VuZGFyaWVzKFVuc2F0U3RhdCkNCg0KIyMjcGxvdHRpbmcgZnVuY3Rpb25zIHNob3dpbmcgdGhlIHdhdGVyIHRhYmxlIGluIHRoZSBwcm9maWxlLA0KIyMjIHRoZSBtb2lzdHVyZSBkaXN0cmlidXRpb24gYW5kIHRyYW5zcGlyYXRpb24gdGhyb3VnaCB0aGUgcm9vdHMNCg0KcGxvdC53dGFibGUgPSBmdW5jdGlvbihtb2RlbCkNCnsNCiAgZGYgPSBkYXRhZnJhbWUuc3RhdGVzKG1vZGVsKQ0KICBwbG90KGRmJHN0YXRlLCBkZiR4LCB0eXBlID0gIm8iLGNvbCA9ICJibHVlIixsd2Q9Mi41LA0KICAgICAgIG1haW4gPSAiUHNpKHopLCBsb2NhdGlvbiB3YXRlciB0YWJsZSIsIHhsYWIgPSAicHNpIGluIGNtIiwgeWxhYiA9ICJkZXB0aCBwcm9maWxlIikNCiAgYWJsaW5lKGggPSBkcm4ubGV2ZWwsIGNvbCA9ICJicm93biIsIGx3ZD0gMywgbHR5ID0gImRhc2hlZCIpDQogIHRleHQoMjAsZHJuLmxldmVsLCAiZHJhaW5hZ2UgbGV2ZWwiLCBjb2wgPSAiYnJvd24iKQ0KICBwc2kuZnVuID0gYXBwcm94ZnVuKGRmJHN0YXRlLGRmJHgpDQogIGFibGluZShoID0gcHNpLmZ1bigwKSxjb2wgPSAicmVkIiwgbHdkID0gMikNCiAgdGV4dCgyMCxwc2kuZnVuKDApLCJ3YXRlciB0YWJsZSIgLCBjb2wgPSAicmVkIikNCiAgZ3JpZCgpDQp9DQoNCnBsb3QudGhldGEgPSBmdW5jdGlvbihtb2RlbCkNCnsNCiAgZGYgPSBkYXRhZnJhbWUuc3RhdGVzKG1vZGVsKQ0KICB0aGV0YS5wc2kgPSBzb2lsdHlwZSR0aGV0YS5mdW4oZGYkc3RhdGUpDQogIHBsb3QodGhldGEucHNpLGRmJHgsIHR5cGUgPSAibyIsIGNvbCA9ICJncmVlbiIsIGx3ZCA9IDIuNSwNCiAgICAgICBtYWluID0gIlRoZXRhKHopLCBtb2lzdHVyZSBjb250ZW50IGluIHByb2ZpbGUiLCB4bGFiID0gInRoZXRhIiwgeWxhYiA9ICJkZXB0aCBwcm9maWxlIikNCiAgYWJsaW5lKGggPSBkcm4ubGV2ZWwsIGNvbCA9ICJicm93biIsIGx3ZD0gMywgbHR5ID0gImRhc2hlZCIpDQogIHRleHQoKG1lYW4odGhldGEucHNpKSksZHJuLmxldmVsLCAiZHJhaW5hZ2UgbGV2ZWwiLCBjb2wgPSAiYnJvd24iKQ0KICBncmlkKCkNCn0NCg0KDQpwbG90Lnd0YWJsZShVbnNhdFN0YXQpDQpwbG90LnRoZXRhKFVuc2F0U3RhdCkNCg0KYGBgDQojIyBBZGRpbmcgdGhlIHZlZ2V0YXRpb24gdG8gdGhlIG1vZGVsDQoNCkV2YXBvdHJhbnNwaXJhdGlvbiBpcyBhbiBpbXBvcnRhbnQgZW52aXJvbm1lbnRhbCBjb25kaXRpb24gZm9yIHByb2Nlc3NlcyBuZWFyIHRoZSBzb2lsIHN1cmZhY2UsIHN1Y2ggYXMgdGhlIGRldmVsb3BtZW50IG9mIHZlZ2V0YXRpb24uIEV2YXBvdHJhbnNwaXJhdGlvbiBkZXBlbmRzIG9uOiAgICAgIA0KDQoqICAgVGhlIHBvdGVudGlhbCBldmFwb3RyYW5zcGlyYXRpb24gcmF0ZSwgd2hpY2ggaXMgZGV0ZXJtaW5lZCBieSB3ZWF0aGVyIGNvbmRpdGlvbnM7ICAgIA0KKiAgIFRoZSByb290aW5nIGRpc3RyaWJ1dGlvbiBvZiB0aGUgdmVnZXRhdGlvbjsgICANCiogICBUaGUgYWN0dWFsIHNvaWwgbW9pc3R1cmUgY29udGVudFxwcmVzc3VyZSBpbiB0aGUgcm9vdGluZyB6b25lLiAgIA0KDQpUaGUgdG90YWwgKnBvdGVudGlhbCogZXZhcG90cmFuc3BpcmF0aW9uIHJhdGUgb2YgdGhlIHZlZ2V0YXRpb24gaXMgZGV0ZXJtaW5lZCBieSB3ZWF0aGVyIGNvbmRpdGlvbnMgYW5kIHdpbGwgYmUgZ2l2ZW4uDQoNClRoZSByb290aW5nIGRpc3RyaWJ1dGlvbiBtYXkgYmUgZGVzY3JpYmVkIGJ5IGEgZnVuY3Rpb24gYHJvb3QuZGlzdHIuZnVuKHopYCB3aGljaCBnaXZlcyB0aGUgZnJhY3Rpb24gb2Ygcm9vdHMgcHJlc2VudCBhdCBlYWNoIGRlcHRoIChnaXZpbmcgMSB3aGVuIGludGVncmF0ZWQgd2l0aCBkZXB0aCkuIA0KVGhlIHJvb3RpbmcgZGVwdGggZm9yIGVhY2ggdmVnZXRhdGlvbiB0eXBlIGNhbiBiZSBkZXNjcmliZWQgaW4gbWFueSB3YXlzLiANCkhlcmUsIGF0IHRoaXMgc3RhZ2Ugd2Ugd2lsbCBhc3N1bWUgZ3Jhc3Mgd2l0aCA3NSUgcm9vdGluZyBpbiB0aGUgZmlyc3QgMjAgY20NCmFuZCAyNSUgdGlsbCAxMTAgY20gYmVsb3cgc3VyZmFjZS4gKEp1c3QgbXkgZ3Vlc3MpLg0KDQpgYGB7cn0NCnJvb3QucGFydGl0aW9uID0gYygwLjAsLTIwLjAsLTIwLjAxLC0xMTAuMCwtMTEwLjAxKQ0Kcm9vdC5yYXRlcy5mcmFjdGlvbnMgPSBjKDAuNzUvMjAsMC43NS8yMCwwLjI1LzkwLDAuMjUvOTAsMC4wKQ0Kcm9vdC5kaXN0ci5mdW4gPSBhcHByb3hmdW4ocm9vdC5wYXJ0aXRpb24scm9vdC5yYXRlcy5mcmFjdGlvbnMscnVsZSA9IDIpDQpwbG90KHJvb3QucGFydGl0aW9uLHJvb3QuZGlzdHIuZnVuKHJvb3QucGFydGl0aW9uKSwgdHlwZSA9ICJvIiwNCiAgIGNvbCA9ICJncmVlbiIsIGx3ZCA9IDIsIG1haW4gPSAiUm9vdCBkaXN0cmlidXRpb24gZnVuY3Rpb24iLA0KICAgeWxhYiA9ICJGcmFjdGlvbiBvZiB0aGUgcm9vdCBkaXN0cmlidXRpb24gKC0pIiwgeGxhYiA9ICJSb290aW5nIGRlcHRoIChjbSkgcmVsLiB0byBzdXJmYWNlIGxldmVsIikNCmdyaWQoKQ0KDQojIyNjaGVjayBmb3IgdGhlIGludGVncmFsDQoyMCowLjc1LzIwICsgOTAqMC4yNS85MA0Kei5yb290ID0gc2VxKC0xMDkuNSwtMC41LGJ5ID0gMSkNCnN1bShyb290LmRpc3RyLmZ1bih6LnJvb3QpKQ0KDQoNCmBgYA0KDQoNCg0KVGhlIHZlZ2V0YXRpb24gY2FuIG9ubHkgdHJhbnNwaXJlIG9wdGltYWxseSB3aGVuIHRoZXJlIGlzIHN1ZmZpY2llbnQgbW9pc3R1cmUgaW4gdGhlIHNvaWwgd2hpY2ggaXMgJ2Vhc2lseScgYXZhaWxhYmxlIGZvciB0aGUgdmVnZXRhdGlvbi4gVGhlIGVmZmVjdCBvZiBzb2lsIG1vaXN0dXJlIGNhbiBiZSBzaW11bGF0ZWQgd2l0aCB0aGUgc28gY2FsbGVkIEZlZGRlcyAocmV0aXJlZCBQcm9mLiBmcm9tIFdVUikgZnVuY3Rpb24uICANCg0KVGhlIGNvZGUgY2h1bmsgYmVsb3cgaWxsdXN0cmF0ZXMgc3VjaCBhIEZlZGRlcyBmdW5jdGlvbi4gVGhpcyBmdW5jdGlvbiByZWxhdGVzIHRoZSByb290IHdhdGVyIHVwdGFrZSBhdCBhIGNlcnRhaW4gZGVwdGggd2l0aCB0aGUgcHJlc3N1cmUgKHN1Y3Rpb24pIGF0IHRoYXQgZGVwdGguIFdoZW4gdGhlIHNvaWwgaXMgdG9vIHdldCBvciB0b28gZHJ5LCBwbGFudCB0cmFuc3BpcmF0aW9uIGFuZCB0aHVzIHJvb3QgdXB0YWtlIGFyZSByZWR1Y2VkIGJlbG93IHRoZWlyIHBvdGVudGlhbCByYXRlcywgYW5kIHRoZSBGZWRkZXMgcmVkdWN0aW9uIGZhY3RvciBpcyBzbWFsbGVyIHRoYW4gMS4gDQoNCmBgYHtyfQ0Kb3B0aW9ucyhzY2lwZW4gPSA5OTkpICNhdm9pZCBzY2llbnRpZmljIG5vdGF0aW9uICh3aXRoIEUgcG93ZXIpIGluIHRoZSBwbG90DQojRmVkZGVzLnBzaWFicyA9IGMoMC4xLDEwMCwxMDAwLDE2MDAwKSAjdHlwaWNhbCBwc2kgdmFsdWVzIHRvIGRlZmluZSBGZWRkZXMgZnVuY3Rpb24uIE5vdGUgdGhlc2UgYXJlIGFic29sdXRlIHBzaSdzIGR1ZSB0byB0aGUgcGxvdC4gTk9UIGluIHRoZSBmdW5jdGlvbg0KRmVkZGVzLnBzaWFicyA9IGMoMC4xLDEwMCwxMDAwLDE1MDAwKQ0KRmVkZGVzLnJlZCA9IGMoMCwxLDEsMCkNCkZlZGRlcy5mdW49IGFwcHJveGZ1bigtRmVkZGVzLnBzaWFicyxGZWRkZXMucmVkLHJ1bGU9MikNCnBsb3QoRmVkZGVzLnBzaWFicyxGZWRkZXMuZnVuKC1GZWRkZXMucHNpYWJzKSx0eXBlPSJsIixsd2Q9Myx4bGltPXJldihyYW5nZShGZWRkZXMucHNpYWJzKSksbG9nPSJ4Iix5bGFiPSJSb290IHVwdGFrZSByZWR1Y3Rpb24iLHhsYWI9IlN1Y3Rpb24gKHxwc2l8KSIpDQphYmxpbmUodj1GZWRkZXMucHNpYWJzLGNvbD0icmVkIikNCnRleHQoMSwwLjksIlRvbyB3ZXQiLCBjb2w9InJlZCIpDQp0ZXh0KDMwMCwwLjcsIk9wdGltYWwiLCBjb2w9InJlZCIpDQp0ZXh0KDUwMDAsMC45LCJUb28gZHJ5Iixjb2w9InJlZCIpDQpncmlkKCkNCmBgYA0KDQoNCiMjIyBBZGRpbmcgdGhpcyBzcGF0aWFsIGZsdXggdG8gdGhlIG1vZGVsDQoNClRoZSBlZmZlY3Qgb2YgdGhlIHZlZ2V0YXRpb24sIHRvb3RzIGV4dHJhY3Rpbmcgd2F0ZXIgZnJvbSB0aGUgc29pbCBmb3IgdHJhbnNwaXJhdGlvbiBkZXBlbmRpbmcgb24gdGhlIHJvb3QgZGlzdHJpYnV0aW9uLCBzdWN0aW9uICgkXHBzaSQpIGFuZCB0aGUgcG90ZW50aWFsIGV2YXBvdHJhbnBzaXJhdGlvbi4NCg0KVGhyZWUgYXNwZWN0cyBhcmUgY29uc2lkZXJlZDogIA0KDQoxLiAgUm9vdGluZyBkaXN0cmlidXRpb24gYHJvb3QuZGlzdHIuZnVuYCAgDQoyLiAgRmVkZGVzIGZ1bmN0aW9uIGBGZWRkZXMuZnVuYCAgDQozLiAgUG90ZW50aWFsIGV2YXBvdHJhbnNwaXJhdGlvbiBgRXBvdGANCg0KVGhlIGZ1bmN0aW9uIHRvIGNhbGN1bGF0ZSB0aGlzIGlzIDogYEVhY3QgPSByb290LmRpc3RyLmZ1bih6KSpGZWRkZXMuZnVuKHBzaSkqRXBvdGANCg0KYGBge3J9DQoNCkVhY3QgPSBmdW5jdGlvbih6LHBzaSkNCnsNCiAgcmV0dXJuKC1yb290LmRpc3RyLmZ1bih6KSAqIEZlZGRlcy5mdW4ocHNpKSAqIEVwb3QpDQp9DQoNCmFkZC5zcGF0aWFsZmx1eChVbnNhdFN0YXQsIHJhdGUgPSBFYWN0LCBuYW1lID0gIlRyYW5zcGlyYXRpb24iKQ0KDQpgYGANCg0KDQojIyBjdXJyZW50IG1vZGVsDQpUaGUgY3VycmVudCBtb2RlbCBjb25zaXN0cyBvZjogIA0KDQoqICBwcmVjaXBpdGF0aW9uIGFwcGxpZWQgdG8gdGhlIHVwcGVyIGJvdW5kYXJ5IGFzIGEgcHJlc2NyaWJlZCBmbHV4IGBwcHRgICANCiogIHBvdGVudGlhbCBldmFwb3RyYW5zcGlyYXRpb24gYGVwb3RgICBhcyBhIGZsdXggcmF0ZSB0byBkZXRlcm1pbmUgRWFjdCAgDQoqICBhY3R1YWwgZXZhcG90cmFuc3BpcmF0aW9uIGBFYWN0YCBhcyBhIHNwYXRpYWwgZmx1eCAoYWxvbmcgdGhlIHJvb3RzKSAgDQoqICBDYXVjaHkgdHlwZSBib3VuZGFyeSBjb25kaXRpb24gYXQgdGhlIGxvd2VyIGJvdW5kYXJ5IHNpbXVsYXRpbmcgdGhlIGVmZmVjdCBvZiBhIGxvY2FsIGRyYWluYWdlIHN5c3RlbSAgDQoNCg0KYGBge3J9DQoNCnBsb3Qucm9vdHMgPSBmdW5jdGlvbihtb2RlbCkNCnsNCiAgZGYgPSBkYXRhZnJhbWUuZXh0ZXJuYWxmbHV4ZXMobW9kZWwpDQogIHBsb3QoZGYkYFN7VHJhbnNwaXJhdGlvbn1keGAsZGYkeCwgdHlwZSA9ICJvIiwgY29sID0gImdyZWVuIiwgbHdkID0gMi41LA0KICAgICAgIG1haW4gPSAiTW9pc3R1cmUgZXh0cmFjdGlvbiBieSByb290cyBpbiBkZXB0aCIsIHhsYWIgPSAicmF0ZSAoY20vZCkiLCB5bGFiID0gImRlcHRoIHByb2ZpbGUiKQ0KICBhYmxpbmUoaCA9IGRybi5sZXZlbCwgY29sID0gImJyb3duIiwgbHdkID0gMywgbHR5ID0gImRhc2hlZCIpDQogIHRleHQobWVhbihkZiRgU3tUcmFuc3BpcmF0aW9ufWR4YCksZHJuLmxldmVsLCJkcmFpbmFnZSBsZXZlbCIsIGNvbCA9ICJicm93biIpDQogIGdyaWQoKQ0KfQ0KDQpzdW1tYXJ5KFVuc2F0U3RhdCkNCg0KIyMjI3Rlc3RpbmcgdGhlIG1vZGVsIHdpdGggbWV0ZW8gZGF0YQ0KcHB0ID0gMC4xMjQNCkVwb3QgPSAwLjA1IzAuMg0KZG8uaW5pdGlhbGl6ZShVbnNhdFN0YXQsIGRybi5wc2kuZXF1aSkNCnNvbHZlLnN0ZXBzKFVuc2F0U3RhdCx2ZXJib3NlbGV2ZWwgPSAxKQ0KDQpwbG90KFVuc2F0U3RhdCxmbHV4cGxvdCA9IFQpDQpwbG90Lnd0YWJsZShVbnNhdFN0YXQpDQpwbG90LnRoZXRhKFVuc2F0U3RhdCkNCnBsb3Qucm9vdHMoVW5zYXRTdGF0KQ0KYmFsLmFzZml4ZWRmbHV4ID0gZGF0YWZyYW1lLmJhbGFuY2UoVW5zYXRTdGF0KQ0KZXh0ZXIuYXNmaXhlZGZsdXggPSBkYXRhZnJhbWUuZXh0ZXJuYWxmbHV4ZXMoVW5zYXRTdGF0KQ0KYGBgDQoNCiMjIEluIGNhc2Ugb2YgcG9kaW5nDQoNCldoZW4gdGhlIGFtb3VudCBvZiBwcmVjaXBpdGF0aW9uIGlzIGxhcmdlciB0aGFuIHRoZSBzb2lsIGNhbiBhYnNvcmIsIHBvbmRpbmcgKG5sOiBwbGFzdm9ybWluZykgd2lsbCB0YWtlIHBsYWNlIA0KVG8gc2ltdWxhdGUgdGhpcywgYW4gYWRkaXRpb25hbCBDYXVjaHkgdHlwZSBib3VuZGFyeSBjb25kaXRpb24gd2lsbCBiZSBhcHBsaWVkIHRvIHRoZSB1cHBlciBub2RlLg0KDQokJA0KIHFfe3BvbmRpbmd9ID0gXGZyYWN7XHBzaV97c3VyZmFjZX0tMH17Q197dG9wfX0NCiQkDQpJdCB3aWxsIG9ubHkgYmUgYWN0aXZhdGVkIHdoZW4gdGhlIHByZXNzdXJlIGF0IHRoZSBzdXJmYWNlICh1cHBlciBub2RlKSBpcyBsYXJnZXIgdGhhbiAwLjAgY20NCg0KYGBge3J9DQpjLnRvcCA9IDUwICNkYXlzDQpwb25kaW5nID0gZnVuY3Rpb24ocHNpKQ0Kew0KICBpZiAocHNpID4gMCkNCiAgew0KICAgIHJldHVybigtcHNpL2MudG9wKSAgDQogIH1lbHNlew0KICAgIHJldHVybigwKQ0KICB9DQp9DQpgYGANCg0KIyMgQWRqdXN0aW5nIHRoZSBtb2RlbCANCldlIG5lZWQgdG8gcmVwbGFjZSB0aGUgYHNldC5CQy5maXhlZGZsdXhgIHVwcGVyIGJvdW5kYXJ5IHRvIGFuIGBhZGQucG9pbnRmbHV4YCBmdW5jdGlvbiBhbmQgYXBwbHkgdGhlIHBvbmRpbmcgZnVuY3Rpb24gdG8gdGhlIHVwcGVyIG5vZGUgd2l0aCBgc2V0LkJDLmZsdXhzdGF0ZWAuDQoNCmBgYHtyfQ0Kc2V0LkJDLmZsdXhzdGF0ZShtb2RlbCA9IFVuc2F0U3RhdCwgd2hlcmUgPSAicmlnaHQiLCBmdW5jID0gcG9uZGluZykNCmFkZC5wb2ludGZsdXgobW9kZWwgPSBVbnNhdFN0YXQsIGF0ID0gMCwgdmFsdWUgPSAicHB0IiwgInByZWNpcGl0YXRpb24iKQ0Kc3VtbWFyeShVbnNhdFN0YXQpDQpzb2x2ZS5zdGVwcyhVbnNhdFN0YXQsIHZlcmJvc2VsZXZlbCA9IDEpDQpwbG90Lnd0YWJsZShVbnNhdFN0YXQpDQpwbG90LnRoZXRhKFVuc2F0U3RhdCkNCmBgYA0KDQoNCg0KDQoNCiMjIFRoZSB0cmFuc2llbnQgbW9kZWwNCg0KV2hlbiB0aGUgZHluYW1pY3MgYXJlIGludHJvZHVjZWQgaW50byB0aGUgbW9kZWwgdGltZXMgc2VyaWVzIG9mIGRpZmZlcmVudCBhc3BlY3RzIGNvdWxkIGJlIGFuYWx5emVkLg0KSW50cm9kdWNpbmcgZHluYW1pY3MgbWVhbnMgYWN0dWFsbHkgdGhhdCBjaGFuZ2VzIGluIHRpbWUgYXJlIGNvbnNpZGVyZWQuIFRoaXMgaW1wYWN0cyBpbiB0aGUgZmlyc3QgcGxhY2UgdGhlIGFtb3VudCBvZiB3YXRlciBzdG9yZWQgaW4gdGhlIG1vZGVsOiAkXGZyYWN7XHBhcnRpYWwgXHRoZXRhfXtccGFydGlhbCB0fSQuDQoNClRoZSBiYXNpYyBlcXVhdGlvbiB0byBiZSBzb2x2ZWQgKFJpY2hhcmQncyBlcXVhdGlvbikgaXM6DQokJA0KXGZyYWN7XHBhcnRpYWwgXHRoZXRhKFxwc2kpfXtccGFydGlhbCB0fSA9IC1rKFxwc2kpXGxlZnQgW1xmcmFje1xwYXJ0aWFsIFxwc2l9e1xwYXJ0aWFsIHp9ICsgMSBccmlnaHRdDQokJA0KDQpUaGUgbW9pc3R1cmUgKHN0b3JhZ2UpIGNoYW5nZSBvdmVyIHRpbWUgKDogJFxmcmFje1xwYXJ0aWFsIFx0aGV0YShccHNpKX17XHBhcnRpYWwgdH0kKSBpcyBpbXBsZW1lbnRlZCBhcyBhIGBzcGF0aWFsIGZsdXhgIG9mIHRpbWUgYW5kIHNwYWNlDQpTaW5jZSB3ZSBhcmUgY2FsY3VsYXRpbmcgJFxwc2kkIGFzIHRoZSBzdGF0ZSBvZiB0aGUgbW9kZWwsIHdlIGFyZSBnb2luZyB0byB1c2UgdGhlIGBzb2lsdHlwZSR0aGV0YS5mdW4ocHNpKWAgZnVuY3Rpb24gdG8gZGV0ZXJtaW5lIHRoZSBhbW91bnQgb2YgbW9pc3R1cmUgY29udGVudCBiYXNlZCBvbiB0aGUgY3VycmVudCBwcmVzc3VyZSBoZWFkICgkXHBzaSQpDQoNClRoaXMgY2FuIGJlIGltcGxlbWVudGVkIG51bWVyaWNhbGx5IGFzOg0KDQokJA0KUV97c3RvcmFnZX0gPSBcZnJhY3tcdGhldGEgXGxlZnQgKFxwc2lee3QrXERlbHRhIHR9XHJpZ2h0ICkgLSBcdGhldGEgXGxlZnQoIFxwc2ledCBccmlnaHQgKSB9e1xEZWx0YSB0fVxEZWx0YSB6DQokJA0KVGhlIG5ldy9jdXJyZW50IHByZXNzdXJlIGhlYWQgaXMgdGhlIG9uZSB3aGljaCBpcyBldmFsdWF0ZWQgYXQgdGhlIG5ldyB0aW1lIHN0ZXA6ICRccHNpXnt0K1xEZWx0YSB0fSQgYnV0IHNpbmNlIHdlIGFsc28gbmVlZCB0byBjb25zaWRlciB0aGUgcHJldmlvdXMgcHJlc3N1cmUgaGVhZCAkXHBzaV50JCwgd2Ugd2lsbCBzYXZlIHRoZXNlIHByZXNzdXJlIGp1c3QgYmVmb3JlIHdlIHN0YXJ0IHdpdGggYSBuZXcgdGltZSBzdGVwIGFuZCBjYWxsIHRoaXMgYHBzaS5vbGRgLiAgDQpXaGVuIG1vcmUgd2F0ZXIgaXMgc3RvcmVkIGluIHRoZSBtb2RlbCwgaXQgbWVhbnMgdGhhdCBsZXNzIGlzIGF2YWlsYWJsZSBmb3IgZmxvdy4gVGhpcyBtZWFucyB0aGF0IGEgbWludXMgc2lnbiBzaG91bGQgYmUgaW50cm9kdWNlZC4gVGhpcyBjb3VsZCBlYXNpbHkgYmUgYWNjb21wbGlzaGVkIGJ5IHN3aXRjaGluZyBib3RoIHByZXNzdXJlIGhlYWRzIGluIHRoZSBhZm9yZW1lbnRpb25lZCBlcXVhdGlvbi4NCg0KYGBge3J9DQpzdG9yYWdlID0gZnVuY3Rpb24oeixzdGF0ZSkNCnsNCiAgc3RvcmFnZSA9IChzb2lsdHlwZSR0aGV0YS5mdW4ocHNpLm9sZCh6KSkgLSBzb2lsdHlwZSR0aGV0YS5mdW4oc3RhdGUpKS9kZWx0YS50DQp9DQpgYGANCg0KQmUgYXdhcmUgdGhlIGBzdG9yYWdlYCBmdW5jdGlvbiB3aWxsIGJlIG11bHRpcGxpZWQgd2l0aCB0aGUgbm9kYWwgZGlzdGFuY2UgaW50ZXJuYWxseSB0byBjb21lIHRvIHRoZSB1bml0IG9mIEwvVCAoY20vZCkuDQoNCiMjIyBTZXR0aW5nIHVwIHRoZSB0cmFuc2llbnQgbW9kZWwNCg0KQSBjb3B5IG9mIHRoZSBzdGF0aW9uYXJ5IG1vZGVsIGlzIHVzZWQgdG8gcXVpY2tseSBzZXQgdXAgdGhlICB0cmFuc2llbnQgbW9kZWwgDQoNCmBgYHtyfQ0KVW5zYXRUcmFucyA9IGNvcHkubW9kZWwoVW5zYXRTdGF0KQ0KYWRkLnNwYXRpYWxmbHV4KFVuc2F0VHJhbnMsIHJhdGUgPSBzdG9yYWdlLCBuYW1lID0gIlN0b3JhZ2UiKQ0KZG8uaW5pdGlhbGl6ZShVbnNhdFRyYW5zLCBpbml0ID0gZHJuLnBzaS5lcXVpKQ0KcHNpLm9sZCA9IHN0YXRlLmZ1bihVbnNhdFRyYW5zKSAjc2V0dGluZyB0aGUgc3RyYXRpbmcgc3RhdGUgb2YgdGhlIHRyYW5zaWVudCBzaW11bGF0aW9uIHRvIGRldGVybWluZSB0aGUgZmlyc3Qgc3RvcmFnZSBjaGFuZ2UNCmBgYA0KDQpDcmVhdGluZyBhIHRpbWUgbG9vcCBmb3IgdGhlIGNhbGN1bGF0aW9ucw0KDQpBIGZ1bmN0aW9uIHRvIHdyaXRlIHRoZSByZXF1aXJlZCBvdXQgcHV0IGluIGRpZmZlcmVudCBkYXRhLmZyYW1lcyBmb3IgYW5hbHlzaXMgYW5kIGluc3BlY3Rpb24gbGF0ZXIgb24uDQoNCmBgYHtyfQ0Kc2F2ZS5kYXRhID0gZnVuY3Rpb24obW9kZWwgPSBtb2RlbG5hbWUpDQp7DQogIGRmLnN0YXRlID0gZGF0YWZyYW1lLnN0YXRlcyhtb2RlbCkNCiAgc3RhdGUudG1wIDw8LSByYmluZChzdGF0ZS50bXAsZGYuc3RhdGUkc3RhdGUpICNhZGRpbmcgdGhlIHN0YXRlcyAocHNpKSB0byB0aGlzIGRhdGEuZnJhbWUNCiAgZGYuZXh0ZXJuID0gZGF0YWZyYW1lLmV4dGVybmFsZmx1eGVzKG1vZGVsKQ0KICB0cmFuc3AudG1wIDw8LSByYmluZCh0cmFuc3AudG1wLGRmLmV4dGVybiRgU3tUcmFuc3BpcmF0aW9ufWR4YCkgI2ZvciB0cmFuc3BpcmF0aW9uKHJvb3RzKSBkYXRhDQogIG1vaXN0LnRtcCA8PC0gcmJpbmQobW9pc3QudG1wLGRmLmV4dGVybiRgU3tTdG9yYWdlfWR4YCkNCiAgI2JhbGFuY2UgdGVybXMgdHJhbnNwOlsyLDNdLGZsb3cyc3RvOlszLDJdLHN0bzJmbG93WzMsM10NCiAgZGYuYmFsYW5jZSA9IGRhdGFmcmFtZS5iYWxhbmNlKG1vZGVsKQ0KICAjdGhlIGRyYWluYWdlL3Vwd2FyZCBzZWVwYWdlIGF0IHRoZSBib3R0b20gd2lsbCBiZSB3cml0dGVuIGluIGJhbGFuY2UuZGF0DQogIGRmLmJvdW5kID0gZGF0YWZyYW1lLmJvdW5kYXJpZXMobW9kZWwpDQogIGJhbGFuY2UudG1wIDw8LSByYmluZChiYWxhbmNlLnRtcCxjKHByZWNpcFtjdXJyZW50LnRpbWVdLGV2YXBvdFtjdXJyZW50LnRpbWVdLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkZi5iYWxhbmNlWzIsM10sZGYuYmFsYW5jZVszLDJdLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkZi5iYWxhbmNlWzMsM10sZGYuYm91bmRbMSw0XSxkZi5ib3VuZFsyLDRdKSkNCn0NCmBgYA0KDQoNClJlYWRpbmcgc29tZSBtZXRlbyBkYXRhIGZyb20gS05NSSwgcHJlY2lwaXRhdGlvbiBhbmQgcG90ZW50aWFsIGV2YXBvdHJhbnNwaXJhdGlvbg0KSW4gdGhlIGZpbGUgaXQgcmVhZHM6DQpSSCAgICAgICAgPSBFdG1hYWxzb20gdmFuIGRlIG5lZXJzbGFnIChpbiAwLjEgbW0pICgtMSB2b29yIDwwLjA1IG1tKSAvIERhaWx5IHByZWNpcGl0YXRpb24gYW1vdW50IChpbiAwLjEgbW0pICgtMSBmb3IgPDAuMDUgbW0pDQpFVjI0ICAgICAgPSBSZWZlcmVudGllZ2V3YXN2ZXJkYW1waW5nIChNYWtraW5rKSAoaW4gMC4xIG1tKSAvIFBvdGVudGlhbCBldmFwb3RyYW5zcGlyYXRpb24gKE1ha2tpbmspIChpbiAwLjEgbW0pDQoNCg0KYGBge3J9DQptZXRlbyA9IHJlYWQudGFibGUoZmlsZSA9ICJwcHRfZXBvdDIwMjAudHh0IiwgaGVhZGVyID0gVFJVRSkNCm5hbWVzKG1ldGVvKSA9IGMoImRhdGUiLCJwcHQiLCJlcG90IikNCiNhZGp1c3RpbmcgdGhlIHVuaXQgdG8gY20vZGF5DQpwcmVjaXAgPSBtZXRlbyRwcHQvMTAwICANCmV2YXBvdCA9IG1ldGVvJGVwb3QvMTAwDQojY29ycmVjdGluZyBmb3IgdGhlIC0xIHN5bWJvbCBhbmQgc2V0dGluZyB0aGVzZSB2YWx1ZXMgdG8gMA0KcHJlY2lwW3doaWNoKHByZWNpcCA8IDApXSA9IDANCmV2YXBvdFt3aGljaChldmFwb3QgPCAwKV0gPSAwDQpgYGANCg0KDQojIyMgQmFzZSBydW4NCg0KYGBge3J9DQoNCg0KYmVnaW4udGltZSAgPSAxDQpjdXJyZW50LnRpbWUgPSBiZWdpbi50aW1lDQplbmQudGltZSA9IGxlbmd0aChtZXRlbyRkYXRlKQ0KZGVsdGEudCA9IDEgI2RheSwgaS5lLiB0aGUgdGltZSBzdGVwcGluZw0KDQoNCiMjIyNjb250YWluZXJzIHRvIHNhdmUgaW50ZXJtZWRpYXRlIHJlc3VsdHMNCnN0YXRlLnRtcCA9IGMoKQ0KdHJhbnNwLnRtcCA9IGMoKQ0KbW9pc3QudG1wID0gYygpDQpiYWxhbmNlLnRtcCA9IGMoKQ0KDQojI3NldCB0aGUgaW5pdGlhbCBzdGF0ZSB0byBlcXVpbGlicml1bQ0KZG8uaW5pdGlhbGl6ZShVbnNhdFRyYW5zLCBpbml0ID0gZHJuLnBzaS5lcXVpKQ0Kc3RhcnQucnVuLnRpbWUgPSBTeXMudGltZSgpDQp3aGlsZShjdXJyZW50LnRpbWUgPD0gZW5kLnRpbWUpDQp7DQogIHBwdCA9IHByZWNpcFthcy5pbnRlZ2VyKGN1cnJlbnQudGltZSldDQogIEVwb3QgPSBldmFwb3RbYXMuaW50ZWdlcihjdXJyZW50LnRpbWUpXQ0KICBwc2kub2xkID0gc3RhdGUuZnVuKFVuc2F0VHJhbnMpDQogIHNvbHZlLnN0ZXBzKFVuc2F0VHJhbnMsbWF4c3RhdGVjaGFuZ2Vjcml0ID0gMWUtMDQsdmVyYm9zZWxldmVsID0gMSkNCiAgcGxvdC53dGFibGUoVW5zYXRUcmFucykNCiAgY2F0KHBhc3RlKCJkYXkgbnIuOiIsY3VycmVudC50aW1lLCAiXG4iKSkNCiAgc2F2ZS5kYXRhKFVuc2F0VHJhbnMpDQogIGN1cnJlbnQudGltZSA9IGN1cnJlbnQudGltZSArIGRlbHRhLnQgDQogICMjIyNwbG90dGluZyBzb21lIHN0dWZmDQogICMgb2xkLnBhciA9IHBhcihuby5yZWFkb25seT1UUlVFKSAjdG8gc2F2ZSBjdXJyZW50IHBsb3Qgc2V0dGluZ3MgYW5kIHN0dWZmDQogICMgbGF5b3V0ID0gbWF0cml4KGMoMiwyKSxuY29sID0gMikjLCBieXJvdyA9IFRSVUUpDQogICMgcGxvdC5yb290cyhVbnNhdFRyYW5zKQ0KICAjIHBsb3QudGhldGEoVW5zYXRUcmFucykNCiAgIyBwbG90Lnd0YWJsZShVbnNhdFRyYW5zKQ0KICAjIHBhciA9IG9sZC5wYXINCiAgIyMjIyNzYXZpbmcgdGhlIHJlc3VsdHMNCn0NCmVuZC5ydW4udGltZSA9IFN5cy50aW1lKCkNCg0KY2F0KHBhc3RlKCJzaW11bGF0aW9uIHRpbWUgaXMgOiIsZW5kLnJ1bi50aW1lIC0gc3RhcnQucnVuLnRpbWUpKQ0KIyMgdG8gbWFrZSBkYXRhZnJhbWVzIGFuZCBjc3YncyBmcm9tIHRoZSBzYXZlZCBkYXRhDQpiYWxhbmNlLnJlcyA9IGRhdGEuZnJhbWUoYmFsYW5jZS50bXApDQp0aW1lLnNlcmllcyA9IHNlcShiZWdpbi50aW1lLGVuZC50aW1lLGJ5ID0gZGVsdGEudCkgI2FkZCB0aGlzIHRvIHRoZSBiYWxhbmMucmVzIGRhdGEuZnJhbWUgZm9yIHVzZSB3aXRoIElKUGxvdA0KYmFsYW5jZS5yZXMgPSBjYmluZCh0aW1lLnNlcmllcyxiYWxhbmNlLnJlcykNCm5hbWVzKGJhbGFuY2UucmVzKSA9IGMoIkRheW5yIiwiUHJlY2lwIiwiRXBvdCIsIkVhY3QiLCJTdG8yZmxvdyIsIkZsb3cyc3RvIiwiUmVjaGFyZ2UiLCJQb25kaW5nIikNCndyaXRlLmNzdih4ID0gYmFsYW5jZS5yZXMsIGZpbGUgPSAiQmFsYW5jZS5jc3YiKQ0Kc3RhdGUucmVzID0gZGF0YS5mcmFtZShzdGF0ZS50bXApDQpuYW1lcyhzdGF0ZS5yZXMpID0gYXMuY2hhcmFjdGVyKG5vZGVzKQ0Kd3JpdGUuY3N2KHggPSBzdGF0ZS5yZXMsIGZpbGUgPSAiU3RhdGVzLmNzdiIpDQp0cmFuc3AucmVzID0gZGF0YS5mcmFtZSh0cmFuc3AudG1wKQ0KbmFtZXModHJhbnNwLnJlcykgPSBhcy5jaGFyYWN0ZXIobm9kZXMpDQp3cml0ZS5jc3YoeCA9IHRyYW5zcC5yZXMsIGZpbGUgPSAiVHJhbnNwLmNzdiIpDQptb2lzdC5yZXMgPSBkYXRhLmZyYW1lKG1vaXN0LnRtcCkNCm5hbWVzKG1vaXN0LnJlcykgPSBhcy5jaGFyYWN0ZXIobm9kZXMpDQp3cml0ZS5jc3YoeCA9IG1vaXN0LnJlcywgZmlsZSA9ICJNb2lzdHVyZS5jc3YiKQ0KDQoNCiMjIGNyZWF0aW5nIHRoZSBiYXNlIG91dHB1dCBkYXRhLmZyYW1lcw0KYmFsYW5jZS5iYXNlID0gYmFsYW5jZS5yZXMNCnN0YXRlLmJhc2UgPSBzdGF0ZS5yZXMNCnRyYW5zcC5iYXNlID0gdHJhbnNwLnJlcw0KbW9pc3QuYmFzZSA9IG1vaXN0LnJlcw0KDQpgYGANCg0KDQoNCiMjIyBWaWV3aW5nIHRoZSByZXN1bHRzDQoNCkZpcnN0IGEgYmFzZSBwbG90IGZ1bmN0aW9uIHRvIGhhdmUgYSBsb29rIGF0IHRoZSB3YXRlciBiYWxhbmNlIHRlcm1zIGluIHRpbWUNCk5leHQgd2UgdXNlIE1JUExPVCBwYWNrYWdlIGZyb20gUGF1bCBUb3JmcyAoZm9ybWVyIGNvbGxlYWd1ZSkgdG8gbG9vayBhdCB0aGUgbW9pc3R1cmUsIHByZXNzdXJlIGFuZCByb290IHVwdGFrZSBkaXN0cmlidXRpb24gaW4gdGhlIG1vZGVsIGluIHRpbWUuIFRoZXNlIGFyZSB0aHVzIHRpbWUgc2VyaWVzDQoNCg0KDQpgYGB7cn0NCnBsb3QudHJhbnMuYmFsYW5jZSA9IGZ1bmN0aW9uKCkgIyB0aGlzIGlzIGEgdmVyeSBiYXNpYyBzaW1wbGUgYW5kIGdlbmVyaWMgUiBwbG90DQp7DQogIHRpbWUuc2VyaWVzID0gYygxOmxlbmd0aChiYWxhbmNlLnJlc1ssMV0pKQ0KICBmbHV4LnJhbmdlID0gcmFuZ2UoYmFsYW5jZS5yZXNbLDI6OF0pDQogIHBsb3QoeCA9IHRpbWUuc2VyaWVzICwgeSA9IGJhbGFuY2UucmVzJEVhY3QsIHlsaW0gPSBmbHV4LnJhbmdlLCB4bGFiID0gInRpbWUgKGQpIiwNCiAgICAgICB5bGFiID0gImZsdXggcmF0ZSBjbS9kIiwgbWFpbiA9ICJXYXRlciBiYWxhbmNlIHRlcm1zIiwgdHlwZSA9ICJsIiwNCiAgICAgICBsd2QgPSAyLCBjb2wgPSAicmVkIikNCiAgbGluZXMoeCA9IHRpbWUuc2VyaWVzLCB5ID0gYmFsYW5jZS5yZXMkU3RvMmZsb3csIGx3ZCA9IDIsY29sID0gImxpZ2h0Ymx1ZSIpDQogIGxpbmVzKHggPSB0aW1lLnNlcmllcywgeSA9IGJhbGFuY2UucmVzJEZsb3cyc3RvLCBsd2QgPSAyLCBjb2wgPSAibGlnaHRncmVlbiIpDQogIGxpbmVzKHggPSB0aW1lLnNlcmllcywgeSA9IGJhbGFuY2UucmVzJFJlY2hhcmdlLCBsd2QgPSAyLCBjb2wgPSAibWFnZW50YSIpDQogIGxpbmVzKHggPSB0aW1lLnNlcmllcywgeSA9IGJhbGFuY2UucmVzJFByZWNpcCwgbHdkID0gMiwgY29sID0gJ2JsdWUnKQ0KICBsaW5lcyh4ID0gdGltZS5zZXJpZXMsIHkgPSBiYWxhbmNlLnJlcyRQb25kaW5nLCBsd2QgPSAyLCBjb2wgPSAnb3JhbmdlJykNCiAgbGluZXMoeCA9IHRpbWUuc2VyaWVzLCB5ID0gYmFsYW5jZS5yZXMkRXBvdCwgbHdkID0gMiwgY29sID0gJ3RvbWF0bzMnKQ0KICBncmlkKCkNCn0NCnBsb3QudHJhbnMuYmFsYW5jZSgpICMgdGhpcyBhIHZlcnkgYmFzaWMgc2ltcGxlIGFuZCBnZW5lcmljIFIgcGxvdA0KDQojYmVsb3cgdGhlIElKUExPVCBwYWNrYWdlIGlzIHVzZWQgdG8gcGxvdCBkaWZmZXJlbnQgdHlwZSBvZiBkYXRhDQoNCmJhbGFuY2UucGxvdCA9IG5ld19taXBsb3QobmFtZSA9ICJiYWxhbmNlIixmeHJpZ2h0ID0gMC45MCwgcmlnaHRtYXJnaW4gPSAwLjApDQpjb2xvci5zZXQgPSByYWluYm93KG4gPSBsZW5ndGgoYmFsYW5jZS5yZXNbMSxdKSkNCmFkZF9saW5lKGJhbGFuY2UucGxvdCwgZGF0YSA9IGJhbGFuY2UucmVzLGRhdGF4bmFtZSA9IG5hbWVzKGJhbGFuY2UucmVzWzFdKSwgZGF0YXluYW1lID0gbmFtZXMoYmFsYW5jZS5yZXNbMl0pLCBjb2xvciA9IGNvbG9yLnNldFsyXSx3aWR0aCA9IDIpDQphZGRfbGluZShiYWxhbmNlLnBsb3QsIGRhdGEgPSBiYWxhbmNlLnJlcyxkYXRheG5hbWUgPSBuYW1lcyhiYWxhbmNlLnJlc1sxXSksIGRhdGF5bmFtZSA9IG5hbWVzKGJhbGFuY2UucmVzWzNdKSwgY29sb3IgPSBjb2xvci5zZXRbM10sd2lkdGggPSAyKQ0KYWRkX2xpbmUoYmFsYW5jZS5wbG90LCBkYXRhID0gYmFsYW5jZS5yZXMsZGF0YXhuYW1lID0gbmFtZXMoYmFsYW5jZS5yZXNbMV0pLGRhdGF5bmFtZSA9IG5hbWVzKGJhbGFuY2UucmVzWzRdKSwgY29sb3IgPSBjb2xvci5zZXRbNF0sd2lkdGggPSAyKQ0KYWRkX2xpbmUoYmFsYW5jZS5wbG90LCBkYXRhID0gYmFsYW5jZS5yZXMsZGF0YXhuYW1lID0gbmFtZXMoYmFsYW5jZS5yZXNbMV0pLGRhdGF5bmFtZSA9IG5hbWVzKGJhbGFuY2UucmVzWzVdKSwgY29sb3IgPSBjb2xvci5zZXRbNV0sd2lkdGggPSAyKQ0KYWRkX2xpbmUoYmFsYW5jZS5wbG90LCBkYXRhID0gYmFsYW5jZS5yZXMsZGF0YXhuYW1lID0gbmFtZXMoYmFsYW5jZS5yZXNbMV0pLGRhdGF5bmFtZSA9IG5hbWVzKGJhbGFuY2UucmVzWzZdKSwgY29sb3IgPSBjb2xvci5zZXRbNl0sd2lkdGggPSAyKQ0KYWRkX2xpbmUoYmFsYW5jZS5wbG90LCBkYXRhID0gYmFsYW5jZS5yZXMsZGF0YXhuYW1lID0gbmFtZXMoYmFsYW5jZS5yZXNbMV0pLGRhdGF5bmFtZSA9IG5hbWVzKGJhbGFuY2UucmVzWzddKSwgY29sb3IgPSBjb2xvci5zZXRbN10sd2lkdGggPSAyKQ0KYWRkX2xpbmUoYmFsYW5jZS5wbG90LCBkYXRhID0gYmFsYW5jZS5yZXMsZGF0YXhuYW1lID0gbmFtZXMoYmFsYW5jZS5yZXNbMV0pLGRhdGF5bmFtZSA9IG5hbWVzKGJhbGFuY2UucmVzWzhdKSwgY29sb3IgPSBjb2xvci5zZXRbOF0sd2lkdGggPSAyKQ0KYWRkX3hsYWJlbChiYWxhbmNlLnBsb3QsIHRleHQgPSAidGltZSBzdGVwcyAoZCkiLCBmb250c2l6ZSA9IDEyKQ0KYWRkX3lsYWJlbChiYWxhbmNlLnBsb3QsIHRleHQgPSAiZmx1eCByYXRlIChjbS9kKSIsIGZvbnRzaXplID0gMTIsIG9mZnNldCA9IDI1KQ0KYWRkX3RpdGxlKGJhbGFuY2UucGxvdCwgdGV4dCA9ICJXYXRlciBiYWxhbmNlIHRlcm1zIikNCmFkZF92bGluZShiYWxhbmNlLnBsb3QsIHggPSBzZXEobWluKGJhbGFuY2UucmVzWywxXSksbWF4KGJhbGFuY2UucmVzWywxXSksbGVuZ3RoLm91dCA9IDUpLCBkYXNoPWMoMiwyKSwgY29sb3IgPSAiZ3JleSIpDQphZGRfaGxpbmUoYmFsYW5jZS5wbG90LCB5ID0gc2VxKG1pbihiYWxhbmNlLnJlc1ssMjo4XSksbWF4KGJhbGFuY2UucmVzWywyOjhdKSxsZW5ndGgub3V0ID0gNSksZGFzaD1jKDIsMiksIGNvbG9yID0gImdyZXkiKQ0KDQpsZWdlbmQucGxvdCA9IG5ld19taXBsb3QobmFtZSA9ICJsZWdlbmQiLHlmaXhlZCA9IFRSVUUsIHhmaXhlZCA9IFRSVUUsZnhsZWZ0ID0gMC45MCxmeWJvdHRvbSA9IDAuNywNCiAgICAgICAgICAgICAgICAgICAgICAgICBib3R0b21tYXJnaW4gPSAwLCBsZWZ0bWFyZ2luID0gMCwgaG9yYXhpc251bWRpdiA9IDAsIHZlcmF4aXNudW1kaXYgPSAwLA0KICAgICAgICAgICAgICAgICAgICAgICAgIGltcG9zZWRtYXh4PSBtYXgobmNoYXIobmFtZXMoYmFsYW5jZS5yZXMpKSkpDQoNCmZvciAoaSBpbiAyOmxlbmd0aChiYWxhbmNlLnJlc1sxLF0pKQ0Kew0KICBhZGRfcG9pbnRzKG1pcCA9IGxlZ2VuZC5wbG90LHggPSAwLCB5ID0gbGVuZ3RoKGJhbGFuY2UucmVzWzEsXSktaSwgbWFya2VyID0gInNxdWFyZSIsc2l6ZSA9IDgsIGZpbGxjb2xvciA9IGNvbG9yLnNldFtpXSkNCiAgYWRkX3RleHQobWlwID0gbGVnZW5kLnBsb3QsIGF0eD0gMCwgYXR5ID0gbGVuZ3RoKGJhbGFuY2UucmVzWzEsXSktaSx0ZXh0ID0gcGFzdGUobmFtZXMoYmFsYW5jZS5yZXNbaV0pKSxmb250c2l6ZSA9IDEyLA0KICAgICAgICAgICBwb3MgPSAiRSIsb2Zmc2V0ID0gNSkNCn0NCg0KIyMjZm9yIHRlc3RpbmcNCndyaXRlLnRhYmxlKG1vaXN0LnJlcyxmaWxlID0gIm1vaXN0dXJlLmRhdCIpDQp3cml0ZS50YWJsZSh0cmFuc3AucmVzLCBmaWxlID0gInRyYW5zcC5kYXQiKQ0KDQojIyMNCmRyYXcoYmFsYW5jZS5wbG90LGxlZ2VuZC5wbG90KQ0KDQpmbHV4LnBsb3QgPSBuZXdfaWpwbG90KG5hbWUgPSAiZmx1eCIsIGZ4cmlnaHQgPSAwLjkwLCByaWdodG1hcmdpbiA9IDAuMCkNCm5yLnRpbWUuc3RlcHMgPSBsZW5ndGgobW9pc3QucmVzWywxXSkNCmFkZF9zZWxlY3RudW0oaWpwPWZsdXgucGxvdCxuYW1lID0gInRpbWUgc3RlcCIsbWF4bnVtID0gbnIudGltZS5zdGVwcykNCmZvciAoaSBpbiAxOm5yLnRpbWUuc3RlcHMpDQp7DQogIGFkZF9saW5lKGZsdXgucGxvdCx4ID0gbm9kZXMsIHkgPSBzb2lsdHlwZSR0aGV0YS5mdW4oc3RhdGUucmVzW2ksXSksIHNlbGVjdG51bSA9IGMoInRpbWUgc3RlcCIsIGkpLCBjb2xvciA9IGNvbG9yLnNldFsxXSkNCiAgYWRkX2xpbmUoZmx1eC5wbG90LHggPSBub2RlcywgeSA9IDEwKnRyYW5zcC5yZXNbaSxdLCBzZWxlY3RudW0gPSBjKCJ0aW1lIHN0ZXAiLCBpKSxjb2xvciA9IGNvbG9yLnNldFsyXSkNCn0NCmFkZF90aXRsZShmbHV4LnBsb3QsdGV4dCA9ICJmbHV4IHJhdGVzIHJvb3QgdXB0YWtlKG1tL2QpIG1vaXN0dXJlIGNvbnRlbnQoY20vY20pIikNCmFkZF94bGFiZWwoZmx1eC5wbG90LHRleHQgPSAidGltZSAoZCkiLGZvbnRzaXplID0gMTIpDQphZGRfeWxhYmVsKGZsdXgucGxvdCx0ZXh0ID0gImZsdXggcmF0ZSByb290IHVwdGFrZSAobW0vZCkgb3IgZnJhY3Rpb24gdGhldGEiLGZvbnRzaXplID0xMikNCmFkZF92bGluZShmbHV4LnBsb3QsIHg9IHNlcShtaW4obm9kZXMpLG1heChub2RlcyksbGVuZ3RoLm91dCA9IDUpLCBkYXNoID0gYygyLDIpLGNvbG9yID0gImdyZXkiKQ0KcHNpLm1heCA9IG1heChzdGF0ZS5yZXMpIyB0byBkZXRlcm1pbmUgdGhlIHJhbmdlIGZvciB0aGUgZ3JpZA0KcHNpLm1pbiA9IG1pbihzdGF0ZS5yZXMpDQp0cmFuc3AucmFuZ2UgPSByYW5nZSh0cmFuc3AucmVzKQ0KdGhldGEucmFuZ2UgPSBjKHNvaWx0eXBlJHRoZXRhLmZ1bihwc2kubWluKSxzb2lsdHlwZSR0aGV0YS5mdW4ocHNpLm1heCkpDQptYXgueS5heGlzID0gbWF4KHRoZXRhLnJhbmdlLCB0cmFuc3AucmFuZ2UpDQptaW4ueS5heGlzID0gbWluKHRoZXRhLnJhbmdlLCB0cmFuc3AucmFuZ2UpDQphZGRfaGxpbmUoZmx1eC5wbG90LCB5ID0gc2VxKG1pbi55LmF4aXMsbWF4LnkuYXhpcyxsZW5ndGgub3V0ID0gNSksIGRhc2g9YygyLDIpLGNvbG9yID0gImdyZXkiKQ0KDQpmbHV4LmxlZ2VuZCA9IG5ld19panBsb3QobmFtZSA9ICJmbHV4LmxlZ2VuZCIseWZpeGVkID0gVFJVRSwgeGZpeGVkID0gVFJVRSxmeGxlZnQgPSAwLjkwLGZ5Ym90dG9tID0gMC44LA0KICAgICAgICAgICAgICAgICAgICAgICAgIGJvdHRvbW1hcmdpbiA9IDAsIGxlZnRtYXJnaW4gPSAwLCByaWdodG1hcmdpbiA9IDIsIGhvcmF4aXNudW1kaXYgPSAwLCB2ZXJheGlzbnVtZGl2ID0gMCwNCiAgICAgICAgICAgICAgICAgICAgICAgICBpbXBvc2VkbWlueCA9IC0xLGltcG9zZWRtYXh4PSBuY2hhcigicm9vdCB1cHRha2UiKSxpbXBvc2VkbWlueSA9LTEsaW1wb3NlZG1heHkgPSAzKQ0KYWRkX3BvaW50cyhmbHV4LmxlZ2VuZCx4ID0gMC41LCB5ID0gMC41LCBtYXJrZXIgPSAidGlsZSIsIHNpemUgPSA4LCBmaWxsY29sb3IgPSBjb2xvci5zZXRbMV0pDQphZGRfcG9pbnRzKGZsdXgubGVnZW5kLHggPSAwLjUsIHkgPSAyLjUsIG1hcmtlciA9ICJ0aWxlIiwgc2l6ZSA9IDgsIGZpbGxjb2xvciA9IGNvbG9yLnNldFsyXSkNCmFkZF90ZXh0KGZsdXgubGVnZW5kLCBhdHggPSAwLjUsIGF0eSA9IDAuNSwgdGV4dCA9ICJtb2lzdHVyZSIsZm9udHNpemUgPSAxMiwgcG9zID0gIkUiLG9mZnNldCA9IDUpDQphZGRfdGV4dChmbHV4LmxlZ2VuZCwgYXR4ID0gMC41LCBhdHkgPTIuNSwgdGV4dCA9ICJyb290IHVwdGFrZSIsZm9udHNpemUgPSAxMiwgcG9zPSJFIixvZmZzZXQgPSA1KQ0KDQpkcmF3KGZsdXgucGxvdCxmbHV4LmxlZ2VuZCkNCg0KcHNpLnBsb3QgPSBuZXdfbWlwbG90KG5hbWUgPSAicHNpIiwgZnhyaWdodCA9IDAuOTAsIHJpZ2h0bWFyZ2luID0gMC4wKQ0KbnIudGltZS5zdGVwcyA9IGxlbmd0aChtb2lzdC5yZXNbLDFdKQ0KYWRkX3NlbGVjdG51bShtaXA9cHNpLnBsb3QsbmFtZSA9ICJ0aW1lIHN0ZXAiLG1heG51bSA9IG5yLnRpbWUuc3RlcHMpDQpmb3IgKGkgaW4gMTpuci50aW1lLnN0ZXBzKQ0Kew0KICBhZGRfbGluZShwc2kucGxvdCx4ID0gbm9kZXMsIHkgPSBzdGF0ZS5yZXNbaSxdLCBzZWxlY3RudW0gPSBjKCJ0aW1lIHN0ZXAiLCBpKSwgY29sb3IgPSBjb2xvci5zZXRbMV0pDQp9DQphZGRfdmxpbmUocHNpLnBsb3QseD1zZXEobWluKG5vZGVzKSxtYXgobm9kZXMpLGxlbmd0aC5vdXQgPSA1KSxkYXNoID0gYygyLDIpLCBjb2xvciA9ICdncmV5JykNCmFkZF9obGluZShwc2kucGxvdCx5ID0gc2VxKG1pbihzdGF0ZS5yZXMpLG1heChzdGF0ZS5yZXMpLCBsZW5ndGgub3V0ID0gNSksZGFzaCA9IGMoMiwyKSxjb2xvciA9ICJncmV5IikNCmFkZF90aXRsZShwc2kucGxvdCwgdGV4dCA9ICJwcmVzc3VyZSBkaXN0cmlidXRpb24gKHBzaSkgb3ZlciB0aW1lIikNCmFkZF94bGFiZWwocHNpLnBsb3QsIHRleHQgPSAidGltZSAoZCkiLGZvbnRzaXplID0gMTIpDQphZGRfeWxhYmVsKHBzaS5wbG90LCB0ZXh0ID0gInByZXNzdXJlIChjbSkiLGZvbnRzaXplID0gMTIpDQoNCmRyYXcocHNpLnBsb3QpDQoNCmBgYA0KDQo=